June 24 2002  11:02:16.398 AM
  
 ONEDFE
  
 Solve the two-point boundary value problem
  
   - d/dX (P dU/dX) + Q U  =  F
  
 on the interval [XL,XR], specifying
 the value of U or U' at each end.
  
 The interval [XL,XR] is broken into  5  subintervals
 Number of basis functions per element is  2
  
 The equation is to be solved for
 X greater than XL  =   0.E+0  and less than XR =  1.
  
 The boundary conditions are:
  
   At X = XL, U= 0.E+0
   At X = XR, U'= 1.
  
 Number of quadrature points per element is  1
  
   Node      Location
  
 0,  0.E+0
 1,  0.200000003
 2,  0.400000006
 3,  0.600000024
 4,  0.800000012
 5,  1.
  
 Subint    Length
  
 1,  0.200000003
 2,  0.200000003
 3,  0.200000018
 4,  0.199999988
 5,  0.199999988
  
 Subint    Quadrature point
  
 1,  0.100000001
 2,  0.300000012
 3,  0.5
 4,  0.700000048
 5,  0.899999976
  
 Subint  Left Node  Right Node
  
 1,  0,  1
 2,  1,  2
 3,  2,  3
 4,  3,  4
 5,  4,  5
  
   Node  Unknown
  
 0,  -1
 2*1
 2*2
 2*3
 2*4
 2*5
  
 Printout of tridiagonal linear system:
  
 Equation  ALEFT  ADIAG  ARITE  RHS
  
 1,  0.E+0,  10.,  -5.,  0.E+0
 2,  -5.,  10.,  -4.99999952,  0.E+0
 3,  -4.99999952,  10.,  -5.00000095,  0.E+0
 4,  -5.00000095,  10.0000019,  -5.00000095,  0.E+0
 5,  -5.00000095,  5.00000095,  0.E+0,  1.
 
Computed solution coefficients:
 
Node    X(I)        U(X(I))
 
     0  0.00   0.00000    
     1  0.20  0.200000    
     2  0.40  0.400000    
     3  0.60  0.600000    
     4  0.80  0.800000    
     5  1.00   1.00000