基本是这样的，用梯形发求定积分，对应于一个积分式就要有一段程序，不过你可以改变程序的一小部分来改变你所要求的积分式。
目前创新互联公司已为数千家的企业提供了网站建设、域名、网页空间、绵阳服务器托管、企业网站设计、台安网站维护等服务，公司将坚持客户导向、应用为本的策略，正道将秉承"和谐、参与、激情"的文化，与客户和合作伙伴齐心协力一起成长，共同发展。
求定积分的近似值常有矩形法与梯形法，其实质都是面积求和。矩形法是把所要求的面积垂直x轴分成n个小矩形，然后把这n个小矩形的面积相加，即为所求的定积分的值。
用梯形法估算，再用辛普森法。fsimpf 积分函数 a，b 积分下上限，eps 精度。
首先解决怎么算，计算机肯定不会积分，所以我开始想用sinx的泰勒展开式，然后选3-4次作为近似，然后积分。听你说梯形法，是数值计算的内容，刚好这学期在学，就把我调试的程序发一个给你吧这是romberg算法，把a 换为0，b换为pi就好了吧。
求定积分 f(x)dx， x=a 到 b 一般用 数值方法计算 计算 f(x) 曲线与x轴之间的面积。最普通的用梯形法，并且用2分法 加密。直到满足精度。
{ double sum=0；int i；for(i=0； iN； i++){ sum += sin((double)(i)/N)/N；} printf(%lf\n%lf，sum，1-cos(1))；} N后面的0有点多了，不过这个数刚好能精确到小数点后6位。
1、这是辛普森积分法。给你写了fun_1( ），fun_2(），请自己添加另外几个被积函数。调用方法 t=fsimp(a，b，eps，fun_i)；a，b --上下限，eps -- 迭代精度要求。
2、例：求函数f(x)=x*x+2*x+1在【0，2】上的定积分。
3、(x*x)在(0，1)上定积分为%lf\n，fun(0，1，1000000))；//区间数自己设越大值越精确 } 结果：数学分析：f(x)=x^2=x*x；定积分：x*x*x/3+c(常数)在区间(0，1)上定积分：1/3=0.333333 结果正确。
这是辛普森积分法。给你写了fun_1( ），fun_2(），请自己添加另外几个被积函数。调用方法 t=fsimp(a，b，eps，fun_i)；a，b --上下限，eps -- 迭代精度要求。
对于一重定积分来说其求解可以使用梯形法进行求解，计算公式如下所示：其中，f(x)为被积函数，为横坐标的两点间的间隔，越小，则计算出的结果越精确。
例：求函数f(x)=x*x+2*x+1在【0，2】上的定积分。
