function z=flv1(t,v)
g=9.81; 
z=g-v.^2;
