clc;clear;
syms u1 t1 v1 bt2 alf th1 th2
%syms th1
dt=pi/180;
v1=3*dt;
i=16/15;
bet1=21.5*dt;
u1=70;
for alf=2.5*dt;
th1=0:pi/60:2*pi;
st=-th1;
t1=atan(-i*sin(alf)*sin(st)/(i*(sin(alf)*cos(st)*sin(bet1)+cos(alf)*cos(bet1))-cos(bet1)));
x=-u1.*sin(v1).*(cos(t1).*((4190228212462809*sin((16*th1)/15).*sin(alf))/4503599627370496+(1650574788106375.*cos((16*th1)/15).*sin(th1))/4503599627370496-(1650574788106375*sin((16*th1)/15).*cos(alf).*cos(th1))/4503599627370496)-sin(t1).*(cos((16*th1)/15).*cos(th1)+sin((16*th1)/15).*cos(alf).*sin(th1)))-u1.*cos(v1).*((1650574788106375.*sin((16*th1)/15).*sin(alf))/4503599627370496-(4190228212462809*cos((16*th1)/15).*sin(th1))/4503599627370496+(4190228212462809*sin((16*th1)/15).*cos(alf).*cos(th1))/4503599627370496);
y=u1.*sin(v1).*(cos(t1).*((1650574788106375*sin((16*th1)/15).*sin(th1))/4503599627370496-(4190228212462809*cos((16*th1)/15).*sin(alf))/4503599627370496+(1650574788106375*cos((16*th1)/15).*cos(alf).*cos(th1))/4503599627370496)-sin(t1).*(sin((16*th1)/15).*cos(th1)-cos((16*th1)/15).*cos(alf).*sin(th1)))-u1.*cos(v1).*((1650574788106375*cos((16*th1)/15).*sin(alf))/4503599627370496+(4190228212462809*sin((16*th1)/15).*sin(th1))/4503599627370496&#