19 January 2020 11:30:35 AM

FEM1D
  C 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.000000
  and less than XR = 1.000000

  The boundary conditions are:

  At X = XL, U = 0.000000
  At X = XR, U' = 1.000000

  Number of quadrature points per element is 1

  Node      Location

         0        0.000000 
         1        0.200000 
         2        0.400000 
         3        0.600000 
         4        0.800000 
         5        1.000000 

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.000000       10.000000       -5.000000        0.000000
         2       -5.000000       10.000000       -5.000000        0.000000
         3       -5.000000       10.000000       -5.000000        0.000000
         4       -5.000000       10.000000       -5.000000        0.000000
         5       -5.000000        5.000000        0.000000        1.000000

  Computed solution coefficients:

  Node    X(I)        U(X(I))

         0  0.000000        0.000000
         1  0.200000        0.200000
         2  0.400000        0.400000
         3  0.600000        0.600000
         4  0.800000        0.800000
         5  1.000000        1.000000

FEM1D:
  Normal end of execution.

19 January 2020 11:30:35 AM