  Different variants of the method of weighted residual finite element method
are used to get a solution for the parabolic heat equation, which is considered
to be the model equation for the steady state Navier-Stokes equations. Results
show that the Collocation and the Least-Squares variants are more suitable for
first order systems. Results also show that the Galerkin/Least-Squares method
is more diffusive than other methods, and hence gives stable solutions for a
wide range of P\'eclet numbers.
