hmm. I am glad this works, thanks for the help. But according to my prof, the values being outputted seem to be off. For instance, it should equals closer to -2 at T=1, and closer to zero at T=4. The function I am trying to emulate is the energy per spin in here
I believe, according to the document, that the function should be:
I think that the derivative may be wrong...Any ideas?
It doesnt leave an answer for me using this code (@ e(1)), just a giant mess...