Convection is a process by which heat is transferred by movement of a heated water. The phenomenon where heat transfer process is driven by natural fluid motion due to temperature and density gradients is called free convection. The effect of slip velocity on thermal behaviours of MHD free convection flow over a vertical parallel two walls is the focus of this research. The mathematical equations which described the phenomenom includes momentum, magnetic field and energy equations. Condition for the existence and uniqueness of solution of the model equations were established using the approach of Lipschitz continuity. The model equations in dimensionless form were solved using Olayiwola’s generalized polynomial approximation method (OGPAM). The results obtained which contains the dimensionless parameters were presented graphically and discussed. It was observed that increase in Reynolds number and Knudsen number both enhanced the velocity distribution at the two walls. Also, magnetic field and temperature field were enhanced by magnetic prandtl number and Eckert number respectively.