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

Compare commits

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

2 commits
57bb3a3368 ... 0d8f8ef6cc

Author SHA1 Message Date
Manuel Thalmann 0d8f8ef6cc Fix algorithm 2023-03-29 13:03:17 +02:00
Manuel Thalmann 2537f32208 Remove duplicated code 2023-03-29 12:35:11 +02:00
1 changed files with 6 additions and 8 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,15 +26,13 @@ def Miguel_S5_Aufg2(x, y, xx):
A[i][0] = 2 * (h[0] + h[1])
A[i][1] = h[1]
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:
A[i][i - 1] = h[i - 1]
A[i][i] = 2 * (h[i - 1] + 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]))
A[i][i - 1] = h[i]
A[i][i] = 2 * (h[i] + h[i + 1])
z[i] = 3 * (((y[i + 2] - y[i + 1]) / h[i + 1]) - ((y[i + 1] - y[i]) / h[i]))
if i < n - 3:
A[i][i + 1] = h[i + 1]
c[1:n - 1,:] = linalg.solve(A, z)