function y=gdes(x)
y = x;

y(1) = cos(x(2))*sin(x(1));
y(2) = sin(x(2))*cos(x(1));