
function res=phi(v)

res=0.5*[cos(v(2)); sin(v(1))];