這是辛普森積分法。給你寫了fun_1( ),fun_2(),請(qǐng)自己添加另外幾個(gè)被積函數(shù)。調(diào)用方法 t=fsimp(a,b,eps,fun_i);a,b --上下限,eps -- 迭代精度要求。
創(chuàng)新互聯(lián)公司是一家專注網(wǎng)站建設(shè)、網(wǎng)絡(luò)營(yíng)銷策劃、小程序設(shè)計(jì)、電子商務(wù)建設(shè)、網(wǎng)絡(luò)推廣、移動(dòng)互聯(lián)開發(fā)、研究、服務(wù)為一體的技術(shù)型公司。公司成立十載以來(lái),已經(jīng)為1000+混凝土攪拌站各業(yè)的企業(yè)公司提供互聯(lián)網(wǎng)服務(wù)?,F(xiàn)在,服務(wù)的1000+客戶與我們一路同行,見證我們的成長(zhǎng);未來(lái),我們一起分享成功的喜悅。
例:求函數(shù)f(x)=x*x+2*x+1在【0,2】上的定積分。
(x*x)在(0,1)上定積分為%lf\n,fun(0,1,1000000));//區(qū)間數(shù)自己設(shè)越大值越精確 } 結(jié)果:數(shù)學(xué)分析:f(x)=x^2=x*x;定積分:x*x*x/3+c(常數(shù))在區(qū)間(0,1)上定積分:1/3=0.333333 結(jié)果正確。
void main(){ double sum=0;int i;for(i=0; iN; i++){ sum += sin((double)(i)/N)/N;} printf(%lf\n%lf,sum,1-cos(1));} N后面的0有點(diǎn)多了,不過這個(gè)數(shù)剛好能精確到小數(shù)點(diǎn)后6位。
1、這是辛普森積分法。給你寫了fun_1( ),fun_2(),請(qǐng)自己添加另外幾個(gè)被積函數(shù)。調(diào)用方法 t=fsimp(a,b,eps,fun_i);a,b --上下限,eps -- 迭代精度要求。
2、H=(B-A)/N這一句不應(yīng)該放在變量聲明中,因?yàn)閯偠x時(shí),A,B,N沒有賦值,是隨機(jī)值 應(yīng)該將 float A,B,N,H=(B-A)/N,sum=0;改為:float A,B,N,H,sum=0;將H=(B-A)/N;放在scanf (%f,&N); 后面。
3、積分 由定義可知:求函數(shù)f(x)的不定積分,就是要求出f(x)的所有的原函數(shù),由原函數(shù)的性質(zhì)可知,只要求出函數(shù)f(x)的一個(gè)原函數(shù),再加上任意的常數(shù)C,就得到函數(shù)f(x)的不定積分。
4、http://blog.itpub點(diǎn)虐 /14139916/viewspace-497879/ 求定積分的近似值常有矩形法與梯形法,其實(shí)質(zhì)都是面積求和。矩形法是把所要求的面積垂直x軸分成n個(gè)小矩形,然后把這n個(gè)小矩形的面積相加,即為所求的定積分的值。
5、i=0,b=10,n=1000,得到result=36833,與精確值比較接近。因?yàn)榫匦畏ㄇ蠓e分本來(lái)就不太準(zhǔn) 另外,樓上的同學(xué)不要誤導(dǎo)人家,積分的起始值可不一定是0,s=n*((n*i)*(n*i)+3)+s;明顯把i默認(rèn)當(dāng)做0了。
6、這個(gè)你用手工積分就行了,當(dāng)然常數(shù)的代入計(jì)算還是可以用計(jì)算器的。
基本是這樣的,用梯形發(fā)求定積分,對(duì)應(yīng)于一個(gè)積分式就要有一段程序,不過你可以改變程序的一小部分來(lái)改變你所要求的積分式。
求定積分 f(x)dx, x=a 到 b 一般用 數(shù)值方法計(jì)算 計(jì)算 f(x) 曲線與x軸之間的面積。最普通的用梯形法,并且用2分法 加密。直到滿足精度。