Google
 

Trailing-Edge - PDP-10 Archives - decuslib20-01 - decus/20-0025/ampbd.for
There is 1 other file named ampbd.for in the archive. Click here to see a list.
	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