Nonlinear heat conduction problems are solved by the Galerkin method with trial functions which are piecewise polynomials. The resulting system of nonlinear ordinary differential equations is solved by linear one-step and two-step A-stable methods. The final schemes are algebraically linear and unconditionally stable.