and yet another method:
(%i2) A : matrix([1,4],[7,4])$
(%i3) b : [1,5]$
(%i4) x: [x1,x2]$
(%i5) solve (makelist (A[i].x - b[i], i, 1, 2));
(%o5) [[x2 = 1/12,x1 = 2/3]]
Jaime
On 05/04/2013 08:37 AM, Mario Rodriguez wrote:
> El s?b, 04-05-2013 a las 11:12 +1000, Alasdair McAndrew escribi?:
>> Given an n x n matrix A, and an n x 1 matrix b, is there a way of
>> solving
>>
>>
>> Ax = b
>>
>>
>> other than invert(A).b?
> A : matrix([1,4],[7,4]);
> B : [1,5];
> linsolve_by_lu(A,B);
>
>
>
>
>>
>> If I could append b to A, then I could apply "echelon" to the
>> augmented matrix. However, I can't find a command which allows me to
>> append, or stack matrices.
>>
> addcol(A,B);
>
>