from numpy import mat
A = mat([[2.04,-1,0,0],[-1,2.04,-1,0],[0,-1,2.04,-1],[0,0,-1,2.04]])
B = mat([[40.8],[0.8],[0.8],[200.8]])
g = A[0,1]
f1 =A[0,0]
e2 = A[1,0]
f2 = A[0,0] - e2 * A[1,0]
e3 = A[1,0]/f2
f3 = A[0,0] - e3 * A[1,0]
e4 = A[1,0]/f3#
f4 = A[0,0] - e4 * A[1,0]
M = mat([[f1,g,0,0],[e2,f2,g,0],[0,e3,f3,g],[0,0,e4,f4]])
L = mat([[1,0,0,0],[M[1,0],1,0,0],[0,M[2,1],1,0],[0,0,M[3,2],1]])
U = mat([[M[0,0],g,0,0],[0,M[1,1],g,0],[0,0,M[2,2],g],[0,0,0,M[3,3]]])
r1 = B[0,0]
r2 = B[1,0] - e2*B[0,0]
r3 = B[2,0] - e3*r2#
r4= B[3,0] - e4*r3#
N = mat([[r1],[r2],[r3],[r4]])
T4 = r4/U[3,3]
T3 = (r3 - g*T4)/U[2,2]
T2 = (r2 - g*T3)/U[1,1]
T1 = (r1 - g*T2)/U[0,0]
print "T1 = ",T1
print "T2 = ",T2
print "T3 = ",T3
print "T4 = ",T4
from numpy import mat,sqrt
A = mat([[6,15,55],[15,55,225],[55,225,979]])
sl = 0
l11 = sqrt(A[0,0])
#for second row
l21 = (A[1,0])/l11
l22 = (A[1,1] - l21**2)**(0.5)
#for third row
l31 = (A[2,0])/l11#
l32 = (A[2,1] - l21*l31)/l22#
l33 = (A[2,2] - l31**2 - l32**2)**(0.5)#
L = mat([[l11,0,0],[l21,l22,0],[l31,l32,l33]])
print "\nL = \n",L
#3x - 0.1y - 0.2z = 7.85
#0.1x + 7y - 0.3z = -19.3
#0.3x - 0.2y + 10z = 71.4
Y = 0#
Z = 0#
x=range(1,3)
y=range(1,3)
z=range(1,3)
for i in range(1,3):
x[i-1]= (7.85 +0.1*Y+0.2*Z)/3#
X = x[i-1]
y[i-1]= (-19.3 - 0.1*X +0.3*Z)/7#
Y = y[i-1]#
z[i-1]= (71.4 - 0.3*X+0.2*Y)/10#
Z = z[i-1]
if i==2:
ex = (x[i-1] - x[(i-2)])*100/x[i-1]
ey = (y[i-1] - y[i-2])*100/y[i-1]
ez = (z[i-1] - z[i-2])*100/z[i-1]
print "x through two iterations =",x[0:1]
print "y through two iterations =",y[0:1]
print "z through two iterations =",z[0:1]
print "error of x = ",ex,"%"
print "error of y = ",ey,"%"
print "error of z = ",ez,"%"
from __future__ import division
from numpy import mat
A = mat([[1,0.5,1/3],[1,2/3,1/2],[1,3/4,3/5]])
B = mat([[1.833333],[2.166667],[2.35]])
U = A**-1
X = U*B#
x = X[0,0]
y = X[1,0]
z = X[2,0]
print "x = ",x
print "y = ",y
print "z = ",z
from __future__ import division
from numpy import mat
A = mat([[1,0.5,1/3],[1,2/3,1/2],[1,3/4,3/5]])
B = mat([[1.833333],[2.166667],[2.35]])
U = A**-1
X = U*B#
x = X[0,0]
y = X[1,0]
z = X[2,0]
print "x = ",x
print "y = ",y
print "z = ",z
from numpy import mat
A = mat([[1,0.5,1/3],[1/2,1/3,1/4],[1/3,1/4,1/5]])
B = mat([[1.833333],[1.083333],[0.783333]])
U = A**-1
X = U*B#
x = X[0,0]
y = X[1,0]
z = X[2,0]
print "x = ",x
print "y = ",y
print "z = ",z