计算机数值分析：复化辛普生公式（python实现

这里并不能够说完全实现了复化辛普生公式，因为这里面涉及到具体的函数，而我们需要事先知道函数表达式，才能够求出来。
故而这里是以书上的一道例题来写
例题以及书上结果如下：

这里我也并没有按照流程图来写（不过大同小异），而是按照表达式来写
流程图和表达式如下


代码如下：

# coding=gbk;
#因为使用复化辛普生公式会涉及到函数的具体形式
# 所以这里我就暂且令 f(x)=sin(x)/x
import math;
def compute_fx(temp):  #用来计算函数值
    if temp==0:
        return 1.0;
    return math.sin(temp)/temp;

def fuhua_simpson(a1,b1,h1,n1):  #复化辛普生
    S=compute_fx(b1)-compute_fx(a1);
    account=0;
    x=a1; 
    while account

运行结果：

遇事不决，可问春风