迭代法求平方根

#include<iostream>
#include<math.h>
using namespace std;
int main()
{
    double a,x0,x1;
    cin>>a;
    if(a<0)
        cout<<"错误"<<endl;
    else
    {
        x0=a/2;
        x1=(x0+a/x0)/2;
        do
        {
            x0=x1;
            x1=(x0+a/x0)/2;
        }
        while(fabs(x0-x1)>=1e-6);
    }
    cout<<x0<<endl;
    return 0;
}

求平方根的迭代公式：x1=1/2*(x0+a/x0）。
算法：1.先自定一个初值x0，作为a的平方根值，在我们的程序中取a/2作为a的初值；利用迭代公式求出一个x1。此值与真正的a的平方根值相比，误差很大。
⒉把新求得的x1代入x0中，准备用此新的x0再去求出一个新的x1.
⒊利用迭代公式再求出一个新的x1的值，也就是用新的x0又求出一个新的平方根值x1，此值将更趋近于真正的平方根值。
⒋比较前后两次求得的平方根值x0和x1，如果它们的差值小于我们指定的值，即达到我们要求的精度，则认为x1就是a的平方根值，去执行步骤5；否则执行步骤2，即循环进行迭代。
迭代法是用于求方程或方程组近似根的一种常用的算法设计方法。设方程为f(x)=0，用某种数学方法导出等价的形式x=g(x），然后按以下步骤执行：
⑴ 选一个方程的近似根，赋给变量x0；
⑵ 将x0的值保存于变量x1，然后计算g(x1），并将结果存于变量x0；
⑶ 当x0与x1的差的绝对值还小于指定的精度要求时，重复步骤⑵的计算。
若方程有根，并且用上述方法计算出来的近似根序列收敛，则按上述方法求得的x0就认