In this article, generalized Kudryashov technique has been implemented to construct new soliton solutions of cubic KdV Boussinesq equation. Also, the OHAM has been employed to compute the numerical solution. The solutions thus acquired by the above methods are illustrated graphically. The techniques considered here are efficacious, plausible and can be employed in mathematical physics to compute new exact and numerical solutions of NPDEs.