a = array([[2*R0*(Q0**2), 2*R1*(Q1**2), 2*R2*(Q2**2), 2*R3*(Q3**2), 2*R4*(Q4**2), 0, 0, 0],
[0, 0, 0, 2*R3*(Q3**2), 2*R4*(Q4**2), 0, 0, 0], [0, 2*R1*(Q1**2),0, 0, 0, 2*R5*(Q5**2), 2*R6*(Q6**2), 2*R7*(Q7**2)], [1, 1, 0, 0, 0, 0, 1], [0, 1, 1, 0, 0, 1, 0, 0, 0], [0, 0, 1, 1, 1, 0, 0, 0],
[0, 0, 0, 0, 0, 1, 0, 1], [0, 0, 0, 0, 0, 0, 0, 1, 1]])
b = array([[95-(R0*(Q0**3))-(R1*(Q1**3))-(R2*(Q2**3))-(R3*(Q3**3))-(R4*(Q4**3))], [-(R4*(Q4**3))-(R3*(Q3**3))],
[-(R6*(Q6**3))-(R7*(Q7**3))-(R1*(Q1**3))-(R5*(Q5**3))], [-Q0+Q1+Q6], [0.6-Q1+Q2+Q5], [-Q2+Q3+Q4], [1.4-Q7-Q5], [Q7-Q6]])
p = solve(a, b)
print ("deltaQ:", p)
Solve (a,b) is a gauss jordan elimination function that solves linear equations where R0-R7 AND Q0-Q7 ARE CONSTANTS! It works with other matrices, however for this one i get an error.
a = np.matrix(a)andb=np.matrix(b)