标签:f1 xn 函数 求解 double eps C语言 return 迭代法
编写函数,用牛顿迭代法求方程f(x)=2x3-4x2+3x-6=0在1.5附近的根。
牛顿迭代公式为:xn+1=xn-f(xn)/f'(xn)
其中,f'(xn) 是f在xn处的导数。
结束条件:|f(xn+1)|< eps与|xn+1-xn|< eps同时成立(eps是一个很小的正数,从键盘输入)
同时编写主函数,在主函数中调用并输出函数值。
函数原型如下:
double root(double (*f)(double),double (*f1)(double),double x,double eps);
其中,形参f接收原函数f(x),形参f1接收导函数f'(x)
如输入:
1e-8
如输出结果为:
2.000000
注意 浮点数的绝对值函数是 fabs,整数的是abs,不能混用。
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
// f(x) = 2x^3 -4x^2 +3x - 6 = 0 在1.5处的根
// f'(x) = 6x^2 - 8x + 3
double fun(double x)
{
x = 2 * x*x*x - 4 * x*x + 3 * x - 6;
return x;
}
double fun1(double x)
{
x = 6 * x*x - 8 * x + 3;
return x;
}
double root(double(*f)(double), double(*f1)(double),double x, double eps)
{
double x0;
while (1)
{
x0 = x;
x = x0 - (*f)(x) / (*f1)(x);
if (fabs((*f)(x)) < eps && fabs(x - x0) < eps)
{
break;
}
}
return x;
}
int main()
{
double x, eps;
x = 1.5;
scanf_s("%lf", &eps);
x = root(fun,fun1,x, eps);
printf("%lf", x);
return 0;
}
标签:f1,xn,函数,求解,double,eps,C语言,return,迭代法 来源: https://blog.csdn.net/qq_36941368/article/details/90578624
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。