ich habe wieder ein Problem beim Programmieren. Ich muss mit der Formel f"(x0) ≈ ( f(x0) − 2f(x0 + h) + f(x0 + 2h)) /h^2 die zweite Ableitung von tan(pi/4) bestimmen. Für die Werte h=1e − 1, 1e − 2, ..., 1e − 9, 1e − 10 einsetzen. Anschließend soll ein Graph gezeichnet werden an der der Unterschied zwischen expected error (O(h))und absolute error für die Werte für h zu sehen sind.
Momentan habe ich diesen Code zusammengestellt:
Code: Alles auswählen
import math
import numpy
x= math.pi/4
fx= math.tan(math.pi/4)
fx2t=0 #2nd derivative theoretically
hs= [1e-1, 1e-2, 1e-3, 1e-4, 1e-5, 1e-6, 1e-7, 1e-8, 1e-9, 1e-10]
for h in hs:
fx2= (fx- 2*math.tan((math.pi+h)/4)+ math.tan((math.pi+2*h)/4)) / h**2 #2nd derivative numerically
print(fx2)
Ich würde mich über eine Antwort freuen.
Viele Grüße