Trailing-Edge
-
PDP-10 Archives
-
decuslib10-01
-
43,50144/rkpbd.f4
There are no other files named rkpbd.f4 in the archive.
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