怎么在C語言中使用矩形法求定積分?針對這個問題,這篇文章詳細介紹了相對應的分析和解答,希望可以幫助更多想解決這個問題的小伙伴找到更簡單易行的方法。
要求:
寫一個用矩形法求定積分的通用函數,分別求:sin(x),cos(x),e^x 。
分析:
矩形法,學過高等數學就知道化曲為直的思想。將定積分化為多個函數連續的和?;舅枷胧菍^間[a,b]化成n等分,當n越大的時候結果越準確。圖形化成一小塊一小塊的矩形。底邊長都為(b-a)/n.高為每個等分點的函數值。然后將每個矩形的面積相加即為所求。
如:
y=x;
可以通過矩形的方法來無限逼近定積分的求解,如下:
因為被分成n等分,就可以認為每一等分是一個矩形,那么每一矩形的面積為: 每一個矩形面積為:***Sn=f(x)(b-a)/n 總面積為:****S=S1+S2+…+Sn
#include <stdio.h> #include <stdlib.h> #include <math.h>//引入sin x,cos x,e^x的庫 //使用指向函數的指針變量來復用一個通用函數 int main() { int i,j,k,n,m; float res; float (*fun)(float);//定義指向函數的指針變量 float integral (float a,float b,float (*fun)(float),int n);//n為將積分區間(b-a)分成n等分,當n的值越大的時候結果越精確,數學上的定義是取無窮大 float fsin(float);//sin x函數的聲明 float fcos(float);//cos x函數的聲明 float fexp(float);//e^x 函數的聲明 printf ("請輸入積分的下限\n"); scanf ("%d",&m); printf ("請輸入積分的上限\n"); scanf ("%d",&n); printf("請輸入你要計算的函數的具體函數\n"); printf("1.sin(x) 2.cos(x) 3.e^x\n"); scanf ("%d",&i); switch(i) { case 1: fun=fsin;//函數地址(入口)交給指針變量,靈活性強 break; case 2: fun=fcos; break; case 3: fun=fexp; } res=(fun)(2.00); //printf("xxxx=%f\n",res); printf ("計算的結果為\n"); res=integral(m,n,fun,200000); printf("res=%f",res); return 0; } float fsin(float x) { // printf("fsinx=%f\n",x); return sin(x); } float fcos(float x) { //printf("fcosx=%f\n",x); return cos(x); } float fexp(float x) { //getchar(); //printf("fexp=%f\n",x); return exp(x); } float integral(float a,float b,float (*fun)(float),int n) { //矩形法計算,定積分轉換為連續求和的形式 int i=0; float x=a,s=0; float h=(b-a)/n; for (i=1;i<=n;i++) { x=x+h; s=s+((*fun)(x)*h); } return s; }
結果:
關于怎么在C語言中使用矩形法求定積分問題的解答就分享到這里了,希望以上內容可以對大家有一定的幫助,如果你還有很多疑惑沒有解開,可以關注億速云行業資訊頻道了解更多相關知識。
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。