c语言用牛顿法编程计算,C语言编程-8_8 公式求解 牛顿迭代法求解方程

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

#include

#include

// 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

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值