frames = 1000 if anim: Ns = (300, 300) side_length = (1 / 3.) square_coaxial_cable = Shape(Ns, 1.5, (0.5, 0.5), 3e-1, side_length, shape='square', filled=True) cable = System(Ns) cable.add(square_coaxial_cable) all_potentials = cable.SOR_anim(w=1.9999999999999999999999999995, tol=0, max_iter=frames) plt.ioff() fig, ax = plt.subplots(figsize=(15, 15)) image = plt.imshow(all_potentials[0].T, vmin=np.min(all_potentials), vmax=np.max(all_potentials), interpolation='none', aspect='equal', extent=None, origin='lower') iter_text = ax.text(0.05, 0.95, '', transform=ax.transAxes, size=20,