Directory of image this file is from
This file as a plain text file

C	SSP PAGE 320
	SUBROUTINE DET3(H,Y,Z,NDIM,IER)
	DIMENSION Y(1),Z(1)
	IF (NDIM-3)4,1,1
1	IF(H)2,5,2
2	HH=.5/H
	YY=Y(NDIM-2)
	B=Y(2)*Y(2)
	B=HH*(B+B-Y(3)-Y(1)-Y(1)-Y(1))
	DO 3 I=3,NDIM
	A=B
	B=HH*(Y(I)-Y(I-2))
3	Z(I-2)=A
	IER=0
	A=Y(NDIM-1)+Y(NDIM-1)
	Z(NDIM)=HH*(Y(NDIM)+Y(NDIM)+Y(NDIM)-A-A+YY)
	Z(NDIM-1)=B
	RETURN
4	IER=-1
	RETURN
5	IER =1
	RETURN
	END