05 March 2020 08:05:15 AM FEM2D_HEAT_RECTANGLE C++ version Compiled on Mar 5 2020 at 07:59:32. 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_RECTANGLE: Wrote an EPS file "rectangle_nodes.eps". containing a picture of the nodes. FEM2D_HEAT_RECTANGLE: Wrote an ASCII node file rectangle_nodes.txt of the form X(I), Y(I) which can be used for plotting. FEM2D_HEAT_RECTANGLE: Wrote an EPS file "rectangle_elements.eps". containing a picture of the elements. FEM2D_HEAT_RECTANGLE: Wrote an ASCII element file "rectangle_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_RECTANGLE: Normal end of execution. 05 March 2020 08:05:15 AM