The propagation of electromagnetic waves is modeled by time-dependent Maxwell's equations coupled with constitutive laws that describe the response of the media. In this work, we examine a nonlinear optical model that describes electromagnetic waves in linear Lorentz and nonlinear Kerr and Raman media. To design efficient, accurate, and stable computational methods, we apply high order discontinuous Galerkin discretizations and finite difference schemes in space. The challenge to achieve provable stability for fully-discrete methods lies in the temporal discretizations of the nonlinear terms. To overcome this, novel modification is proposed for the second-order leap-frog and implicit trapezoidal time integrators. The performance of the method is demonstrated via numerical examples.