Google
 

Trailing-Edge - PDP-10 Archives - decuslib20-01 - decus/20-0025/polyv.for
There is 1 other file named polyv.for in the archive. Click here to see a list.
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