c語言 求定積分的通用函數

對于一重定積分來說其求解可以使用梯形法進行求解，計算公式如下所示：

成都創新互聯是一家專注于成都做網站、成都網站設計、成都外貿網站建設與策劃設計,仁和網站建設哪家好?成都創新互聯做網站,專注于網站建設十年,網設計領域的專業建站公司;建站業務涵蓋:仁和等地區。仁和做網站價格咨詢:18980820575

其中，f(x)為被積函數，為橫坐標的兩點間的間隔，越小，則計算出的結果越精確。

對于求解此類問題可以使用C語言中的回調函數編寫通用的計算函數，代碼如下：

#include?stdio.h

#include?stdlib.h

#includemath.h

//功能：返回f(x)在積分區間[a,b]的值

//參數：FunCallBack?指向用于計算f(x)的函數

//??????a??積分區間的起始值

//??????b??積分區間的結束值

//??????dx?橫坐標的間隔數，越小計算結果越準確

double?Calculate(double?(*FunCallBack)(double?x),

double?a,double?b,double?dx)

{

double?doui;

double?total?=?0;????????//保存最后的計算結果

for?(doui?=?a;?doui?=?b;?doui?+=?dx)

{

total?+=?FunCallBack(doui)*dx;

}

return?total;

}

double?f2(double?x)

{

return?x*x;

}

double?f(double?x)

{

return?x;

}

double?f3(double?x)

{

return?x*x*x?;

}

int?main()

{

double?total;

total?=?(Calculate(f,?2,?3,?0.000001));

printf("total?=?%lf\n",?total);

total?=?(Calculate(f2,?2,?3,?0.000001));

printf("total?=?%lf\n",?total);

total?=?(Calculate(f3,?2,?3,?0.000001));

printf("total?=?%lf\n",?total);

return?0?;

}

其中，函數f，f2，f3為自行編寫的關于x的被積函數。

運行結果：

total?=?2.500000

total?=?6.333331

total?=?16.249991

用C語言求積分

基本是這樣的，用梯形發求定積分，對應于一個積分式就要有一段程序，不過你可以改變程序的一小部分來改變你所要求的積分式。

以c為例：求f(x)=xsinx從1到2的積分

#include math.h

float integral(float(*fun)(float x),float a,float b,int,n)

{float 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;/*返回積分值*/

}

float f(float x)

{return(x*sinx) /*修改此處可以改變被積函數*/

}

main()

{float y;

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

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

}

用C語言求定積分

實際問題描述：

求定積分近似值

程序代碼如下：

#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語言積分函數

積分分為兩種，數值積分，公式積分。

公式積分：部分函數可以直接用公式求得其不定積分函數。C語言中可以直接用積分公式寫出其積分函數。

數值積分：按照積分的定義，設置積分范圍的步長，用梯形面積累加求得其積分。

以【f(x)=x*sin(x) 從1到2的積分】為例：

#include?math.h

#include?stdio.h

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;

}

當前文章：c語言加分函數 c語言分數加減乘除
分享URL：http://www.yijiale78.com/article28/dopsdjp.html

成都網站建設公司_創新互聯，為您提供用戶體驗、企業建站、建站公司、品牌網站建設、靜態網站、域名注冊

廣告

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