Numerical computation of the interaction between electromagnetic waves and nonlinear superconducting materials