17 September 2021   9:31:12.965 AM
 
FEM1D
  FORTRAN90 version.
 
  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 NSUB =        5 subintervals
  Number of basis functions per element is NL =        2
 
The equation is to be solved for
X greater than XL =    0.00000    
 and less than XR =    1.00000    
 
The boundary conditions are:
 
  At X = XL, U=   0.00000    
  At X = XR, U'=   1.00000    
 
Number of quadrature points per element is        1
 
  Node      Location
 
         0     0.00000    
         1    0.200000    
         2    0.400000    
         3    0.600000    
         4    0.800000    
         5     1.00000    
 
Subint    Length
 
         1    0.200000    
         2    0.200000    
         3    0.200000    
         4    0.200000    
         5    0.200000    
 
Subint    Quadrature point
 
         1    0.100000    
         2    0.300000    
         3    0.500000    
         4    0.700000    
         5    0.900000    
 
Subint  Left Node  Right Node
 
         1         0         1
         2         1         2
         3         2         3
         4         3         4
         5         4         5
 
  Number of unknowns NU =        5
 
  Node  Unknown
 
         0        -1
         1         1
         2         2
         3         3
         4         4
         5         5
 
Printout of tridiagonal linear system:
 
Equation  ALEFT  ADIAG  ARITE  RHS
 
         1   0.00000         10.0000        -5.00000         0.00000    
         2  -5.00000         10.0000        -5.00000         0.00000    
         3  -5.00000         10.0000        -5.00000         0.00000    
         4  -5.00000         10.0000        -5.00000         0.00000    
         5  -5.00000         5.00000         0.00000         1.00000    
 
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    
 
FEM1D:
  Normal end of execution.
 
17 September 2021   9:31:12.966 AM