math.h 의 sqrt()를 사용하지 않고 제곱근 구하기

<math.h>의 sqrt()를 사용하지 않고 제곱근을 구하는 소스
NUM_REPEAT값은 16~20 이면 충분하다.
큰수의 제곱근을 구할경우 NUM_REPEAT값을 더 큰수로 정하면 된다.

원리는 뉴튼-랩슨법을 이용한 근사값 구하기

double mysqrt(unsigned int src)
{
    unsigned int NUM_REPEAT = 16;
    unsigned int k;
    double t;
    double buf = (double)src;
    for(k=0,t=buf;k<NUM_REPEAT;k++)
    {
        if(t<1.0)
            break;
        t = (t*t+buf)/(2.0*t);
    }
    return t;
}