Google
 

Trailing-Edge - PDP-10 Archives - decuslib20-01 - decus/20-0025/rkpbd.for
There is 1 other file named rkpbd.for in the archive. Click here to see a list.
	COMMON U(3),F(3)
	DIMENSION TEMP(16)
	T=0.
	DT=0.0625
	N=3
	U(1)=0.
	U(2)=2.0
	U(3)=1.0
	TF=2.0
	TYPE 1
1	FORMAT(/8X,1HT,13X,1HX,13X,1HY,13X,1HZ)
10	CALL RKPB1(TEMP,T,DT,U,F,N)
	TYPE 15,T,U
15	FORMAT(1P4E14.6)
	IF(T-TF) 20,100,100
20	CALL RKPB2(TEMP,T,DT,U,F,N)
	GO TO 10
100	STOP
	END
C	DERIVATIVE  EVALUATION SUBROUTINE
	SUBROUTINE DERIV
	COMMON U(3),F(3)
	F(1)=U(2)
	F(2)=-4.*U(1)
	F(3)=2.*U(3)
	RETURN
	END