A=A0=matrix([[6,4,1,1],[4,6,1,1],[1,1,5,2],[1,1,2,5]])
print A
print A0
v11=0 
a21= A0[1,0]
print A
print a21
a31=A0[2,1]
a41=A0[3,1]
print a31,a41
x1= a21*a21 + a31*a31 + a41*a41
print x1
s1= x1^0.5
print s1
x2=(1+a21/s1)/2
print x2
v21= x2^0.5
print v21
v31=(a31*a21)/(2*v21*s1)
print v31
v41=(a41*a21)/(2*v21*s1)
print v41
v1=matrix([[0],[v21],[v31],[v41]])
print v1
I= identity_matrix(4)
print I
p1= I-(2*(v1)*(v1.transpose()))
print p1
A1= p1*A0*p1
print A1

Kreyszig-18.9-1 (last edited 2010-12-17 09:03:45 by punchagan)