10 ! ! COMPLEX SINH ! ! KEYED IN BY TIMOTHY SHANE MCCOY ! AT CENTRAL STATE UNIVERSITY ! ON PDP 11/45 ! VERSION RSTS/E V06C ! 20 DEF FNC(X) = (EXP(X) + EXP(-X))/2 30 DEF FNS(X) = (EXP(X) - EXP(-X))/2 40 INPUT X1,X2 50 PRINT 60 PRINT "INTERATION","REAL TERM","IMAG TERM","REAL SUM","IMAG SUM" 70 S1 = 0 75 A1,A2 = 0 80 S2 = 0 90 S0 = 1 100 FOR N = 1 TO 100 120 F = 1 130 P1 = 1 140 P2 = 0 150 FOR J = 1 TO 2*N-1 160 F = F * J 170 T1 = X1 * P1 - X2 * P2 180 T2 = X1 * P2 + X2 * P1 190 P1 = T1 200 P2 = T2 210 NEXT J 220 T1 = S0 * P1/F 230 T2 = S0 * P2/F 240 S1 = S1 + T1 250 S2 = S2 + T2 260 PRINT N,T1,T2,S1,S2 270 IF ABS(S1-A1) > 1E-6 THEN 290 280 IF ABS(S2-A2) < 1E-6 THEN 330 290 A1 = S1 300 A2 = S2 310 NEXT N 320 PRINT "DID NOT CONVERGE." 325 GOTO 340 330 PRINT "DID CONVERGE." 340 PRINT "SUM ="S1" +"S2"I" 350 PRINT "SINH ="EXP(X1) * COS(X2)/2 - EXP(-X1) * COS(-X2)/2" +"; \ PRINT EXP(X1) * SIN(X2)/2 - EXP(-X1) * SIN(-X2)/2"I" 360 PRINT 370 GOTO 10 400 END