An old thread, but I catch it anyways.
While the quick demonstration runs without error in GNU Octave (no matter what version I use), it doesn't result in a proper fit and the fitting constants are increadibly large (6.9415e+21 2.4425e+11 -7.7388e+21). Any idea why that is happening?