Hi Martina
I calculate it with a numerical - algorithm.
regarding LM algorithm, you may have a look at fitOO i wrote several
years ago as a macro tool
http://oooconv.free.fr/fitoo/fitoo_en.html
i did not have time to make it as a real extension but it may contain
some things. feel free to take things/ideas if needed
@mathias, IIRC i calculated it using approximation as
f is my function, p as the unknowns
h=0.01
for i=1 to nbval
derive(i,0)=f(i ,p())
for j=1 to nbparam
memval=p(j)
p(j)=p(j)+h/2
f1=f(i,p())
p(j)=memval-h/2
f2=f(i,p())
derive(i,j)=(f1-f2)/h
p(j)=memval
next j
next i
i had quite some good results with this but maybe weak regarding the
constant h
Laurent
---------------------------------------------------------------------
To unsubscribe, e-mail: [email protected]
For additional commands, e-mail: [email protected]
