K-te Ableitung in Python [geschlossen] ⇐ Python

[Guest]

Mein Code für die k-te Ableitung unter Verwendung der Vorwärtsdifferenz zur numerischen Differenzierung funktioniert nicht. Es gibt mir -38, nicht -2,8, und ich kann das Problem nicht finden

Code: Select all

#kth order derivative

#f'(x) = 1/h * (yn+1 - yn)
#looks nice,untested

import math as mt
import numpy as np

k = 4
A = 0
B = 2.5
h = 0.5

nodes =int((B - A) / h) + 1
xn = np.linspace(A, B, nodes)
yn = (np.cos(xn))**2

print(xn)
print(yn)

def KthDeriv (A, B, k, h):
nodes =int((B - A) / h) + 1
xn = np.linspace(A, B, nodes)
yn = (np.cos(xn))**2
derivative = 0
for i in range(k + 1):  #if we want first deriv we need two points
if i + k < len(yn):
term = (-1)**i * mt.comb(k, i) * yn[i + k]
derivative += term
derivative /= h**k
return derivative
test = KthDeriv (A, B, k, h)
print(test)