ADAS Subroutine nspij3

       SUBROUTINE NSPIJ3(N,H,W)
       IMPLICIT REAL*8(A-H,O-Z)
C-----------------------------------------------------------------------
C  PURPOSE: CALCULATE SPLINES WITH VARIOUS END CONDITIONS.
C
C  USES LABELLED COMMON  /SPL3/
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 APPROPRIAT
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  *********                                           *************
C-----------------------------------------------------------------------
C IDL-UNIX CONVERSION:
C
C  VERSION: 1.1                                        DATE: 01/10/96
C  MODIFIED:  WILLIAM OSBORN
C           - FIRST WRITTEN. NO CHANGES.
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-----------------------------------------------------------------------
C-----------------------------------------------------------------------
      INTEGER             N
      REAL*8              H(10),       W(10,10)