function [x3]=z(u)
x3=-13*cos(u(2,:)+u(3,:))-20*sin(u(2,:))