![351d14034f019d3cd45a0c5b3c81b6a8.png](https://i-blog.csdnimg.cn/blog_migrate/6627733dae52792274558d3505af8458.jpeg)
我上第一节数值分析课时,老师着重解释着数学原理,对于如何实现这些算法只是留给我们一些伪代码。这些伪代码由基础运算和各种for loop
组成。我想一段好的数值分析代码有无数种替代for loop
的方法,更简洁和易读。
Sum
求和符号频繁的出现在各种公式里面。比如 Composite Simpson Rule:
我之前看到求和符号的第一反应是这儿又得用 for loop
了。例如
sum=0
for i in range(1,n//2+1):
sum = sum+f(a+(2*i-1)*h)
sum
这是一个有点冗长,不清楚的写法。
Map
我们可以用np.sum
和map
函数来实现。
import numpy as np
np.sum(np.array(list(map(lambd