请编写一个程序horner. py,实现如下功能:编写函数evaluate(x,a),计算多项式a(x)的值,其中,a(x)的系数为数组a[]中的各元素。 a(x)=a0+a1x1+a2x2+⋯+an−2xn−2+an−1xn−1 使用霍纳法,一种有效的计算方法是使用如下建议的括号表达式: a0+x(a1+x(a2+⋯+x(an−2+xan−1)⋯)) 请编写一个函数exp(),调用函数evaluate()以求解( ˣ eˣ 的近似值,使用泰勒级数展开式的前n项: 从命令行接收一个参数x,并把计算结果与math.exp(x)的结果进行比较。
时间: 2023-05-21 11:02:17 浏览: 121
运用Horner规则计算多项式的值
以下是horner.py的代码:
```python
def evaluate(x, a):
n = len(a)
result = a[n-1]
for i in range(n-2, -1, -1):
result = result * x + a[i]
return result
def exp(x, n):
a = [1.0]
for i in range(1, n+1):
a.append(a[i-1] / i * x)
return evaluate(1.0, a)
x = float(input("请输入x的值:"))
n = int(input("请输入泰勒级数展开式的项数n:"))
result = exp(x, n)
print("exp({})的近似值为:{}".format(x, result))
print("math.exp({})的值为:{}".format(x, math.exp(x)))
```
如果你输入x=1,n=10,那么输出结果应该是:
```
请输入x的值:1
请输入泰勒级数展开式的项数n:10
exp(1.0)的近似值为:2.7182818011463845
math.exp(1.0)的值为:2.718281828459045
```
注意,由于泰勒级数展开式只是e^x的近似值,因此当n不够大时,计算结果可能与math.exp(x)的值有一定误差。
阅读全文