Correct an serious error in the previous code. Previous code doesn’t
change the volume at different steps. This is wrong for NPT ensemble
where the volume is constantly changing. This will cause a non-zero
g(r) at small distance.