C语言编程-8

    xiaoxiao2025-04-14  34

    编写函数,用牛顿迭代法求方程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; }

     

    最新回复(0)