亲宝软件园·资讯

展开

C语言求解定积分 C语言求解定积分的方法

xiaocaidayong 人气:2

题目要求:

求下面函数的定积分:

思路:

求一个函数的定积分,其实就是求它的面积,如对函数

求积分。

也就是要求出绿色部分的面积,如下:

我们可以通过矩形的方法来无限逼近定积分的求解,如下:

因为被分成n等分,就可以认为每一等分是一个矩形,那么每一矩形的面积为:

每一个矩形面积为:***Sn=f(x)(b-a)/n
总面积为:****S=S1+S2+…+Sn

通过这样的一个思路就可以完成,定积分的求解。
这样这三个定积分的求解代码,如下:

#include <stdio.h>
#include <math.h>

void main( )
{
 float f1(float x);
 float f2(float x);
 float f3(float x);
 float s1 = 0;
 float s2 = 0;
 float s3 = 0;
 int n; // 被分成的份数
 float a,b; // 积分范围
 int i;  // 循环变量

 printf("请输入极限下限:");
 scanf("%f",&a);
 printf("请输入极限上限:");
 scanf("%f",&b);
 printf("请输入被分成的份数:");
 scanf("%d",&n);
 for(i=1; i<=n;i++)
 {
  s1 = s1 + f1(a+((b-a)/n)*i)*((b-a)/n);
  s2 = s2 + f2(a+((b-a)/n)*i)*((b-a)/n);
  s3 = s3 + f3(a+((b-a)/n)*i)*((b-a)/n);
 }
 printf("积分的结果为:%0.4f\n",s1);
 printf("积分的结果为:%0.4f\n",s2);
 printf("积分的结果为:%0.4f\n",s3);
}

/* y=x */
float f1(float x)
{
 float y=0;
 y = x;

 return y;

}

/* y = x*x */
float f2(float x)
{
 float y=0;
 y = x*x;
 return y;
}

/* y = sin(x) */
float f3(float x)
{
 float y=0;
 y = (float)sin(x);
 return y;
}

结果为:

这样就把定积分求完了,当然还有其他的方法。

加载全部内容

相关教程
猜你喜欢
用户评论