A = matrix([[0.49,0.02,0.22],[0.02,0.28,0.20],[0.22,0.20,0.40]])
print A
x0 = matrix([[1],[1],[1]])
print x0
p= x0.transpose()
print p
x1= matrix([[0.890244],[0.609756],[1]])
print x1
x2= matrix([[0.931193],[0.541284],[1]])
print x2
x5 =matrix([[0.990663],[0.504682],[1]])
print x5
x10 =matrix([[0.99707],[0.500146],[1]])
print x10
x15 =matrix([[0.999991],[0.500005],[1]])
print x15
print x15
print p
print x15
print x0
print  A
m1= p*A
print m1
m1 =m1*x0
print m1
m0= p*x0
print m0
z=m1[0,0]
y=m0[0,0]
w= z/y
print w
v1= A*x0
v2= v1.transpose()
m2= v2*v1
print m2
v3= m2[0,0]
v4=v3/y
print v4
v5=w*w
print v5
v6=v4-v5
print v6
v7=v6^(1/2)
print v7

Kreyszig-18.8-1 (last edited 2010-12-17 09:03:57 by punchagan)