File DET3.FT (FORTRAN source file)
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
Feel free to contact me, David Gesswein djg@pdp8online.com
with any questions, comments on the web site, or if you have related equipment,
documentation, software etc. you are willing to part with. I am
interested in anything PDP-8
related, computers, peripherals used with them, DEC or third party, or
documentation.
PDP-8 Home Page
 PDP-8 Site Map
 PDP-8 Site Search