Babylonian-SQRT.cpp

Babylonian-SQRT.cpp

#include <iostream>

using namespace std;

//Divide Et Impera
double babylonianSquareRoot(double number) 
{
	if(number < 0) {
        
        cerr<<"Error: Cannot compute Square Root of a negative number."<<endl; 
		return -1;
	}


	double x = number; //initial
	double y = 1.0;

	double e = 0.000001; //tolerance for convergence

	while((x - y) > e) {

		x = (x + y) / 2.0; //update X as the average of x and y;

		y = number / x;
	}

	return x;
}

int main(int argc, char const *argv[])
{
	
	double number;

    cout<<"Enter a number to compute its Square Root: ";
    
    cin>>number;

    double result = babylonianSquareRoot( number );

    if( result != -1 ) {

    	  cout<<"The Square Root of "<<number<<" is approximately: "<<result<<endl; 
    }
 

	return 0;
}