用Python实现一个数值计算方面炫技性的程序案例,要求用尽可能短的代码实现一种数值积分方法
时间: 2024-02-06 18:12:23 浏览: 82
Python实现数值积分方式
5星 · 资源好评率100%
好的,我会回答你的问题。这里提供一个基于Simpson方法的数值积分程序案例,将函数f(x)从a到b进行积分:
```python
def simpson_integrate(f, a, b, n=4):
'''用Simpson方法对函数进行数值积分'''
h = (b-a)/n
x = [a+i*h for i in range(n+1)]
y = [f(x[i]) for i in range(n+1)]
s = y[0] + y[-1] + 4*sum(y[1:-1:2]) + 2*sum(y[2:-1:2])
return h/3 * s
```
这个程序使用了Python内置的列表操作和简化的Simpson方法公式来进行数值积分。你可以调用这个函数来计算函数在任何区间上的数值积分:
```python
def f(x):
return x**2 + 2*x + 1
a = 0
b = 10
approx = simpson_integrate(f, a, b, 100)
print(f"The numerical integral of f(x) from {a} to {b} is approximately {approx:.6f}")
```
我希望这个程序能够满足你的需求!
阅读全文