我正在使用辛普森一家的规则创建一个用于解决积分的C程序,我编写并运行它,但在给程序赋值后,它总是返回一个0.0000的定积分的值。我重新检查了它的每一行似乎都很好,这里的代码可以帮助解决这个问题
#include<stdio.h>
#include<math.h>
float ad(int a, int b, int n)
{
float f4, f2, z, v, q4=0, q2=0, d, m, fa, fb;
int w=2, g=1, k=1, j=1;
z=(b-a)/n;
fa=6*pow(a,2)+16*pow(a,3);
fb=6*pow(b,2)+16*pow(b,3);
f4=6*pow(a+z*w,2)+16*pow(a+z*w,3);
f2=6*pow(a+z*g,2)+16*pow(a+z*g,3);
v=fa+fb;
m=v*z;
while(k<=n/2)
{
q4=q4+(z/3)*(4*f4);
w=w+2;
k++;
}
while(j<=(n-2)/2)
{
q2=q2+(z/3)*(2*f2);
g=g+2;
j++;
}
d=m+q4+q2;
return d;
}
main()
{
int x, y, l;
float o;
printf("Enter number x: ");
scanf("%d", &x);
printf("Enter number y: ");
scanf("%d", &y);
printf("Enter an even number: ");
scanf("%d", &l);
if(l%2!=0)
{
printf("The number is odd!\n");
return 1;
}
o=ad(x, y, l);
printf("The aprox integral is es: %f\n", o);
return 0;
}
它也给了我这两个错误:
--------------------Configuration: mingw5 - CUI Debug, Builder Type: MinGW--------------------
Checking file dependency...
Compiling E:\anti simpson\ad.cpp...
[Warning] E:\anti simpson\ad.cpp:29: warning: converting to `int' from `float'
[Warning] E:\anti simpson\ad.cpp:50:2: warning: no newline at end of file
Linking...
Complete Make ad: 0 error(s), 2 warning(s)
Generated E:\anti simpson\ad.exe
答案 0 :(得分:3)
这一行存在一个问题:
int ad(int a, int b, int n)
将其更改为:
float ad(int a, int b, int n)
这一行和上面类似的一行让我感到困惑:
q2=0;
q2=q2+(z/3)*(2*f2);
为什么将其设置为零然后将其设置为值。我希望它应该在while循环之前设置为零。
答案 1 :(得分:3)
你声明你的函数返回int,但你返回一个浮点数, 结果是: 你要返回的浮动被截断,你只得到int。
我猜你所有的积分都有0到1之间的值,所以函数只返回0
只需将int ad(int a, int b, int n)
更改为float ad(int a, int b, int n)
修改强>:
z=(b-a)/n;
所有a,b和n都是整数,你不会在这个分区中获得小数部分。
尝试z=(b-a)/(n * 1.0);
只是为了让其中一个opreands浮动,所以你也得到了小数部分