function [tetha3]=Fcn2(u) s3=((u(1,:)^2+u(2,:)^2+u(3,:)^2)-569)/520 c3=sqrt(1-s3^2) tetha3=-atan2(s3,c3)*(180/pi)