ich habe einen kleinen Schnipsel für die Taylorreihe des Sinus:
Code: Alles auswählen
def S(x, n):
j = np.reshape(np.arange(n + 1), (-1, 1))
Series = (-1)**j*x**(2*j + 1)/factorial(2*j + 1)
return np.sum(Series, axis=0)
def main():
x = np.array([0, np.pi/4, np.pi/2, np.pi])
n = 10
test = abs(S(x, n) - np.sin(x)) < 1e-3
print(test.all())
import numpy as np
from scipy.special import factorial
if __name__ == '__main__':
main()
Meine Versuche endeten dabei, z. B. für
Code: Alles auswählen
n = np.array([1, 2])
Code: Alles auswählen
j = [[[0], [1]], [[0], [1], [2]]]
Viele Grüße