from scipy import *
from scipy import linalg
A = mat("[1 3 5; 2+1j 5 1; 2 3 8]")
b = mat("[10;8;3]")
print A.I*b # .I means to take the inverse
print linalg.solve(A, b)
[[-8.21529745+2.95750708j]
[ 5.39660057+0.6572238j ]
[ 0.40509915-0.98583569j]]
[[-8.21529745+2.95750708j]
[ 5.39660057+0.6572238j ]
[ 0.40509915-0.98583569j]]
and we have the solutionsolve -3*(2 A1) - 2 A1 - (A1 - A2) = 0 , -1 + 3*(2 A1) - (A2 - A1) - 4 A2 = 0
With or without Mathematica 8, wolframalpha is a great engine! With Mathematica 8 it's awesome!!!This site can solve simultaneous equations and much much more
http://www.wolframalpha.com/
For example
-3*(2 I1) - 2 I1 - (I1 - I2) = 0 (1)
-1 + 3*(2 I1) - (I2 - I1) - 4 I2 = 0 (2)
So we type
and we have the solution
http://www.wolframalpha.com/input/?...2)+=+0+,+-1+++3*(2+A1)+-+(A2+-+A1)+-+4+A2+=+0