Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-01 - 43,50144/polyc.f4
There are no other files named polyc.f4 in the archive.
00100	      SUBROUTINE POLYC(RR,CR,NARG,A)
00110		DIMENSION RR(100),CR(100),A(3),B(3),Q(26)
00120	      J=1
00130	      IF (CR(J))  10, 20, 10
00140	10    NA=3
00150	      A(1)=RR(J)**2+CR(J)**2
00160	      A(2)=-2.*RR(J)
00170	      A(3)=1.
00180	      J=J+2
00190	      GO TO 29
00200	20    NA=2
00210	      A(1)=-RR(J)
00220	      A(2)=1.
00230	      J=J+1
00240	29    IF (J-NARG)  30, 30, 1000
00250	30    IF (CR(J))  40, 60, 40
00260	40    NB=3
00270	      B(1)=RR(J)**2+CR(J)**2
00280	      B(2)=-2.*RR(J)
00290	      B(3)=1.
00300	      J=J+2
00310	      GO TO 80
00320	60    NB=2
00330	      B(1)=-RR(J)
00340	      B(2)=1.
00350	      J=J+1
00360	80    NC=NA+NB-1
00370	      DO 120 I=1,NC
00380	      TEMP=0.0
00390	      DO 110 JP=1,I
00400	      IF (JP-NA)  90, 90, 110
00410	90    N2=I-JP+1
00420	      IF (N2-NB)  100, 100, 110
00430	100   TEMP=TEMP+A(JP)*B(N2)
00440	110   CONTINUE
00450	      Q(I)=TEMP
00460	120   CONTINUE
00470	      NA=NA+NB-2
00480	      NAPL1=NA+1
00490	      DO 130 I=1,NAPL1
00500	130   A(I)=Q(I)
00510	      IF (NA-NARG)  140, 1000, 1000
00520	140   NA=NA+1
00530	      GO TO 30
00540	1000  END