求一段matlab求积分程序

画图程序：

x=linspace(-02,07,100);

f1=x;

f2=31055x^5-515x^4+22945x^3-04345x^2+11845x;

plot(x,f1,x,f2,'r');

首先是求交点，也就是求方程的根，如下：

f=@(x)31055x^5-515x^4+22945x^3-04345x^2+01845x;

z=fzero(f,07)

z =

06704

两个根为06704和0；

然后就是分别求积分，也就是与X轴的面积，再相减就是所围的面积

>> data=@(x)31055x^5-515x^4+22945x^3-04345x^2+11845x;

y1=quad(data,0,06704)

y1 =

02459

>> f=@(x)x;

y2=quad(f,0,06704)

y2 =

02247

>> ff=y1-y2

ff =

00212

ff即为面积！！！

1、算法原理：将sin(x)+e^x转化为用梯形法积分计算。2、程序流程图：i=1,j=3Double sum=0,x,y;int m;x=(j-i)/m;y=in=0n<m n=n+1 sum=sum+(x/2)(sin(y)+sin(y+x)+exp(y)+exp(y+x));y=y+x 输出 sin(x)+e^x在积分区间[1，3]的积分是 sum 3、变量说明：i,j控制区间的变量，求和变量sum，过程变量y，分成n等份的梯形的高变量x,控制分段数的变量m4、使用方法：确定要求的区间的大小；用循环逐步进行计算其中注意变量的条件的改变进而影响的循环条件；循环逐步进行计算5、源程序：#include<iostream>#include<cmath>using namespace std;void main(){const double i=1, j=3;</p><p>double sum=0,x,y;</p><p>int m,n=0;</p><p>cout<<"请输入分段数:";</p><p>cin>>m; </p><p>x=(j-i)/m;</p><p>y=i;</p><p>for(n=0;n<m;n++)</p><p>{</p><p>sum=sum+(x/2)(sin(y)+sin(y+x)+exp(y)+exp(y+x));</p><p>y=y+x;</p><p>}cout<<"sin(x)+e^x在积分区间[1，3]的积分是"<<sum<<endl;} 6、程序运行结果：

#include <stdioh>

float integral(float a)

{

return (aaaa-(a-1)(a-1)(a-1)(a-1))/4;

}

int main()

{

float i,sum = 0;

for(i=10;i>=1;i--)

{

sum+=integral(i);

}

printf("sum =%f",sum);

getchar();

return 0;

}

记得采纳

呵呵，给你一个完整的，我已经在VC++60上编译运行确认了：

问题描述：

在数值积分法中，利用求单位圆的1/4的面积来求得Pi/4从而得到Pi。单位圆的1/4面积是一个扇形，它是边长为1单位正方形的一部分。只要能求出扇形面积S1在正方形面积S中占的比例K=S1/S就立即能得到S1，从而得到Pi的值。怎样求出扇形面积在正方形面积中占的比例K呢？一个办法是在正方形中随机投入很多点，使所投的点落在正方形中每一个位置的机会相等看其中有多少个点落在扇形内。将落在扇形内的点数m与所投点的总数n的比m/n作为k的近似值。P落在扇形内的充要条件是x^2+y^2<=1。

程序描述：

#include<iostream>

#include<cmath>

#include<ctime>

#define COUNT 500000 //循环取样次数

using namespace std;

bool InCircle(double x,double y)//是否在1/4圆范围之内

{

if((xx+yy)<=1)return true;

return false;

}

void main()

{

double x,y;

int num=0;

int i;

srand((unsigned)time(NULL));

for(i=0;i<COUNT;i++)

{

x=rand()10/RAND_MAX;

y=rand()10/RAND_MAX;

if(InCircle(x,y)) num++;

}

cout<<"PI:"<<(num40)/COUNT<<endl;

}

以上就是关于求一段matlab求积分程序全部的内容，包括:求一段matlab求积分程序、4． 编写使用梯形法计算定积分的程序。被积函数可取sin(x)+ex，积分区间[1，3]。、用c语言编一个程序求定积分 ∫ab x3（其中ab为上下限，a=10，b=1，x3表示x的3次方）等相关内容解答，如果想了解更多相关内容，可以关注我们，你们的支持是我们更新的动力！