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