ADAS Subroutine b2spij3

       SUBROUTINE B2SPIJ3(N,H,W)                                        
C
       IMPLICIT REAL*8(A-H,O-Z)                                         
C-----------------------------------------------------------------------
C
C  ****************** FORTRAN77 SUBROUTINE: B2SPIJ3 ********************
C
C PURPOSE:
C  SUBROUTINE TO CALCULATE SPLINES WITH VARIOUS END CONDITIONS.
C
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 
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 NOTES: THIS ROUTINE IS NOT YET PROPERLY ANNOTATED
C
C UNIX-IDL PORT:
C
C VERSION: 1.1                          DATE: 06-03-96
C MODIFIED: TIM HAMMOND (TESSELLA SUPPORT SERVICES PLC)
C               - PUT UNDER S.C.C.S. CONTROL
C
C-----------------------------------------------------------------------
C
      INTEGER             N
      REAL*8              H(10),       W(10,10)