function phi = fcn(u,v,r,r_v_dot)

phi = [-u*v u*v 0 r_v_dot 0 0 0 abs(r)*r 0 0 0 r 0 0 0 0 0 0];