牛顿迭代法

最新推荐文章于 2024-03-21 11:08:33 发布

霜刃未曾试

最新推荐文章于 2024-03-21 11:08:33 发布

阅读量597

点赞数

分类专栏：牛顿迭代法

本文链接：https://blog.csdn.net/discreeter/article/details/77163252

版权

牛顿迭代法专栏收录该内容

3 篇文章 0 订阅

订阅专栏

今天粗略看了一下牛顿迭代法。。。
这篇生动形象简单易懂：https://www.zhihu.com/question/20690553
这篇用牛顿迭代法求平方根：http://blog.csdn.net/wumuzi520/article/details/7026808
牛顿迭代法求平方根：

#include <bits/stdc++.h>

using namespace std;

const double eps = 1e-5;
double y;

double df(double x) //函数的导数
{
    return 2*x;
}
double f(double x) //要求根的函数
{
    return x*x - y;
}
double newton_iteration(double x)
{
    int tot = 0;
    //这里不用eps控制循环的原因是当输入的数很大时，经过一定次数的迭代求出的根已经足够精确了，但是用f(x)放大x后造成了巨大的误差，导致不能结束循环。。。
//    while(fabs(f(x) - 0) > eps)
//    {
//        x = x - f(x) / df(x);
//        if(++tot > 100) return -1;
//    }
    while(++tot <= 100)
    {
        x = x - f(x) / df(x);
    }
    return x;
}
int main()
{//不能处理负数的喔
    while(~ scanf("%lf", &y))
    {
        //printf("%.6f %.6f\n", sqrt(y), pow(y, 0.5));
        double ans = newton_iteration(y / 2.0);
        if(fabs(-1 - ans) < eps) puts("无解");
        else printf("%.6f\n", ans);
    }
    return 0;
}