>> fem1d

10-Nov-2006 07:43:29

FEM1D
  MATLAB version

  Solve the two-point boundary value problem:

    -d/dx (p(x) du/dx) + q(x)*u  =  f(x)

  on an interval [xl,xr], with the values of
  u or u' specified at xl and xr.

  The interval is broken into 5 subintervals.
  The number of basis functions per element is 2

The equation is to be solved for
X greater than XL  =   and less than XR = 
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

  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                   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                    1.000000

Computed solution:

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.

10-Nov-2006 07:43:30
>>