manuth/ZHAWNotes
1
0
Fork 0
Code Issues Pull requests Projects Releases Wiki Activity

Compare commits

base: manuth:0d8f8ef6cc7b54eeca789cdc5283918a4a081215
Branches Tags
manuth:main
..
compare: manuth:57bb3a3368472c4facdbaffc940e6567d9ff5bd8
Branches Tags
manuth:main

No commits in common. "0d8f8ef6cc7b54eeca789cdc5283918a4a081215" and "57bb3a3368472c4facdbaffc940e6567d9ff5bd8" have entirely different histories.
0d8f8ef6cc ... 57bb3a3368

1 changed files with 8 additions and 6 deletions
Show stats Expand all files Collapse all files

14
Notes/Semester 4/HM2 - Höhere Mathematik 2/Week 5/Miguel_S5_Aufg2.py
View file

@ -26,13 +26,15 @@ def Miguel_S5_Aufg2(x, y, xx):
A[i][0] = 2 * (h[0] + h[1]) A[i][0] = 2 * (h[0] + h[1])
A[i][1] = h[1] A[i][1] = h[1]
z[i] = 3 * (((y[2] - y[1]) / h[1]) - ((y[1] - y[0]) / h[0])) z[i] = 3 * (((y[2] - y[1]) / h[1]) - ((y[1] - y[0]) / h[0]))
elif i == n - 3:
A[i][i - 1] = h[i - 1]
A[i][i] = 2 * (h[i - 1] + h[i])
z[i] = 3 * (((y[i + 1] - y[i]) / h[i]) - ((y[i] - y[i - 1]) / h[i - 1]))
else: else:
A[i][i - 1] = h[i] A[i][i - 1] = h[i - 1]
A[i][i] = 2 * (h[i] + h[i + 1]) A[i][i] = 2 * (h[i - 1] + h[i])
z[i] = 3 * (((y[i + 2] - y[i + 1]) / h[i + 1]) - ((y[i + 1] - y[i]) / h[i])) A[i][i + 1] = h[i]
z[i] = 3 * (((y[i + 1] - y[i]) / h[i]) - ((y[i] - y[i - 1]) / h[i - 1]))
if i < n - 3:
A[i][i + 1] = h[i + 1]
c[1:n - 1,:] = linalg.solve(A, z) c[1:n - 1,:] = linalg.solve(A, z)