t,h,g,x,a=var('t,h,g,x,a')
tn2(t,h,g)=g^2*t^2 - h^2*(1-g^2)
zo(t,h,g,x)=sqrt(tn2(t,h,g)+g^2*x^2/(1-g^2))
zop(t,h,g,a)=sqrt(tn2(t,h,g)+a^2*(1-g^2)*g^2)
xop(a,g)=a*(1-g^2)
def Push(t,g):
    return parametric_plot([xop(a,g),zop(t,1,g,a)],(a,-3.5*t,3.5*t))
p = Push(0.8,0.8)+Push(1.2,0.8)+Push(1.6,0.8)+Push(2.,0.8)+Push(2.4,0.8)+Push(2.8,0.8)+Push(3.2,0.8)+Push(3.6,0.8)
p.axes_range(xmin=-4.75,xmax=4.75,ymin=7,ymax=0)
p.save(frame=True,axes=False,axes_labels=["Midpoint (km)","Pseudo-depth (km)"],fontsize=12,aspect_ratio=9.5/(7*golden_ratio),filename='junk_sage.pdf')