Numerical Investigations of the Coupled Nonlinear Non-homogeneous Partial Differential Equations