求导问题
求大家帮忙看一下这段代码有什么问题,总是死循环。求函数导数
#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);
}
#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这代码是什么玩意?
ba21 发表于 2019-10-29 18:08
#include
#include
?
本帖最后由 superbe 于 2019-10-31 12:06 编辑
最后那个 while(i=1) 应该是 while(i==1),改了测试结果正确。
也可以用break跳出循环。 superbe 发表于 2019-10-31 12:05
最后那个 while(i=1) 应该是 while(i==1),改了测试结果正确。
也可以用break跳出循环。
感谢老哥
页:
[1]