(*
Solution to Problem 4 of the Bigger Problem Set - Nonspherical thingys.
At the Mathematica prompt, type
<< spher.math
*)

r[k2_,ph_,th_] := r/. Solve[r^2 + k2 (3/2 Cos[th]^2 - 1/2) == ph r^3 , r][[3]];

ParametricPlot[{
	{r[.2,.5,th] Cos[th], r[.2,.5,th] Sin[th]},
	{r[.2,.7,th] Cos[th], r[.2,.7,th] Sin[th]},
	{r[.2,.7,th] Cos[th], r[.2,.7,th] Sin[th]},
	{r[.2,.8,th] Cos[th], r[.2,.8,th] Sin[th]},
	{r[.2,.9,th] Cos[th], r[.2,.9,th] Sin[th]},
	{r[.2,1,th] Cos[th], r[.2,1,th] Sin[th]},
	{r[.2,1.1,th] Cos[th], r[.2,1.1,th] Sin[th]},
	{r[.2,1.2,th] Cos[th], r[.2,1.2,th] Sin[th]},
	{r[.2,1.3,th] Cos[th], r[.2,1.3,th] Sin[th]},
	{r[.2,1.4,th] Cos[th], r[.2,1.4,th] Sin[th]},
	{r[.2,1.5,th] Cos[th], r[.2,1.5,th] Sin[th]},
	{r[.2,2,th] Cos[th], r[.2,2,th] Sin[th]}}
	, {th,0,Pi/2}]