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