Matlab作业（一）
作业要求：用梯形法或者辛普森法数值积分，分别用Matlab和c语言实现。
C语言
1. 程序代码：
#include #include double fun(double x)
{
return x*x;
}
double definfresult1(double (*pfun)(double),double a,double b,double eps)
{
int n=1;
--
double h,k,tn,tn1,fh,fh1=0;
double fa=pfun(a);
double fb=pfun(b);
tn=(b-a)*(fa+fb)/2;
do {
for(k=0,fh1=0;kh=(b-a)/(n);fh=pfun(a + (2*k + 1)*(b-a)/(2*n));
fh=fh+fh1;
fh1=fh;
}
tn1=tn;
tn=(tn1+fh*h)/2;
--
n=2*n;
} while(fabs(tn- tn1) >= eps);
return tn;
}
int main()
{
double a,b,eps,definfresult;
printf(\"积分下限a=\");
scanf(\"%lf\
printf(\"积分上限b=\");
scanf(\"%lf\
printf(\"精度eps=\");
scanf(\"%lf\
--
definfresult=definfresult1(fun,a,b,eps);
printf(\"\\n计算结果=%.7lf\\n\
}
2.运行结果：
MatLab
1. 程序代码：
function y=fun1(x)
y=x*x;
--
function result=definf1(fhandle, a, b, eps)
fa=feval(fhandle, a);
fb=feval(fhandle, b);
tn=(b-a)*(fa+fb)/2;
tn1=0;
n=1;
fh1=0;
while abs(tn- tn1) > eps
fh1=0;
for i=0:n-1
h=(b-a)/n;
fh=feval(fhandle,a + (2*i + 1)*(b-a)/(2*n) );
fh=fh+fh1;
--
fh1=fh;
end;
tn1=tn;
tn=(tn1+fh*h)/2;
n=2*n;
end
result=tn;
>> result=definf1(@fun1, 0, 10,0.001 )
result =
333.3335
2.运行结果
--
--
THANKS
致力为企业和个人提供合同协议，策划案计划书，学习课件等等
打造全网一站式需求
欢迎您的下载，资料仅供参考
--