#1133 mnewton & Polynomial quotient is not exact

closed
nobody
5
2007-07-18
2007-02-22
Barton Willis
No

(%i1) load(mnewton);
(%i2) mnewton([x^y - %i,x+y-2],[x,y],[1,1.5]);
Polynomial quotient is not exact

I can't eliminate the error by setting gcd : spmod.

Discussion

  • Logged In: YES
    user_id=1179910
    Originator: NO

    Fixed in cvs by replacing the call to linsolve with linsolve_by_lu.

    In current cvs:

    (%i1) load(mnewton)$
    (%i2) mnewton([x^y - %i,x+y-2],[x,y],[1,1.5]);
    (%o2) [[x=0.563519172421*%i+0.451201761049,y=1.5487982389505-0.563519172421*%i]]
    (%i3) abs(rectform(subst(%o2[1], part(%i2, 1))));
    (%o3) [5.58478238525471*10^-16,0.0]

    Andrej

     
    • status: open --> closed