In this article, a numerical scheme is introduced for solving the fractional partial differential equation (FPDE) arising from electromagnetic waves in dielectric media (EMWDM) by using an efficient class of finite difference methods. The numerical scheme is based on the Hermite formula. The Caputo’s fractional derivatives in time are discretized by a finite difference scheme of order O((k^(4-alpha)) & O(k^(4-beta)), 1