01 May 2006 11:25:46 AM

FEM2D_HEAT (C++ version):

  Compiled on May  1 2006 at 11:25:39.

  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 corner nodes of the triangles are generated by an
  underlying grid whose dimensions are

  NX =                 5
  NY =                 5

  Number of nodes    = 81
  Number of elements = 32

  The matrix half bandwidth is 18
  The matrix row size is       55

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

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

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

FEM2D_HEAT:
  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.05    0.00423576      0.128476
           0.1    0.00401111      0.122263
          0.15    0.00380908      0.116322
           0.2    0.00362096       0.11066
          0.25    0.00344343      0.105269
           0.3    0.00327508      0.100137
          0.35    0.00311516     0.0952552
           0.4    0.00296314     0.0906103
          0.45    0.00281859     0.0861916
           0.5     0.0026811     0.0819881

FEM2D_HEAT:
  Normal end of execution.

01 May 2006 11:25:47 AM