problem with fitting data using lsquares



On 6/11/10, Hugo Coolens <coolens at kahosl.be> wrote:

> m:matrix([100,16000],[200,8000],[500,3636],[1000,1818],[2000,870],[5000,323],[10000,160],[20000,83],[50000,33],[10000,20]);
> load (lsquares)$
> lsquares_estimates
> (m,[x,y],y=sqrt((a^2)+((10e7/(2*pi*x*b))^2)),[a,b],initial=[10,10]);

I think you want %pi instead of pi there. Otherwise the expression
won't evaluate to a number.

best,

Robert Dodier