Higher-order solver for cylindrical shock problem

A fourth-order solver is developed by using the Runge-Kutta method and cubic spline function approximation technique for solving the hyperbolic system of equations with source term. The cylindrical shock wave problem is solved by using the solver. The proposed solver gives better results compared to the previous works as the results are free from more oscillations.