previous: integration
next: read data file

11: Solve matrix equation Ax=b

Solve Ax=b (and check the solution) given that

A = 
1   2   0
2   5  -1
4  10  -1
      b = 
7
8
9

Mathematica:

A = { {1,2,0}, {2,5,-1}, {4,10,-1} }
b = {7,8,9}
x = LinearSolve[ A, b ]
A . x     (does equal b)

Matlab:

A = [ 1 2 0 ; 2 5 -1 ; 4 10 -1 ]
b = [ 7 8 9 ]
x = A\transpose(b)
A * x     (does equal b)

Maple:

xxx

IDL:

A = [ [1,2,0], [2,5,-1], [4,10,-1] ]
b = [7,8,9]
svd, transpose(A), w,u,v
svbksb, u,w,v,b,x
print, x
print, x # A   (does equal b)