(*Plots an ellipse*)

h = 1000. ;
b = 100. ;
v = 1500. ;
vt0 = v*t0 ;

hn1 = Sqrt[(h^2*(-b^2 + h^2 - b*vt0)*(-b^2 + h^2 + b*vt0))/
  (b^2 + b*h - h^2)^2] ;

Plot[hn1,{t0,0,2}]