function [y]=func(x)
   y=x;
   y(3)=x(3)-sin(1.5*(x(3)+x(4)));
   y(4)=x(4)-cos(1.5*(x(3)+x(4)));
endfunction