The time-dependent Ginzburg-Landau model of superconductors consists of coupled nonlinear partial differential equations, which presents difficulties in the numerical solution. We present an ...