Matlab实验报告四(矩形法梯形法抛物线法求定积分).doc

《Matlab实验报告四(矩形法梯形法抛物线法求定积分).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数学与信息科学系实验报告实验名称定积分的近似计算所属课程数学软件与实验实验类型综合型实验专业信息与计算科学班级学号姓名指导教师一、实验概述【实验目的】【实验原理】【实验环境】MatlabR2010b二、实验内容分别用中点矩形法、梯形法、抛物线法计算，其中的值分别为.1.分析问题2.问题求解矩形法：functionf=zuoye(fun,a,b,n)h=(b-a)/n;x=a:h:b;y=x;fori=2:n+1y(i)=fun((x(i)+x(i-1))/2);endf=h*sum(y(1:end))梯形法：functionf=zuoye1(fun,a,b,n)h=(b-a)/n

2、;x=a:h:b;y=x;fori=1:n+1y(i)=fun(x(i));endf=h*(sum(y)-1/2*y(1)-1/2*y(end));抛物线法：functionf=zuoye3(fun,a,b,n)h=(b-a)/(2*n);x=a:h:b;y=x;fori=1:2*n+1y(i)=fun(x(i));ends1=0;s2=0;fori=1:ns1=s1+4*y(2*i);endfori=2:ns2=s2+2*y(2*i-1);ends=y(1)+y(end)+s1+s2;f=(b-a)/(6*n)*s;x=0.1的方程：functionf=f1(x)f=sqrt(

3、1+0.1^2*cos(x).^2);x=0.5的方程：functionf=f1(x)f=sqrt(1+0.5^2*cos(x).^2);x=0.9的方程：functionf=f1(x)f=sqrt(1+0.9^2*cos(x).^2);3.结果f=zuoye(@f1,0,2*pi,100)f=6.2989f=zuoye1(@f1,0,2*pi,100)f=6.2989f=zuoye3(@f1,0,2*pi,100)x=6.29894.结论及分析三、实验小结