I am using linsolve() for the solution of linear equations. In the example
below I want to express the solution as:
xsol = A . r
where r = matrix([%r1],[%r2])
and A = matrix([-a13/a11,-a12/a11],[0,1],[1,0])
Its not clear how to construct or extract the matrices A and r from the
list ysol.
Any assistance would be appreciated.
TIA
C. Frangos.
(%i513) ysol : linsolve([a11*x1+a12*x2+a13*x3],[x1,x2,x3]);
(%o513) [x1 = -(%r1*a13+%r2*a12)/a11,x2 = %r2,x3 = %r1]
(%i514)