PY
py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
#Stack Over Flow
#thanks to izidor
# link https://goo.gl/zRmH9m,
def linearsolver(A,b):
n = len(A)
M = A
i = 0
for x in M:
x.append(b[i])
i += 1
for k in range(n):
for i in range(k,n):
if abs(M[i][k]) > abs(M[k][k]):
M[k], M[i] = M[i],M[k]
else:
pass
for j in range(k+1,n):
q = M[j][k] / M[k][k]·
for m in range(k, n+1):
M[j][m] += q * M[k][m]
x = [0 for i in range(n)]
print "n = ", n
Enter to Rename, Shift+Enter to Preview
OUTPUT
Run