High order FDTD methods for electromagnetic systems in dispersive inhomogeneous media