【正文】
0 for j=1:n/2 s=s+f(2*j) end t=0 for j=1:(n/21) t=t+f(2*j1) end S(i)=1/3/n*(f(1)+4*s+2*t+f(n+1))end五.實(shí)驗(yàn)內(nèi)容復(fù)化梯形公式和復(fù)化辛普森公式的引入復(fù)化梯形公式:;復(fù)化辛普森公式:;根據(jù)題意和復(fù)化梯形公式、復(fù)化辛普森公式的原理編輯程序求解代碼如下:Matlab代碼clcs=quad(39。p2。既糾結(jié)了自己,又打擾了別人。2. 若不是心寬似海,哪有人生風(fēng)平浪靜。 end for j=1:(k1) s2=s2+(1/6)*(x(j+1)x(j))*(y(j)+y(j+1)+4*sz(j))。二、實(shí)驗(yàn)內(nèi)容:對(duì)于函數(shù),試?yán)孟卤碛?jì)算積分。 sz=sin(z)./z。由于以上的算法只算了結(jié)點(diǎn)個(gè)數(shù)為9的情況,只能橫向比較兩公式的精確程度,而不能分別比較兩公式隨節(jié)點(diǎn)個(gè)數(shù)變化精度的變化,故而將以上程序重新編(以上程序?yàn)樽罱K