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

数学与信息科学系实验报告实验名称定积分的近似计算所属课程数学软件与实验实验类型综合型实验专业信息与计算科学班级学号姓名指导教师1一、实验概述【实验目的】【实验原理】【实验环境】MatlabR2010b二、实验内容分别用中点矩形法、梯形法、抛物线法计算22201cossd，其中的值分别为0.1,0.5,0.9.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;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;2s2=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(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.结论及分析三、实验小结