plot(sqrt(.177-0.59*(x^0.8-0.61)^2)&&-sqrt(.177-0.59*(x^0.8-0.61)^2),x=0..1.3)