Julia LsqFit is giving wrong results whereas I am able to get correct results using scipy.optimize (almost equal to exp results)

Question

I tried to fit Magnetization vs Magnetic field using law of approach to saturation magnetization using Julia, but not able to get correct results. The parameters, p[2] is showin