Web pdp-10.trailing-edge.com

Trailing-Edge - PDP-10 Archives - decuslib20-01 - decus/20-0025/rkpbx.for
There is 1 other file named rkpbx.for in the archive. Click here to see a list.
```C	THE PROGRAM RKPBD.F4 IS A SAMPLE PROGRAM,WHEN USED IN
C	CONJUNCTION WITH RKPBX(LOADED AND SAVED AS RKPBD.SAV)
C	WILL PRODUCE A RUNNING EXAMPLE OF THE USE OF RKPBX.
SUBROUTINE RKPB2 (      TEMP,X,DX,Y,F,N)
DIMENSION TEMP(1),Y(1),F(1),RKCON(4)
401	M=N+1
RKCON(1)=1.0
RKCON(2)=2.0
RKCON(3)=2.0
RKCON(4)=1.0
X=TEMP(2*M+1)
DO 421 I=2,M
IP2M=I+2*M
Y(I-1)=TEMP(IP2M)
IP3M=IP2M+M
F(I-1)=TEMP(IP3M)
421	TEMP(I)=0.0
DO 491 J=1,4
DO 441 I=2,M
IPM=I+M
TEMP(IPM)=F(I-1)*DX
441	TEMP(I)=TEMP(I)+TEMP(IPM)*RKCON(J)
IF(J-4) 451,501,5001
451	X=TEMP(2*M+1)+DX/RKCON(J+1)
DO 461 I=2,M
IPM=I+M
IP2M=IPM+M
461	Y(I-1)=TEMP(IP2M)+TEMP(IPM)/RKCON(J+1)
CALL DERIV
491	CONTINUE
501	X=TEMP(2*M+1)+DX
DO 521 I=2,M
IP2M=I+2*M
521	Y(I-1)=TEMP(IP2M)+TEMP(I)/6.0
5001	RETURN
END
SUBROUTINE RKPB1 (      TEMP,X,DX,Y,F,N)
DIMENSION TEMP(1),Y(1),F(1)
101	M=N+1
TEMP(2*M+1)=X
DO 121 I=2,M
IP2M=I+2*M
121	TEMP(IP2M)=Y(I-1)
CALL DERIV
TEMP(3*M+1)=DX
DO 141 I=2,M
IP3M=I+3*M
141	TEMP(IP3M)=F(I-1)
5001	RETURN
END
CCCCCCCCCCCCCCCCCCCCCCCCCC

```