syms g v aA bA cA zA u t1 t2; + P$ e. d: N! D/ U( k5 L2 Fint(g); ; E Z* u N0 ^" K. Xint(v); : J. E2 ^3 D/ P6 |; Ah=0.856;6 X- L, Y2 k; N/ s/ h( { j4 L& o
aA=(1+tan(g)*tan(g));3 w- s u- u( Q7 {7 I
bA=-(2*1.625+2*tan(g)*(h+30*tan(g))-30*tan(g));+ q/ m6 S2 S& }9 D( r) Z7 _
cA=(h+30*tan(g)-15)*(h+30*tan(g)-15);; Y8 t% p& q# ?, b: c: z9 V1 D
zA=(-bA-sqrt(bA*bA-4*aA*cA))/(2*aA); " _ u7 x8 T5 [6 Y [' Rya=h+30*tan(g)-tan(g)*zA;& v2 R5 g. W, s' X1 ]$ n9 F
s1=(1.625*1.625-(z-1.625)*(z-1.625))*cos(v)*sqrt(tan(v)*tan(v)+1)*(sin(2*((asin(((h+30*tan(g)-z*tan(g))-15)/(sqrt(1.625*1.625-(z-1.625)*(z-1.625))*sqrt(tan(v)*tan(v)+1))))))/2-sin(2*((asin(-1/(cos(v)*sqrt(tan(v)*tan(v)+1))))))/2+((asin(((h+30*tan(g)-z*tan(g))-15)/(sqrt(1.625*1.625-(z-1.625)*(z-1.625))*sqrt(tan(v)*tan(v)+1)))))+((asin(-1/(cos(v)*sqrt(tan(v)*tan(v)+1))))));3 D% |- C' D1 r& w% c) g: Y# a
p=1.625-sqrt(1.625^2-((15-ya)/sqrt(tan(v)*tan(v)+1))); # f) J% N9 v6 L* ?* v+ c; d8 {H=int(s1,z,p,10);1 E% V1 s t, J; K- V