plot rho*(-6378.1*sin(0.1)+sqrt(6378.1^2*(sin(0.1))^2+2*6378.1*7.746*1.3/rho+(7.746*1.3/rho)^2)) for rho from 0 to 1.3