牛顿迭代法c代码

㈠ 牛顿迭代法的C语言代码

doublefunc(doublex)//函数{returnx*x*x*x-3*x*x*x+1.5*x*x-4.0;}doublefunc1(doublex)//导函数{return4*x*x*x-9*x*x+3*x;}intNewton(double*x,doubleprecision,intmaxcyc)//迭代次数{doublex1,x0;intk;x0=*x;for(k=0;k<maxcyc;k++){if(func1(x0)==0.0)//若通过初值，函数返回值为{printf(迭代过程中导数为0! );return0;}x1=x0-func(x0)/func1(x0);//进行牛顿迭代计算if(fabs(x1-x0)<precision||fabs(func(x1））<precision)//达到结束条件{*x=x1;//返回结果return1;}else//未达到结束条件x0=x1;//准备下一次迭代}printf(迭代次数超过预期！ );//迭代次数达到，仍没有达到精度return0;}intmain(){doublex,precision;intmaxcyc;printf(输入初始迭代值x0:);scanf(%lf,&x);printf(输入最大迭代次数：);scanf(%d,&maxcyc);printf(迭代要求的精度：);scanf(%lf,&precision);if(Newton(&x,precision,maxcyc)==1）//若函数返回值为1printf(该值附近的根为：%lf ,x);else//若函数返回值为0printf(迭代失败！ );getch();return0;}

㈡ C语言程序设计 迭代法

给你一点提示。
牛顿迭代法要计算
(1)
y1=f(x)
在
x
的函数值
(2)
d1=f(x)
的一阶导数
在
x
的值
你可以写两个函数，分别计算y1,d1
如果一阶导数有解析解，则可用赋值语句，否则要写数值解子程序。
步骤：
设解的精度,例
float
eps=0.000001;
设x初值，x1;
算y1=f(x1);
迭代循环开始
算一阶导数
在
x1
的值
d1
用牛顿公式
算出
x2;
[x2
=
x1
-
y1
/
d1]
如果
fabs(x2-x1)
>
eps
则从新迭代
--
用新的函数值和一阶导数值推下一个
新x.