matlab的數值求解實驗報告,matlab計算方法實驗報告5(數值積分)

計算方法實驗報告(5)

學生姓名楊賢邦學号指導教師吳明芬實驗時間2014.4.16地點綜合實驗大樓203

實驗題目數值積分方法

實驗目的●利用複化梯形、辛普森公式和龍貝格數值積分公式計算定積分的

近似植。

實驗内容●梯形、辛普森、柯特斯法及其Matlab實作；

●變步長的梯形、辛普森、柯特斯法及其Matlab實作。

●題目由同學從學習材料中任意選兩題

算法分析梯形：function y=jifeng_tixing(a,b,n,fun)

fa=feval(fun,a);

fb=feval(fun,b);

s=0;

h=(b-a)/n;

for k=1:n-1

xk=a+k*h;

s=feval(fun,xk)+s;

end

y=(h/2)*(fa+fb+2*s);

辛普生：function y=jifeng_xingpu(a,b,n,fun) fa=feval(fun,a);

fb=feval(fun,b);

h=(b-a)/n;

s=0;

s2=feval(fun,a+0.5*h);

for k=1:n-1

xk=a+k*h;

s=feval(fun,xk)+s;

s2=feval(fun,xk+(h/2))+s2;

end