Explore how the results vary when parameter Po varies. Create an array of 100 Po values covering the range from Po = 1 to Po = 100. For  each value of Po in the array repeat the calculation above and find the corresponding values of R, and M. You thus should have 1000 values of R, and  M. Plot M, as a function of R, in another plot window. Using the worked example: exoplanets" in Section 08, make a power-law fit to the M, vs R,relationship, and display it on the same plot.
int: making some plots logarithmic may be a good idea.