Python中使用Maclaurin系列逼近e ^ x

时间:2014-07-23 15:49:11

标签: python series expansion exponential

我试图在一个名为my_exp(x)的函数中使用Maclaurin系列逼近e ^ x,我相信到目前为止我所做的一切都是正确的,但我得到了不正确的近似值无论我尝试什么数字。

import math
for i in range (x):
    exp = 1 + ((x**i)/math.factorial(i))
print(exp)

例如,每当我尝试my_exp(12)时,我会得到18614.926233766233而不是162754.79141900392 帮助

2 个答案:

答案 0 :(得分:3)

为了累积系列的条款,您需要将作业替换为exp,例如:

exp = exp + ((x**i)/math.factorial(i))

答案 1 :(得分:1)

你的问题是e ^ x系列是一个无限系列,因此只对该系列的前x项进行求和是没有意义的。

def myexp(x):
  e=0
  for i in range(0,100): #Sum the first 100 terms of the series
    e=e+(x**i)/math.factorial(i)
  return e

您还可以定义结果的精确度并获得更好的解决方案。

def myexp(x):
    e=0
    pres=0.0001
    s=1
    i=1
    while s>pres:
        e=e+s
        s=(x**i)/math.factorial(i)
        i=i+1
    return e