This project provides how the method of characteristics is applied to derive the fundamental solution of linear homogeneous wave equation and how shock waves and shock solutions are formed by a discontinuous initial condition. Moreover, it provides how the solution of the linear non-homogeneous wave equation is determined by combining the split equations, the homogeneous wave equation with non-zero initial data and the non-homogeneous wave equation with zero initial data using method of characteristics.