Subversion Repositories tcp_socket

/* sqrt.c */

/* Jonathan P Dawson */

/* 2013-12-23 */

/* find absolute value of a floating point number*/

float fabs(float n){

    if (n < 0.0) {

        return - n;

    } else {

        return n;

    }

}

/* approximate sqrt using newton's method*/

float sqrt(float n){

    float square, x, old;

    x = 10.0;

    old = 0.0;

    while(fabs(old - x) > 0.000001){

        old = x;

        x -= (x*x-n)/(2*x);

    }

    return x;

}

/* test sqrt function*/

void main(){

    float x;

    for(x=0.0; x <= 10.0; x+= 0.1){

        file_write(x, "x");

        file_write(sqrt(x), "sqrt_x");

    }

}