PROGRAM CALCULATE INTEGER PH (2),PHI (2),PP (2),BET (2),BT (2),VELC (2),VELB (2),ACCC (2),ACCB (2),C1 (2),C2 (2),C3 (2),C4 (2),B1 (2),B2 (2),B3 (2),B4 (2) READ (*,*) A,B,C,D, VELA, ACCA, THETA PI=4.0 * ATAN (1.0) THET = 0 IHT = 180/THETA DTHET= PI/IHT DO

PROGRAM CALCULATE
INTEGER PH (2),PHI (2),PP (2),BET (2),BT (2),VELC (2),VELB (2),ACCC (2),ACCB (2),C1 (2),C2 (2),C3 (2),C4 (2),B1 (2),B2 (2),B3 (2),B4 (2)
READ (*,*) A,B,C,D, VELA, ACCA, THETA
PI=4.0 * ATAN (1.0)
THET = 0
IHT = 180/THETA
DTHET= PI/IHT
DO 10 J = 1, 2 * IHT
THET = (J - 1) * DTHET
AK =( (A * A - B * B + C * C + D * D) * 0.5)
TH = THET * 180/PI
AA = AK - A * (D - C) * COS (THET) - (C * D)
BB = -2.0 * A * C * SIN (THET)
CC = AK - A * (D + C) * COS (THET) + (C *D)
AB = BB * * 2 - 4 * AA * CC
IF (AB . LT . 0) GOTO 10
PHH = SQRT (AB)
PH (1) = - BB + PHH
PH (2) = - BB - PHH
DO 9 I =    1 ,2
PHI (1) = ATAN (PH (I) * 0.5/AA) * 2
PP (I) = PHI (I) * 180/PI
BET (I) = ASIN ((C * SIN (PHI (I) - A * SIN ( THET)) / B)
BT (I) = BET (I) * 180/PI
VELC (I) = A * VELA * SIN (BET (I) - THET) / (C * SIN (BET (I) - PHI (I)))
VELB (I) = (A * VELA * SIN (PHI (I) - THET) ) / (B * SIN (BET (I) - PHI (I))))
C1 (I) = A * ACCA * SIN (BET (I) - THET)
C2 (I) = A * VELA * * 2 * COS (BET (I) - THET) + B * VELB (I) * * 2
C3 (I) = C * VELC (I) * * 2 * COS (PHI (I) - BET (I) )
C4 (I) = C * SIN (BET (I) - PHI (I))
ACCC (I) = (C1 (I) - C2 (I) + C3 (I) / C4 (I)
B1 (I) = A * ACCA * SIN (PHI (I) -THET )
B2 (I) = A * VELA * * 2 * COS (PHI (I) - THET )
B3 (I) = B * VELB (I) * * 2 * COS (PHI (I) - BET (I) ) - C * VELC (I) * * 2
B4 (I) = B * (SIN (BET (I) - PHI (I))))
ACCB (I) = (B1 (I) - B2 (I) - B3 (I)) / B4 (I)
IF (J .NE .1) GO TO 8
WRITE (*, 7)
FORMAT (4X, ’  THET’ ,4X, ’ PHI’ , 4X, ’ BETA’, 4X, ’  VELC’ , 4X, ’  VELB’ , 4X , ’ ACCC , 4X,’ACCB’)
WRITE (*, 6) TH, PP (1), BT (1), VELC (1), VELB (1), ACCC (1), ACCB (1)
FORMAT (8F8 . 2)
WRITE (*, 5) PP (2), BT (2), VELC (2), VELB (2), ACCC (2), ACCB (2)
FORMAT (8X, 8F8 . 2)
CONTINUE
STOP 
END PROGRAM CALCULATE