Hi,
Folks in physics are probably familiar with the scaling ansatz for finite-size scaling phase transition $y=L^{\xi/\nu}f[(x-x_c)L^{1/\nu}]$. The function $f$ is generally unknown but the parameters can be obtained by tunning so that the data $y(x)$ at different $L$ collapse onto each other. Does anyone know how to obtain these parameters numerically? A Mathematica package will be much appreciated. Thanks.