集成功能时意外的结果

作者: 雪浴冰灵
发布时间: 2024-02-21 03:07:42 (3月前)
转自：

3 条回复

0#
回复此人
记忆短浅♡思念不变 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> 改革 <code> for() </code> 循环迭代和精确 整数 次数。 OP的;代码正在改变 <code> a </code> 每个循环如此 <code> i < ((b - a) * 1000 </code> 不会导致正确的迭代计数。 <a href="https://stackoverflow.com/questions/55081244/unexpected-result-in-integration-of-a-function/55081408#comment96908438_55081244"> @Michail </A> <pre> <code> #include <stdio.h> #include <math.h> int main() { double a, b, c; int i; c = 0; printf("Insert limit 1:\n"); a = 0; //scanf("%lf",&a); printf("Insert limit 2:\n"); b = 1.0; /// scanf("%lf",&b); #if 0 for (i = 0; i < ((b - a) * 1000); i = i + 1) { c = c + 0.001 * pow(a, a); a = a + 0.001; } #else double delta = (b - a) / 1000; for (i = 0; i < 1000; i = i + 1) { c = c + delta * pow(a, a); a = a + delta; } #endif printf("The area is %lf\n", c); return 0; } </code> </pre> 产量 <pre> <code> The area is 0.783431 </code> </pre> </DIV>

编辑
1#
回复此人
产品你是狗 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> 循环体修改 <code> a </code> 然后在那里重新评估 <code> for </code> 循环应该是循环不变的。更改循环体： <pre> <code> c=c+0.001*pow(a,a); a=a+0.001; </code> </pre> 至 <pre> <code> double x = i * 0.001 ; c += 0.001 * pow( x, x ) ; </code> </pre> 实际上，您可以使用x作为循环控制变量来进一步简化此操作： <pre> <code> const double dx = (b - a) / 1000 ; for( double x = a; x < b; x += dx ) { c += 0.001 * pow( x, x ) ; } </code> </pre> 此外，每次迭代乘以0.001是不必要的;你可以在循环之后将其应用于相同的效果，在某些情况下可以减少累积误差。 <pre> <code> const double dx = (b - a) / 1000 ; for( double x = a; x < b; x += dx ) { c += pow( x, x ) ; } c *= dx ; </code> </pre> 所以完整的代码变成： <pre> <code> #include <stdio.h> #include <math.h> int main() { double a = 0 ; double b = 1.0 ; printf( "Enter limit 1:"); scanf( "%lf", &a ) ; printf( "Enter limit 2:"); scanf( "%lf", &b ) ; const double dx = (b - a) / 1000 ; double c = 0 ; for( double x = a; x < b; x += dx ) { c += pow( x, x ) ; } c *= dx ; printf( "The area is %lf", c ) ; return 0; } </code> </pre> 例： <pre> <code> Enter limit 1:0 Enter limit 2:1 The area is 0.783431 </code> </pre> </DIV>

编辑

登录后才能参与评论