C CALCULATE SUMS OF SQUARES AND PRODUCTS 



at" r »0» 



SFX=0. 



SFY=0. 



SFQ=o. 



SFZ=0, 



DT=, 31 *15926536 



A«-OT 



DO 15 I«1,*0 



A«A*DT 



XII)=COS (A) 



e=i.*A 



Y< I >»Cr)S<B» 



G( I):>SINtB> 



CxSlN< A) 



Z<I>»C*ABS (C) 



F( I >aFP( I )-SF 



SFF»SFF*F( I )*F(I J 



SFX«SFX*F(I |*X( I » 



SFY=SFV*F( I >*Y< I ) 



SFQ=SFQ«-F( I )*Q( I ) 



SFZsSFZ*F(I )*2( I) 

 IS CONTINUE 



AX«SFX/20. 



AYaSFy/20» 



AQ»SFQ/20. 



AZ«SF2/15, 



VX=AX*SFX 



VY=AY«3FY 



VQ=Aa«SFQ 



VZxAZ^SFZ 



VB »SFF-VX-VY-VZ-VQ 

 C CALCULATE COEFFICIENTS AND PARAMETERS PHI AND K 



PHl328.64789*ATAN2(AQ,AY) 



IF(PHI,LT,-45. )7999,a999 

 7999 PMIxPHI+lPO» 

 8999 CONTINUE 



YAsSQRT(»a*AQ+AY*Ay) 



CLV=2,«YA/FLV 



ANG3( ANGLE*PI l/lSO* 



CLVAxCLV/COS< ANG) 



CLVU=CLVA/COS( ANG) 



CMV3-AX/FMV 



CDV=-AZ/FOV 



XK=.5-<SF/( CLV*FLV) ) 

 C PRINT RESULTS OF ANALYSIS 



PRINT 200 

 200 FCRMATi IHl ) 



PRINT A 

 ♦ FORMAT ( 1 OX, 6HT (SEC) , 7X , 6HHT (F T ) , 5X, 9HWA VEL ( FT ) , 6X , 7M0EP(F T > ,4X,9HU 

 IMAX(FPS) ,4X,9HCLEAR(FT) ,6X,7HDIA<FT),3X,I2HCYL LGTMIFT ) t4X, BHANG (O 

 2EG) ) 



PRINT 300, T (H,XL,D,U,CL,OI A,XC, ANGLE 

 300 F0RMAT(2X,F 13. 2,F13.3,F13.2,5F13, J,F13. 1////) 



PRINT 307 

 307 FaRMAT(6X,12HT0TAL SUM SO t6X , SHC aS2 A , 1 OX, 5 HS I N2 A , 1 1 X , 4HC0S A , SX , 1 OH 



142 



