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

Fix incorrect algorithm

Browse source
This commit is contained in:
Manuel Thalmann 2023-03-29 08:40:25 +02:00
parent 483f4d1c25
commit 57bb3a3368
1 changed files with 1 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -28,7 +28,7 @@ def Miguel_S5_Aufg2(x, y, xx):
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: elif i == n - 3:
A[i][i - 1] = h[i - 1] A[i][i - 1] = h[i - 1]
A[i][i] = h[i - 1] + h[i] 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])) 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 - 1] A[i][i - 1] = h[i - 1]
@ -63,6 +63,5 @@ y = array([6, 3, 9, 0])
xx = array([4.5, 6.5, 9.0]) xx = array([4.5, 6.5, 9.0])
Miguel_S5_Aufg2(x, y, xx) Miguel_S5_Aufg2(x, y, xx)
pyplot.scatter(xx, interpolate.CubicSpline(x, y)(xx))
# %% # %%