The program computes the velocity and temperature profiles for free convection past an isothermal vertical plate. Results are obtained after solving two coupled nonlinear differential equations. The shooting method is used to obtain the boundary conditions, required to solve the coupled ODEs, which match the conditions in the outer fluid. We plot velocity and temperature profiles versus a similarity variable for different Prandtl numbers (Pr=0.01, 1, 10 and 100). The figures obtained are in agreement with those given in Analysis of Transport Phenomena, William M. Deen, OUP, 1998 (figures 12-5 and 12-6 pages 500-501). This problem was solved originally by Ostrach, 'An analysis of laminar free-convection flow and heat transfer about a flat plate parallel to the direction of the generating body force', NACA Tech. Rept. No 1111, 1953. Derivation of the governing equations can be found in Deen's book.
Please visit this link for a similar treatment using Mathematica: