Introduce error calculation
Put error calculation in effective mass
Divide by sqrt(n)
Remove debug output
Add χ² test to fit function