Laplace transform has been used for solving differential equations of fractional order either PDEs or ODEs. However, using the Laplace transform sometimes leads to solutions in Laplace space that are not readily invertible to the real domain by analytical techniques. Therefore, numerical inversion techniques are then used to convert the obtained solution from Laplace domain into time domain. Various famous methods for numerical inversion of Laplace transform are based on quadrature approximation of Bromwich integral. The key features are the contour deformation and the choice of the quadrature rule. In this work, the Gauss–Hermite quadrature method and the contour integration method based on the trapezoidal and midpoint rule are tested and evaluated according to the criteria of applicability to actual inversion problems, applicability to different types of fractional differential equations, numerical accuracy, computational efficiency, and ease of programming and implementation. The performance and efficiency of the methods are demonstrated with the help of figures and tables. It is observed that the proposed methods converge rapidly with optimal accuracy without any time instability.