Google
 

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