Can reproduce Sun's results within truncation error.
Now fit intensity instead of F. The weight is given by sqrt(I)=F.