the boundary;
        </li>
        <li>
          <b>FUNCTION H = INITIAL_CONDITION ( NODE_NUM, NODE_XY, TIME )</b>
          evaluates H(x,y,t), and is only called for TIME = TINIT.
        </li>
      </ul>
    </p>

    <p>
      The program writes out a file containing the
      solution value at every node.  This file may be used to create
      contour plots of the solution.
    </p>

    <h3 align = "center">
      Program Usage:
    </h3>

    <p>
      <blockquote>
        <b>fem2d_heat</b> ( <i>'