这是辛普森积分法。给你写了fun_1( ），fun_2(），请自己添加另外几个被积函数。调用方法 t=fsimp(a，b，eps，fun_i)；a，b --上下限，eps -- 迭代精度要求。
成都创新互联公司主要从事做网站、网站制作、网页设计、企业做网站、公司建网站等业务。立足成都服务新城,十余年网站建设经验,价格优惠、服务专业,欢迎来电咨询建站服务:18982081108
对于一重定积分来说其求解可以使用梯形法进行求解，计算公式如下所示：其中，f(x)为被积函数，为横坐标的两点间的间隔，越小，则计算出的结果越精确。
例：求函数f(x)=x*x+2*x+1在【0，2】上的定积分。
(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 结果正确。
基本是这样的，用梯形发求定积分，对应于一个积分式就要有一段程序，不过你可以改变程序的一小部分来改变你所要求的积分式。
求定积分的近似值常有矩形法与梯形法，其实质都是面积求和。矩形法是把所要求的面积垂直x轴分成n个小矩形，然后把这n个小矩形的面积相加，即为所求的定积分的值。
首先解决怎么算，计算机肯定不会积分，所以我开始想用sinx的泰勒展开式，然后选3-4次作为近似，然后积分。听你说梯形法，是数值计算的内容，刚好这学期在学，就把我调试的程序发一个给你吧这是romberg算法，把a 换为0，b换为pi就好了吧。
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 结果正确。