Rogers, Raymond escribi?:
>> I want to get something like this:
>>
>> (%i58) find_root(f(x),x,0.3,0.7);
>> (%o58) 0.52360846201579
>>
>> --
>>
> You might try implementing Newton-Raphson.
> http://numericalmethods.eng.usf.edu/ebooks/newtonraphson_03nle_ebook.htm
> I'm really surprised that Maxima doesn't have a package to do this.
load(mnewton)$
numer:true$
mnewton([f(x)=0],[x],[1]);
also works.
If numer equals false, the above mnewton call gives an error. I think we
could set 'numer:true' locally in solve_by_lu in mnewton.mac.
Or is there a better way for fixing this?
--
Mario Rodriguez Riotorto
www.biomates.net