Trailing-Edge
-
PDP-10 Archives
-
decuslib10-01
-
43,50144/polyv.f4
There are no other files named polyv.f4 in the archive.
00100 SUBROUTINE POLYV(N,A,RHO,PHI,M,R,C)
00110 DIMENSION A(1)
00120 V1=1.0
00130 V2 = 0.0
00135 C = 0.
00140 R=A(1)
00150 IF (M) 10, 10, 20
00160 10 W1=RHO*COS(PHI)
00170 W2=RHO*SIN(PHI)
00180 GO TO 30
00190 20 W1=RHO
00200 W2=PHI
00210 30 NN=N+1
00220 DO 40 I=2,NN
00230 T1=W1*V1-W2*V2
00240 V2=W2*V1+W1*V2
00250 V1=T1
00260 R=R+A(I)*V1
00270 C=C+A(I)*V2
00280 40 CONTINUE
00290 RETURN
00300 END