function [tetha2]=Fcn1(u)
d=sqrt(u(1,:)^2+u(2,:)^2)
a=atan2((u(3,:)),d)*(180/pi)
e=sqrt(u(1,:)^2+u(2,:)^2+u(3,:)^2)
f=40*e
c2=((u(1,:)^2+u(2,:)^2+u(3,:)^2)+231)/f
s2=sqrt(1-c2^2)
b=atan2(s2,c2)*(180/pi)
tetha2=-(a+b)