The paper investigates the invariance and integrability properties of Hirota-Satsuma equations. Painleve analysis for the general  similarity reduced ordinary differential equation is performed.  Using Rung-Kutta-Merson method in shooting and matching  technique, the nonlinear ordinary differential equations are solved  that were numerically converted from similarity reduction.