[已解决]求导问题

jiuhu

求大家帮忙看一下这段代码有什么问题，总是死循环。
求函数导数
#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);
}

最佳答案

月排行榜 / 总排行榜

superbe

2019-10-31 12:05:15

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

跳转到最佳答案楼层

ba21

#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

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

？

superbe

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

jiuhu

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

感谢老哥