FUNCTION zvelocity,z

;given a redshift z, determines the velocity using the relativistic
;doppler formula
 
c = double(3.0e8) ;m/s

return, (-2.0*z*c-c*z^2)/(-2-2*z-z^2)

end