求一个数的算术平方根

写一段程序求出一个数的算术平方根：

#include<iostream>
using namespace std;
#define e 0.001
double sqrt(double a)
{
		double x,y;
		x=a/2;
		y=x+1+e;
		while(x-y>e||y-x>e)
		{
				y=x;
				x=(x+a/x)/2;
		}
		return x;
}
int main()
{
		double a;
		cin>>a;
		cout<<sqrt(a);
		return 0;
}