用C語言求定積分

實際問題描述：

成都創新互聯公司-專業網站定制、快速模板網站建設、高性價比民勤網站開發、企業建站全套包干低至880元,成熟完善的模板庫,直接使用。一站式民勤網站制作公司更省心,省錢,快速模板網站建設找我們，業務覆蓋民勤地區。費用合理售后完善，10年實體公司更值得信賴。

求定積分近似值

程序代碼如下：

#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++)

{

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);

}

程序運行結果如下：

矩形法算得積分值：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) /*修改此處可以改變被積函數*/

}

int main(){

double y;

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

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

return 0;

}

int main()

用C語言編寫一個求定積分的程序，用它分別求以下5個函數的定積分：

#includestdio.h

#includemath.h

float f1(float x)

{

return(1.0+x);

}

float f2(float x)

{

return(2.0*x+3.0);

}

float f3(float x)

{

return(exp(x)+1);

}

float f4(float x)

{

return(pow(1+x,2));

}

float f5(float x)

{

return(pow(x,3));

}

float fsimp(float a,float b,float (*p)(float))

{

float c,s;

c=(a+b)/2;

s=(b-a)/6*(p(a)+4*p(c)+p(b));

return s;

}

int main()

{

float a,b;

printf("請輸入積分下限a的值：");

scanf("%f",a);

printf("請輸入積分上限b的值：");

scanf("%f",b);

printf("%f\n",fsimp(a,b,f1));

printf("%f\n",fsimp(a,b,f2));

printf("%f\n",fsimp(a,b,f3));

printf("%f\n",fsimp(a,b,f4));

printf("%f\n",fsimp(a,b,f5));

}

分享標題：c語言求函數的定積分,C語言求定積分
URL鏈接：http://www.yijiale78.com/article40/dsiopeo.html

成都網站建設公司_創新互聯，為您提供虛擬主機、品牌網站設計、營銷型網站建設、網頁設計公司、手機網站建設、網站設計公司

廣告

聲明：本網站發布的內容（圖片、視頻和文字）以用戶投稿、用戶轉載內容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網站立場，如需處理請聯系客服。電話：028-86922220；郵箱：631063699@qq.com。內容未經允許不得轉載，或轉載時需注明來源： 創新互聯