Question

大一菜鸟求教C++大神，矩形法求定积分，为何程序结果总是0.

#include<iostream>
#include<cmath>
using namespace std;
float fun(int a,int b,float(*f)(float x))
{
float d,x,s=0;
d=(b-a)/100;
for(int i=0;i<100;i++)
{
s=s+f(x)*d;
x+=d;
}
return s;
}
float fcos(float x)
{
return cos(x);
}
void main()
{
cout<<fun(0,2,fcos)<<endl;
}

Answer 1

修改无误，注意//的注释

#include<iostream>

#include<cmath>

using namespace std;

float fun(int a,int b,float(*f)(float x))

{

 float d,x=a,s=0;//x没有赋初值，他的初值应该是积分下限a

 //最关键的问题：a,b都是int类型的 (b-a)/100是整除，只要b-a的值小于100，他的值一直是0

 //当然是错的，必须转成浮点数做除法

 d=(b-a)*1.0/10000;//这儿最好设大一点，这样精度会高些 100太小了，乘个1.0转成浮点数。

 for(int i=0;i<10000;i++)

 {

      s=s+f(x)*d;

      x+=d;

 }

 return s;

}

float fcos(float x)

{

  return cos(x);

}

void main()

{

  cout<<fun(0,2,fcos)<<endl;

}

结果：

Answer 2

#include<iostream>
#include<cmath>
using namespace std;
float fun(int a,int b,float(*f)(float x))
{
float d,x,s=0;
d=(b-a)/100.0; //这里不能是整数除整数 ，把100改为100.0
for(int i=0;i<100;i++)
{
s=s+f(x)*d;
x+=d;
}
return s;
}
float fcos(float x)
{
x=x*3.1415926/180; //这里转成弧度 ，就ok
return cos(x);
}
void main()
{
cout<<fun(0,2,fcos)<<endl;
}