定積分c語言
⑴ c語言 求定積分的通用函數
對於一重定積分來說其求解可以使用梯形法進行求解，計算公式如下所示：
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
//功能：返回f(x)在積分區間[a,b]的值
//參數：FunCallBack指向用於計算f(x)的函數
//a積分區間的起始值
//b積分區間的結束值
//dx橫坐標的間隔數，越小計算結果越准確
doubleCalculate(double(*FunCallBack)(doublex),
doublea,doubleb,doubledx)
{
doubledoui;
doubletotal=0;//保存最後的計算結果
for(doui=a;doui<=b;doui+=dx)
{
total+=FunCallBack(doui)*dx;
}
returntotal;
}
doublef2(doublex)
{
returnx*x;
}
doublef(doublex)
{
returnx;
}
doublef3(doublex)
{
returnx*x*x;
}
intmain()
{
doubletotal;
total=(Calculate(f,2,3,0.000001));
printf("total=%lf
",total);
total=(Calculate(f2,2,3,0.000001));
printf("total=%lf
",total);
total=(Calculate(f3,2,3,0.000001));
printf("total=%lf
",total);
return0;
}
其中，函數f，f2，f3為自行編寫的關於x的被積函數。
運行結果：
total=2.500000
total=6.333331
total=16.249991
⑵ 有關c語言用梯形法求定積分的一個程序，請幫忙修改
#include<stdio.h>
#include<math.h>
doublef(doublex)
{
returnsqrt(4-x*x);
}
main()
{
doublea,b,s1,s2,s,h,m=1E-6;
intn,i;
printf("請輸入積分下限a和積分上限b ");
scanf("%lf%lf",&a,&b);
s2=0;
n=10000;
// do
{
s1=s2;
n=2*n;
h=(b-a)/n;
i=1;
s2=0;
do
{
s=(f(a+i*h)+f(a+(i-1)*h))*h/2;
s2=s2+s;
i=i+1;
}
while(i<=n);
}
// while((fabs(s1-s2))>(m*(fabs(s2)-fabs(s1))));
printf("f(x)在[a,b]上的定積分為%.10lf",s2);
return0;
}
⑶ c語言用蒙特卡洛方法求定積分 大佬幫忙看看哪裡不對
#include<stdio.h>
#include<time.h>
#include<stdlib.h>
doublejifen(int,int,int,int);
intmain()
{
inta,b,d,n;
doublec;
srand((unsigned)time(NULL));
printf("請輸入ab ");
scanf("%d%d",&a,&b);
printf("請輸入n ");
scanf("%d",&n);
d=1+b*b;
c=jifen(a,b,d,n);
printf("積分為%lf ",c);
return0;
}
doublejifen(inta,intb,intd,intn)
{
inti,m=0;
doublex,y,result;
for(i=1;i<=n;i++){
x=(double)rand()/(double)(RAND_MAX)*(b-a)+a;
y=(double)rand()/(double)(RAND_MAX)*d;
if(y<=1+x*x)m++;
}
result=(double)m*(b-a)*d/n;
returnresult;
};
如上 使用double 在最後加入result的強制轉換
另外在x定義那裡，改為(rand/rand_max)*(b-a)+a
因為是a,b 之間的距離隨機加a，當a不為0時原來的式子是錯的。
⑷ 用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.
",s1);
printf("梯形法算得積分值：%f.
",s2);
}
程序運行結果如下：
矩形法算得積分值：0.855821
梯形法算得積分值：0.855624
由上面的比較可知，梯形法的精度要高於矩形法。
⑸ 怎麼用C語言用累積梯形面積求定積分
#include <stdio.h>
typedef double Fun(double x);
double integral(Fun f, double x1, double x2, int n)
{
double delta = (x2 - x1) / n;
double I = 0.5 * (f(x1) + f(x2))* delta;
int i = 1;
for( ; i < n; i++)
I += f(x1 + delta * i)* delta;
return I;
}
double fun(double x)
{
return x * x;
}
int main()
{
double a, t, ret;
int n;
printf("輸入積分左右端點a，t，梯形個數n\n");
scanf("%lf %lf %d", &a, &t, &n);
ret = integral(fun, a, t, n);
printf("x^2在[%lf,%lf]上的定積分為%lf\n", a, t, ret);
return 0;
}
運行結果
輸入積分左右端點a，t，梯形個數n
0 1 2000
x^2在[0.000000,1.000000]上的定積分為0.333333
⑹ 用「累計梯形面積」法求定積分求高手！！！！ （c語言）
#include "math.h"
double f(double x) /*定義被積函數*/
{
return x*x+3*x+2;
}
double integrate(double a,double b) /*定義積分公式函數*/
{double t,h;
int i,n=100;
t=0.5*(f(a)+f(b));
h=fabs(a-b)/n;
for(i=1;i<n;i++)
t+=f(a+i*h);
t=t*h;
return t;
}
void main() /*定義主函數*/
{printf("\n%f\n",integrate(0,2));}
⑺ C語言求定積分的問題。
根據梯形法求積分的原理，設間隔h= (b-a)/n，則積分近似計算公式為：
s = h/2 *[f(a)+f(a+h)] +h/2 *[f(a+h)+f(a+2h)] +...+h/2 *[f(b-h)+f(b)]
=h/2 *[f(a)+f(b)] + h*[ f(a+h) + f(a+2h) +f(a+3h) + ... +f(b -h)]
令積分s初始值為h/2 *[f(a)+f(b)] ，後面令i=1,...,n-1來迭代s = s+h*f(a+ih)。
⑻ C語言實慣用梯形法或辛普森法求解定積分的值
//梯形法求定積分
#include<stdio.h>
#include<math.h>
//定義被積函數
double func(double x){
return sin(x)*cos(x);
}
void main(){
double a,b,h,x,sum;
int i,n;
printf("Input a b and n: ");
scanf("%lf%lf%d",&a,&b,&n);
h=(b-a)/n;
x=a;
sum=(func(a)+func(b))/2;
for(i=1; i<n; i++){
x += h;
sum += func(x);
}
sum *= h;
printf("sum=%.4lf\n",sum);
}