plot -6378.1*sin(theta)+sqrt(6378.1^2*(sin(theta))^2+2*6378.1*7.746+7.746^2) for theta from 0 to Pi/2