A numerical algorithm is presented for studying Marangoni convection flow over a flat plate with an imposed temperature distribution. Plate temperature varies with x in the following prescribed manner: where A and k are constants. By means of similarity transformation, the original nonlinear partial differential equations of flow are transformed to a pair of nonlinear ordinary differential equations. Subsequently they are reduced to a first order system and integrated using Newton Raphson and adaptive Runge-Kutta methods. The computer codes are developed for this numerical analysis in Matlab environment. Velocity profiles for various values of k, and temperature profiles for various Prandtl number and k are illustrated graphically. The results of the present simulation are then compared with the previous works available in literature with good agreement.