07-Jan-2022 20:00:50

fem2d_heat_rectangle_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Test fem2d_heat_rectangle().

FEM2D_HEAT_RECTANGLE
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2

  Solution of the time dependent heat equation on a unit
  box in 2 dimensions.

  Ut - Uxx - Uyy = F(x,y,t) in the box
        U(x,y,t) = G(x,y,t) for (x,y) on the boundary.
        U(x,y,t) = H(x,y,t) for t = T_INIT.

  The finite element method is used, with piecewise
  quadratic basis functions on 6 node triangular
  elements.

  The backward Euler formula is used for the time derivative.

  The corner nodes of the triangles are generated by an
  underlying grid whose dimensions are

  NX =                         7
  NY =                         7

  Number of nodes =            169
  Number of elements =         72
  The matrix half bandwidth is 26
  The matrix row size is       79

FEM2D_HEAT_RECTANGLE:
  Wrote an EPS node file
    "nodes.eps".
  containing a picture of the nodes.

FEM2D_HEAT_RECTANGLE:
  Wrote an ASCII node file
    "nodes.txt"
  of the form
    X(I), Y(I)
  which can be used for plotting.

FEM2D_HEAT_RECTANGLE:
  Wrote an EPS triangulation file
    "elements.eps".
  containing a picture of the triangulation.

FEM2D_HEAT_RECTANGLE:
  Wrote an ASCII element file
    "elements.txt"
  of the form
    Node(1) Node(2) Node(3) Node(4) Node(5) Node(6)
  which can be used for plotting.

     Time        L2 Error       H1 Error

      0.000000      0.000000      0.000000
      0.050000      0.023838      0.121133
      0.100000      0.046462      0.214592
      0.150000      0.067813      0.306985
      0.200000      0.087952      0.395284
      0.250000      0.106941      0.479017
      0.300000      0.124838      0.558184
      0.350000      0.141700      0.632917
      0.400000      0.157579      0.703389
      0.450000      0.172527      0.769789
      0.500000      0.186589      0.832303

fem2d_heat_rectangle_test():
  Normal end of execution.

07-Jan-2022 20:00:52