Solutions such as symmetric, periodic, and solitary wave solutions play a significant role in the field of partial differential equations (PDEs), and they can be utilized to explain several phenomena in physics and engineering. Therefore, constructing such solutions is significantly essential. This article concentrates on employing the improved exp(−ϕ(η))-expansion approach and the method of lines on the variant Boussinesq system to establish its exact and numerical solutions. Novel solutions based on the solitary wave structures are obtained. We present a comprehensible comparison between the accomplished exact and numerical results to testify the accuracy of the used numerical technique. Some 3D and 2D diagrams are sketched for some solutions. We also investigate the L2 error and the CPU time of the used numerical method. The used mathematical tools can be comfortably invoked to handle more nonlinear evolution equations.