On 13/08/2010 20:18, Dieter Kaiser wrote:
> Am Freitag, den 13.08.2010, 18:54 +0100 schrieb Leo Butler:
>
>> On Wed, 11 Aug 2010, Sergio Callegari wrote:
>>
>> < Hi,
>> <
>> < just a short bug report about mnewton
>> <
>> < %i1) load(mnewton);
>> < (%o1) /usr/share/maxima/5.22.0/share/contrib/mnewton.mac
>> < (%i2) mnewton([2*a^a-5],[a],[1]);
>> <
>> < Improper argument to coefmatrix:
>> < (2.0 transpose(h ) + [ - 3.0 ])
>> < 1
>> < 1
>> < #0:
>> < solve_by_lu(eqs=(2.0*'transpose(h[1])+matrix([-3.0]))[1],vars=
>> < [h[1]],field=complexfield)(mnewton.mac
>> < line 51)
>> < #1: mnewton(funclist=[2*a^a-5],varlist=[a],guesslist=[1])
>> < (mnewton.mac line 113)
>> < -- an error. To debug this try: debugmode(true);
>> <
>> < Since mnewton([2*a^a-5],[a],[1]) is part of the examples for mnetwon
>> < in the documentation, I guess there is probably some regression.
>> <
>> < Note that
>> <
>> < mnewton([2*a^a-5,x],[a,x],[1,0]);
>> <
>> < i.e. adding a dummy variable works.
>>
>> Try:
>>
>> mnewton([2*a^a-5],[a],[1]), scalarmatrixp=false;
>>
> This problem is not present in Maxima 5.22.1. A patch to mnewton.mac has
> been committed some days ago.
>
> Dieter Kaiser
>
>
Thanks,
I did not notice that 5.22.1 was coming out with the fix.
Sergio