1 Jan 09:01 2008

### Re: Solving linear algebraic eqns. in vector (matrix) notation ?

Radovan Omorjan <omorr <at> uns.ns.ac.yu>

2008-01-01 08:01:59 GMT

2008-01-01 08:01:59 GMT

I've figured out myself how to solve this. It took me some time to realise that equations and unknowns have to be in lists in order to solve them with "solve" or "linsolve". I also realised that I do not need "eigen". Something like transforming matrix m into equation list like: (%i9) eqns : makelist(m[i][1]=0,i,1,6); (%o9) [- n2 + n1 - eps = 0, - n2 + n1 - eps + eps = 0, 1 1 1 2 2 2 1 eps eps 2 1 - n2 + n1 + eps = 0, - n2 + n1 + eps = 0, - n2 + n1 - ---- - ---- = 0, 3 3 2 4 4 1 5 5 2 2 n1 - n2 = 0] 6 6 and solving by: (%i10) linsolve(eqns,n2); (%o10) [n2 = n1 - eps , n2 = n1 - eps + eps , n2 = n1 + eps , 1 1 1 2 2 2 1 3 3 2 - 2 n1 + eps + eps 5 2 1 n2 = n1 + eps , n2 = - ---------------------, n2 = n1 ] 4 4 1 5 2 6 6 will do the job. I should have read the manual more carefuly, sorry. Anyway, I've learned some new things.(Continue reading)