Maxima not solving a linear system



>>>>> "Sergio" == Sergio Rossi <deltaquattro at hotmail.it> writes:

    Sergio> 1.  (*) text/plain          ( ) text/html           
    Sergio> Hello,
 
    Sergio> apparently this simple linear system cannot be solved with Maxima 5.9.3 under WinXP :
 
    Sergio> (%i1) eq1: a0 + x1*a1 + y1*a2 + x1^2*a3+y1^2*a4+x1*y1*a5 = f1;
    Sergio> (%i2) eq2: a0 + x2*a1 + y2*a2 + x2^2*a3+y2^2*a4+x2*y2*a5 = f2;
    Sergio> (%i3) eq3: a0 + x3*a1 + y3*a2 + x3^2*a3+y3^2*a4+x3*y3*a5 = f3;
    Sergio> (%i4) eq4: a0 + x4*a1 + y4*a2 + x4^2*a3+y4^2*a4+x4*y4*a5 = f4;
    Sergio> (%i5) eq5: a0 + x5*a1 + y5*a2 + x5^2*a3+y5^2*a4+x5*y5*a5 = f5;
    Sergio> (%i6) eq6: a0 + x6*a1 + y6*a2 + x6^2*a3+y6^2*a4+x6*y6*a5 = f6;

    Sergio> (%i7) linsolve([e1,eq2,eq3,eq4,eq5,eq6],[a0,a1,a2,a3,a4,a5]);

Perhaps you meant to say eq1 instead of e1 here?  CMUCL produces a
solution when you replace e1 with eq1.  (With e1 instead, CMUCL gets
an error trying to print out message saying the equations are
inconsistent or something.  Need to fix that.)


The Windows version comes with gcl, so if it still can't solve this
equation, it's because gcl has too small a heap to run this.
 
    Sergio> PS is it possible to insert a linear system to be solved
    Sergio> in the form CoefficientMatrix*VectorOfUnknowns, rather
    Sergio> than in the form [Eq1, Eq2, ...], [Unk1, Unk2,...]?
 
You can probably do it with Maxima's matrix support, but I'm not
familiar with that.

Ray