牛頓迭代法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.