jiuhu 发表于 2019-10-29 16:47:33

求导问题

求大家帮忙看一下这段代码有什么问题,总是死循环。
求函数导数
#include <stdio.h>
#include <math.h>

float f(float x)
{
        float s;
        s=x*x*x*x*x*x-x-1;
        return s;
}
int main(void)
{
        float dx=1,dy,dd1,dd2;
int i=1;
float x0=1.0;

dy = f(x0) - f(x0-dx);

dd1=dy/dx;
do
{
        dx = 0.5 * dx;
dy = f(x0) - f(x0-dx);
dd2=dy/dx;


if ( fabs(dd1-dd2) < 1e-06 )
{
        printf("jieguo%f",dd2);
        i=0;
}
else dd1=dd2;
}while(i=1);
}

ba21 发表于 2019-10-29 18:08:09

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

float f(float x)
{
      float s;
      s=x*x*x*x*x*x-x-1;
      return s;
}
int main(void)
{
      float dx=1,dy,dd1,dd2;
int i=1;
float x0=1.0;

dy = f(x0) - f(x0-dx); // 0.00000

dd1=dy/dx; // 0.00000

}

多了不想说,自己用点耐心先一句一句的调试。这几处的0这代码是什么玩意?



jiuhu 发表于 2019-10-29 19:21:39

ba21 发表于 2019-10-29 18:08
#include
#include




superbe 发表于 2019-10-31 12:05:15

本帖最后由 superbe 于 2019-10-31 12:06 编辑

最后那个 while(i=1) 应该是 while(i==1),改了测试结果正确。
也可以用break跳出循环。

jiuhu 发表于 2019-11-5 19:47:03

superbe 发表于 2019-10-31 12:05
最后那个 while(i=1) 应该是 while(i==1),改了测试结果正确。
也可以用break跳出循环。

感谢老哥
页: [1]
查看完整版本: 求导问题