nes(1,n-2) 1]; I=h/2*sum(a.*y); end %计算并作图 n=1:100; t1=ones(1,100); for k=1:100 t1(k)=t1(k)*tq(k); end f=inline( 'x./(4+x.^2)' ); I=quadl(f,0,1); t2=I-t1; plot(n,t2, '*k' ,n,t2, '-y' ) (2) %构造复合辛普森积分公式 function I=simpsion(x,y) m=length(x); n=length(y); if m~=n error; end if rem(n-1,2)~=0 I=tquad(x,y); return ; end ; N=(n-1)/2; h=(x(n)-x(1))/N; a=zeros(1,n); for k=1:N a(2*k-1)=a(2*k-1)+1; a(2*k)=a(2*k)+4; a(2*k+1)=a(2*k+1)+1; end I=h/6*sum(a.*y); End %分别计算积分 format long x=0:0.1:1; y=x./(4+x.^2); i simosion = simpsion (x,y) itquad=tquad(x,y) 三、结果通过 Matlab 程序运行结果如下: 1.the population in1930 is169.649000 the population in1965 is191.767359 the population in2010 is171.351000 由于 lagrange 插值不能准确估计范围外的数值,因此 1930 年和 2010 年的误差较大。 2.ais0.681361 bis0.230620 3.(1)I1 =0.111463380815167 I=0.111571775612031
全文预览
