Tính căn bậc 2

Yêu cầu bài toán: tính căn bậc 2 theo CT:
x[0]=(a+1)/2
x[n+1]=(1/2) * (x[n] + a/x[n]) voi n=1,2,3... với độ chính xac 0.0001

Mã:

#include<stdio.h>
#include<conio.h>
#include<math.h>
void main()
{
	clrscr();
	double x,a,tg;
	printf("Nhap so can tinh can bac 2: ");
	scanf("%lf",&a);
	tg=x=(a+1)/2;
	do
	{
	     tg=x;
	     x=0.5*(tg+a/tg);
	}
	while(fabs(tg-x)>0.0001);
	printf("\nKet qua :  %lf",x);
	printf("\nHam sqrt:  %lf",sqrt(a));
	getch();
}