polar plot (-sin(theta) + sqrt(1+63*cos^2(theta)))/(2*cos^2(theta)) from 0 to pi