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

Fix algorithm

Browse source
This commit is contained in:
Manuel Thalmann 2023-03-29 13:03:17 +02:00
parent 2537f32208
commit 0d8f8ef6cc
1 changed files with 6 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -27,10 +27,12 @@ def Miguel_S5_Aufg2(x, y, xx):
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]))
else: else:
A[i][i - 1] = h[i - 1] A[i][i - 1] = h[i]
A[i][i] = 2 * (h[i - 1] + h[i]) A[i][i] = 2 * (h[i] + h[i + 1])
A[i][i + 1] = h[i] z[i] = 3 * (((y[i + 2] - y[i + 1]) / h[i + 1]) - ((y[i + 1] - y[i]) / 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)