die Sünde mit der taylor-Reihe in python

Dies ist mein 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)

aber es funktioniert nicht und zeigt diesen Fehler:

"OverflowError: long int too large to convert to float"
  • Was ist die Frage? math.factorial(1999) ist zu groß, um eine Konvertierung in eine float. Es ist etwa 10^5733. Der max-Wert von float ist sys.float_info.max, die ich Wette, der ist etwa 10^308 auf Ihrem system.
  • Sie können die Berechnung der as mit Rekursion: a[i] = -a[i-1] x**2 / 2i / (2i + 1)
  • für examaple, wenn ich will berechnen Sie sin 30 mein input ist 0.523 . es zeigt aber, das "OverflowError: long int zu groß, zu konvertieren, zu schweben"
  • habt Ihr meinen letzten Kommentar hier ?!
  • Die Eingabe ist irrelevant, dein code berechnet factorial(1999) (sowie ein paar andere zahlen, die sind zu groß, um zu konvertieren, um float) unabhängig von der input-Wert. Wenn Sie teilen ein float durch eine ganze Zahl, Python versucht zu konvertieren integer zu float. Es fehlschlägt.
InformationsquelleAutor Hippo | 2013-11-14
  1. 2

    Dies ist eine mögliche Umsetzung:

    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

    Dies ist nur für das zeichnen:

    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()

    die Sünde mit der taylor-Reihe in python

    Leidet Sie unter numerische Instabilität und Unterlauf, da nach einer Weile (~100 loops, dependig auf x) a wird, wird zu 0.

    InformationsquelleAutor Ruggero Turra
  2. 1

    richtige Antwort von Sinus taylor-Reihe

    # 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
    • Die taylor-Reihe ist eine Summe, so müssen Sie zu initialisieren, Ihr y zu null, und fügen Sie y += ... in der Schleife. Ansonsten, abgesehen von ein paar Tippfehler, es sieht gut aus.
    InformationsquelleAutor Mehdi raza
Schreibe einen Kommentar