The Burnett equations with slip boundary conditions are used to simulate the compressible gas flow and heat transfer in micro Poiseuille flow in the slip and transition flow regime. A relaxation method on Burnett terms is proposed in the present study and the thermal creep effect is considered. Convergent results at Knudsen number up to 0.4 are achieved and the results agree very well with experimental data. It is found that with the increase of Knudsen number, the Poiseuille number decreases while Nusselt number increases. The local Poiseuille number decreases along the whole channel while the local Nusselt number decreases rapidly first and then increases slowly afterwards. (C) 2008 Elsevier Ltd. All rights reserved.