c語言求定積分

問題就是出在數(shù)據(jù)類型上的選用上，precision=0.0000001時(shí)已經(jīng)超過了float的數(shù)據(jù)范圍，所以導(dǎo)致數(shù)據(jù)截?cái)嗪髉recision=0.000000，從而程序在計(jì)算積分時(shí)可能陷入死循環(huán)，應(yīng)該采用double型數(shù)據(jù)類型。其實(shí)不推薦樓主用如此多的define語句，程序的可讀性和風(fēng)格應(yīng)該重于編程員的勞動度。。。

我們提供的服務(wù)有：網(wǎng)站設(shè)計(jì)制作、成都做網(wǎng)站、微信公眾號開發(fā)、網(wǎng)站優(yōu)化、網(wǎng)站認(rèn)證、富源ssl等。為成百上千企事業(yè)單位解決了網(wǎng)站和推廣的問題。提供周到的售前咨詢和貼心的售后服務(wù)，是有科學(xué)管理、有技術(shù)的富源網(wǎng)站制作公司

還有樓主對自然對數(shù)e的define也已經(jīng)超過了計(jì)算機(jī)的可識別范圍。。您那樣精確的定義e并不會在結(jié)果上獲得更加精確地結(jié)果，其實(shí)反倒會起到相反的作用，要知道與其用一個(gè)這樣可寬穗能導(dǎo)致內(nèi)存出錯(cuò)以及必定會導(dǎo)致數(shù)據(jù)截?cái)嗟淖兞縼韺?shí)現(xiàn)精度的提高遠(yuǎn)遠(yuǎn)不如采用一個(gè)更精確的積分算法，而且c語言提供了自然數(shù)e為底的指數(shù)函數(shù)~而且貌似您的積分算法是不準(zhǔn)確的，梯形積分的定義并非如此，其再兩端的函數(shù)值應(yīng)該只取1/2.希望您多加細(xì)心~

如果不介意的話，就是你的precision應(yīng)該改為step~這樣會能更加準(zhǔn)備的表達(dá)了這個(gè)變量的作用，在你的程序中precision變量其實(shí)是積分步長~在數(shù)值計(jì)算方法中積分精度的控制往往不是通過細(xì)化步長來表達(dá)，而是通過后一個(gè)積分值-前一個(gè)積分值precision

這樣來實(shí)缺仔現(xiàn)精度控制~呵呵伏巧汪

用C語言編寫一個(gè)求定積分的程序

這是辛普森積分法。

給你寫了fun_1( ），fun_2(），請自己添加另外幾個(gè)被積函肆神罩?jǐn)?shù)。

調(diào)用方法 t=fsimp(a,b,eps,fun_i);

a,b --上下限，eps -- 迭代精度要瞎禪求。

#includestdio.h

#includestdlib.h

#include math.h

double fun_1(double x)

{

return 1.0 + x ;

}

double fun_2(double x)

{

return 2.0 * x + 3.0 ;

}

double fsimp(double a,double b,double eps, double (*P)(double))

{

int n,k;

double h,t1,t2,s1,s2,ep,p,x;

n=1; h=b-a;

t1=h*(P(a)+P(b))/2.0;

s1=t1;

ep=eps+1.0;

while (ep=eps)

{

p=0.0;

for (k=0;k=n-1;k++)

{

x=a+(k+0.5)*h;

p=p+P(x);

}

t2=(t1+h*p)/裂鬧2.0;

s2=(4.0*t2-t1)/3.0;

ep=fabs(s2-s1);

t1=t2; s1=s2; n=n+n; h=h/2.0;

}

return(s2);

}

void main()

{

double a,b,eps,t;

a=0.0; b=3.141592653589793238; eps=0.0000001;

// a definite integral by Simpson Method.

t=fsimp(a,b,eps,fun_1);

printf("%g\n",t);

t=fsimp(a,b,eps,fun_2);

printf("%g\n",t);

// ...

printf("\n Press any key to quit...");

getch();

}

用C語言求定積分

實(shí)際問題描述：

求定積分近似值

程序代碼如下：

#include

#include

void main()

{

int i,n=1000;

float a,b,h,t1,t2,s1,s2,x;

printf("請輸入積分限a,b:");

scanf("%f,%f",a,b);

h=(b-a)/n;

for(s1=0,s2=0,i=1;i=n;i++)

{

猛譽(yù)純 x=a+(i-1)*h;

t1=(float)exp(-x*x/2);t2(float)=exp(-(x+h)*(x+h)/2);

s1=s1+t1*h;??????? /*矩形面積累加*/

s2=s2+(t1+t2)*h/2;??????? /*梯形面積累加*/

}

printf("矩形法算得積分值：%f.\n",s1);

printf("梯形法算得積分值：%f.\n",s2);

}

程序運(yùn)行結(jié)果如下：

虛輪?? 矩形法枝咐算得積分值：0.855821

梯形法算得積分值：0.855624

由上面的比較可知，梯形法的精度要高于矩形法。

用c語言求定積分

#include

#include

double integral(double(*fun)(double x),double a,double b,int,n){

double s,h,y;

int i;

s=(fun(a)+fun(b))/2;

h=(b-a)/n; /*積分步長*/

for(i=1;in;i++)

s=s+fun(a+i*h);

y=s*h;

return y;/*返回積分值*/

}

double f(double x){

return(x*sinx) /*修改此耐宴處可以改變被積函數(shù)*/

}

int main(){

double y;

y=integral(f,1.0,2.0,150);/*修改此處昌沒銀可以改變積分上下限和步數(shù),步長=(上限察鍵-下限)/步數(shù)*/

printf("y=%f\n",y);

return 0;

}

int main()