`algsys([expressions],[symbolic variables])`

.
First, all linear equations are solved using the Gauss-method,
then each equation is fed through `solve()`

and the solution used to eliminate one variable in all other
expressions. The equations are treated in the order they are
supplied. This method only works for simple systems. The solution is
provided as vector of solutionvectors, each individual solution in
as linear factor: In the first example below there is one solution
with `xs=-2/3, a2=3/4, a0=2, a1=0`

, the second example has two solutions.
>> syms xs,a0,a1,a2 >> algsys([2-a0,a1-0,a2*xs^2+a1*xs+a0-3-xs, > 2*a2*xs+a1+1],[a2,a1,a0,xs]) ans = [ [ xs+2/3 a2-3/4 a0-2 a1 ] ] >> syms a,xs >> algsys([a*xs+3*a-(3-xs^2),a+2*xs],[a,xs]) ans = [ [ -sqrt(6)+(xs+3) 2*sqrt(6)+(a-6) ] [ sqrt(6)+(xs+3) -2*sqrt(6)+(a-6) ] ] >> float(ans) ans = [[ xs+0.55051 a-1.101 ] [ xs+5.4495 a-10.899 ]]

2009-03-15