Modelling most of the problems caused by physical, economic, and biological processes leads to differential equations where the exponential function of a matrix, e A is used. For example the matrix exponential appeared in solving the linear first-order ordinary differential equation. Therefore, it is important to present a proper numerical method for approximating the solutions. In 2005, Higham presented an algorithm based on scaling and squaring technique [11].