WillPugs/stat-mech-monte-carlo

2D <m> H=0

WillPugs opened this issue · 0 comments

Having problems fitting the <m> function to the data points.

  • <m> uses conditional statements to return the proper value, the curve_fit function does not work with this sort of function.
  • There is a runtime error when raising certain values to a fractional power. Numpy does not allow raising negative numbers to a fractional power but all the numbers in our array should be positive.
  • There is a runtime/zero-division error when we give sigma values as input to curve_fit.