ADAS Subroutine lspij3

       SUBROUTINE LSPIJ3(N,H,W)
       IMPLICIT REAL*8(A-H,O-Z)
C-----------------------------------------------------------------------
C                                                                       
C  ************  FORTRAN 77 SUBROUTINE: LSPIJ3 *************************
C                                                                       
C  PURPOSE: CALCULATE SPLINES WITH VARIOUS END CONDITIONS.
C
C  EXTENDED ARRAY DIMENSION VERSION OF NSPIJ3
C
C  USES LABELLED COMMON  /LSPL3/
C
C  CONDITIONS AT 1ST NODE AND NTH NODE CONTROLLED BY IEND1 AND IENDN
C      IEND=1 : SPECIFIED D LOG(Y) IE. DY/Y AT NODE STORED IN APPROPRIAT5070000
C               APPROPRIATE VECTOR
C          =2 : ZERO CURVATURE
C          =3 : CONSTANT CURVATURE
C          =4 : MATCHED TO SPECIFIED FUNCTIONAL FORM IN TERMS OF
C               TWO PARAMETERS A AND B SUCH THAT
C                        FUNCT = P(1)*A+Q(1)*B
C                   1ST DERIV. = P(2)*A+Q(2)*B
C                   2ND DERIV. = P(3)*A+Q(3)*B
C               WHERE A1,B1,P1,Q1 ARE USED FOR 1ST NODE AND
C               AN,BN,PN,QN FOR NTH NODE
C
C  INPUT
C      N=NUMBER OF KNOTS
C      H(I)=INTERVALS BETWEEN KNOTS
C  OUTPUT
C      W=SPLINE MATRIX
C
C AUTHOR:
C
C  ******* H.P. SUMMERS, JET           7 FEB 1989    ***************
C
C UNIX-IDL CONVERSION:
C
C VERSION: 1.1                          DATE: 07-10-96
C MODIFIED: WILLIAM OSBORN
C               - FIRST CONVERTED.
C
C VERSION: 1.2                          DATE: 15-05-07
C MODIFIED: Allan Whiteford
C               - Updated comments as part of subroutine
C                 documentation production.
C
C----------------------------------------------------------------------
      INTEGER             N
      REAL*8              H(40),       W(40,40)