daknuett/scientific-programming-exercises
2
1
Fork 0
Code Issues Pull Requests Releases Wiki Activity

fixed ex42

Browse Source
This commit is contained in:
Daniel Knüttel 2019-01-23 16:21:46 +01:00
parent 8c2e6145a8
commit c5a686cbed
1 changed files with 7 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
ex_42.py
View File

@ -10,10 +10,13 @@ A = np.array(
, [0, 0, 0, 0, 1, 0] , [0, 0, 0, 0, 1, 0]
]) ])
eigvalues, eigvectors = np.linalg.eig(A) eigvalues, eigvectors = np.linalg.eigh(A)
print(eigvalues)
print(eigvectors)
l_max_i = eigvalues.argmax() l_max_i = eigvalues.argmax()
l_max = eigvalues[l_max_i] l_max = eigvalues[l_max_i]
print(l_max_i, l_max)
v_max = eigvectors[l_max_i] v_max = eigvectors[l_max_i]
def some_norm(M): def some_norm(M):
@ -22,8 +25,8 @@ def some_norm(M):
B = A B = A
for k in range(1, 21): for k in range(1, 21):
B = B.dot(A) B = B.dot(A)
print(B) print("some kind of error for k =", k, ":", some_norm(B - l_max**k * np.outer(v_max, v_max)) / l_max**k)
print(l_max**k * np.outer(v_max, v_max))
print("some kind of error for k =", k, ":", some_norm(B - l_max**k * np.outer(v_max, v_max)))
print(v_max.argmax())