First , I was obtained density of magnetization and density of energy curves of isotropic Heisenberg model on square lattice (2x2) using worm and dirloop_sse.

After I was obtained these curves with wolfram mathematica.

Curves for magnetization density(density of energy) obtained with mathematica and worm(sse) is strongly differ one from another.They must match.

I was attached the files of mathematica and Heis2x2.py.