The magnetic field returned by get_sphere_field_z does not, at present, take into account the radial component of the field in cylindrical coordinates.

After 210 steps some of the can radii are negative and others are on the order of 1e18. This causes find_mutual_inductance to conk out.

Define and implement a function that returns the total resistance of driving coil.

The labels displayed by the method have no bearing on what is actually being simulated.

Allow the user to specify an Hc field for the superconductor under test.

Some of the currents calculated by calculate_current blow up exponentially. Here's a sample output

array([ 3.74652503e-05, 1.56855283e+02, -1.56188786e+02,
-7.84928314e+00, -7.84928314e+00, -1.56188786e+02,
1.56855283e+02])
array([ 2.58594657e-03, 1.20628550e+04, -1.46167635e+04,
2.04974350e+03, 2.04974350e+03, -1.46167635e+04,
1.20628550e+04])
array([ 4.63675506e-01, 2.58176042e+06, -3.63077052e+06,
9.47461236e+05, 9.47461236e+05, -3.63077052e+06,
2.58176042e+06])
array([ 2.95760429e+06, 2.40968714e+13, -3.82953652e+13,
1.32516296e+13, 1.32516296e+13, -3.82953652e+13,
2.40968714e+13])
array([ 3.02247784e+27, 2.99921103e+34, -5.20647322e+34,
2.09618219e+34, 2.09618219e+34, -5.20647322e+34,
2.99921103e+34])
array([ 6.04807121e+90, 7.29115901e+97, -1.33921318e+98,
5.84036823e+97, 5.84036823e+97, -1.33921318e+98,
7.29115901e+97])
main:101: RuntimeWarning: overflow encountered in double_scalars
main:126: RuntimeWarning: overflow encountered in square
array([ 8.51065105e+280, 1.17303217e+288, -2.23322018e+288,
1.01943075e+288, 1.01943075e+288, -2.23322018e+288,
1.17303217e+288])
array([ nan, nan, nan, nan, nan, nan, nan])
Traceback (click to the left of this block for traceback)
...
ZeroDivisionError: Symbolic division by zero

Placing a comment after a simple name-only variable output in a cell will cause the variable output to not happen.

The simulator hangs up in the sagecloud. There were two changes, unfortunately. I agreed to a sagecloud update yesterday and didn't immediately run the sim. I also modified the array that holds the current for graphing to initialize to zeroes.

This needs to be done to make the code more readable before the difference between 3 and 6 is obscured. It will also make the code more usable.

All arrays in IDL are initialized to zero. The original code may depend on this. Learn how to initialize arrays in numpy to zero easily.