Doing a run with the default parameters yields the following figures and the following gain of
1.51%, which is small. Also, what is the MHD term in the ODE function? Also, how to properly
calculate gain? I tried doing E_plasm Initial-E_plasm final but got a low number, and did it again
with the total energy and got a similar number. I also did it with the yield and got like 0.49%.