迭代法求平方根C语言

#include<stdio.h>

#include<math.h>

int main()

{

long int n;

int i=0;

scanf("%ld",&n);

double a[1000];

a[0]=(1+n/1)/2;

a[1]=(a[0]+n/a[0])/2;

while(fabs(a[i+1]-a[i])>0.00001)

{

a[i+1]=(a[i]+n/a[i])/2;

i++;

}

printf("%.3f",a[i]);

return 0;

}