Nonlinear shock waves in plasma was modeled and studied using Ramani equation of sixth order and its coupled form representing the interaction between two waves. A new combined methodology of both Lie infinitesimal transformation and singular manifold methods was exploited to create analytical solutions. The method was extended to investigate a coupled system of evolution equations as same as a single evolution equation. Soliton solutions were created for both models.