le péché avec la série de taylor en python

C'est mon code :

import math    
x=float(( input ('x ? ' )))  
n  = 1000   #a big number 
b=0      
for i in range (n):    
   a=(((((-1)**i))*(x**((2*i)+1)))/(math.factorial((2*i)+1)))   
   b+=a     
print (b)

mais il ne fonctionne pas et affiche cette erreur:

"OverflowError: long int too large to convert to float"
  • Quelle est la question? math.factorial(1999) est trop grande pour convertir un float. Il est à environ 10^5733. La valeur max de float est sys.float_info.max, que je vais pari est d'environ 10^308 sur votre système.
  • vous pouvez calculer le as en utilisant la récursivité: a[i] = -a[i-1] x**2 / 2i / (2i + 1)
  • pour examaple quand je veux calculer le péché 30 mon entrée est performer 0.523 . mais il les "OverflowError: long int trop grand pour convertir flotter"
  • avez-vous vu mon dernier commentaire ici ?!
  • L'entrée est hors de propos, votre code calcule factorial(1999) (ainsi qu'un tas d'autres numéros qui sont trop gros pour convertir float), indépendamment de la valeur d'entrée. Lorsque vous divisez un float par un nombre entier, Python essaie de convertir l'entier à un flotteur. Il échoue.
InformationsquelleAutor Hippo | 2013-11-14
  1. 2

    C'est une implémentation possible:

    def mysin(x, order):
    a = x
    s = a
    for i in range(1, order):
        a *= -1 * x**2 / ((2 * i) * (2 * i + 1))
        s += a
    return s

    C'est juste pour le traçage:

    import numpy as np
vmysin = np.vectorize(mysin, excluded=['order'])

x = np.linspace(-80, 80, 500)
y2 = vmysin(x, 2)
y10 = vmysin(x, 10)
y100 = vmysin(x, 100)
y1000 = vmysin(x, 1000)
y = np.sin(x)

import matplotlib.pyplot as plt
plt.plot(x, y, label='sin(x)')
plt.plot(x, y2, label='order 2')
plt.plot(x, y10, label='order 10')
plt.plot(x, y100, label='order 100')
plt.plot(x, y1000, label='order 1000')
plt.ylim([-3, 3])
plt.legend()
plt.show()

    le péché avec la série de taylor en python

    Il souffre d'instabilité numérique et de dépassement de capacité, car après un certain temps (~100 boucles, dependig sur x) a devient 0.

    InformationsquelleAutor Ruggero Turra
  2. 1

    réponse correcte des sinus de la série de taylor

    # calculate sin taylor series by using for loop in python
from math import*

print "sine taylor series is="

x=float(raw_input("enter value of x="))

for k in range(0,10,1):
    y=((-1)**k)*(x**(1+2*k))/factorial(1+2*k)

    print y
    • La série de taylor est une somme de sorte que vous devez initialiser votre y à zéro et ajouter y += ... dans la boucle. Sinon, en dehors de quelques fautes de frappe, il semble bon.
    InformationsquelleAutor Mehdi raza