Analysis of the mixed convection heat transfer in a uniformly heated horizontal pipe, first considered by Morton[Q. J. Mech. & Appl, Math.12(1959) 410], is extended with the aid of algebraic or symbolic computation packages. A set of four coupled, nonlinear partial differential equations that describe the effect of buoyancy-driven secondary flow on the pressure-driven axial flow are solved using a regular perturbation series expansion. A scale analysis suggests the product of Rayleigh and Reynolds numbers as the natural perturbation parameter and the series is constructed around the basic Poiseuille flow. In addition, the solutions also depend on the Prandtl number. The nonlinear problem is thus converted into an infinite series of linear partial differential equations to be solved sequentially. Using a further coordinate transformation into complex coordinates, solution of the infinite series of linear problems is reduced to one of evaluating quadratures, a task which is ideally suited for symbolic computation packages such as REDUCE, MAPLE, or MATHEMATICA. As an example the listing of a MAPLE program to solve the present problem to any desired order is appended. Although the symbolic computation program developed here is capable of solving the problem up to any desired order, the speed and memory limitations of a SUN workstation allowed us to generate analytical expressions for stream function, axial velocity and temperature only up to 10th order. The results provided by Morton for the first 3 terms are in agreement (except for some minor typographical errors) with our computer generated results. Primarily, we illustrate how the tedium of algebraic calculations can be relegated to algebraic computing machines. With the availability of faster CPU and larger memory machines, such problems can be routinely solved analytically to any desired accuracy.