用牛顿迭代法求某正数n的平方根

#include<iostream>
#include<math.h>
using namespace std;
int main()
{
	double a,x0,x1;
	cout<<"输入数字A:";
	cin>>a;
	x0=a/2;
	x1=(x0+a/x0)/2;
	do
	{
		x0=x1;
		x1=(x0+a/x0)/2;
	}while(fabs(x0-x1)>1e-10);
	cout<<x1<<endl;
	return 0;
}