function [x2]=y(u)
x2=sin(u(1,:))*(-13*sin(u(2,:)+u(3,:))+20*cos(u(2,:)))