Trailing-Edge
-
PDP-10 Archives
-
decuslib10-01
-
43,50144/ampbd.f4
There are no other files named ampbd.f4 in the archive.
COMMON U(3),F(3)
DIMENSION TEMP(40)
T=0.
DT=0.0625
N=3
U(1)=0.
U(2)=2.0
U(3)=1.
TF=2.
IND=0
TYPE 1
1 FORMAT(/8X,1HT,13X,1HX,13X,1HY,13X,1HZ)
NITER=0
MTST=0
10 CALL AMPB1(IND,TEMP,T,DT,U,F,N,ICOUNT,NITER,MTST)
TYPE 15,T,U
15 FORMAT(1P4E14.6)
IF(T-TF) 20,100,100
20 CALL AMPB2(IND,TEMP,T,DT,U,F,N,ICOUNT,NITER,MTST)
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