The aim of this paper is to analyze the conjugate problems of heat conduction in solid walls coupled with laminar free convection flow adjacent to a vertical flat plate under boundary layer approximation. Using the similarity transformations the governing boundary layer equations for momentum and energy are reduced to a system of partial differential equations and then solved numerically using Finite Difference Method(FDM) known as the Keller-box scheme. Computed solutions to the governing equations are obtained for a wide range of non-dimensional parameters that are present in this problem, namely the coupling parameter P. the Prandtl number Pr and the heat generation parameter Q. The variations of the local heat transfer rate as well as the interface temperature and the friction along the plate and typical velocity and temperature profiles in the boundary layer are shown graphically. Numerical solutions have been consider for the Prandtl number Pr=0.70