sub sqrt ($) {
  my $x = shift;
  return sqrt_search( $x, 0, $x, 1, 50);
}

# Tests

sub test
{
    my $n = shift;
    print "SQRT $n: ", &sqrt( $n ), ' ?= ', eval { CORE::sqrt( $n ) } 
+, $/, $/;
}

test( $_ ) for( 16, 9, 4, 2, 0.5, -1, 0, 0.001 );
[download]

SQRT 16: 4.0000114440918 ?= 4

SQRT 9: 3 ?= 3

SQRT 4: 1.99998474121094 ?= 2

SQRT 2: 1.4141845703125 ?= 1.4142135623731

SQRT 0.5: 0.7071533203125 ?= 0.707106781186548

SQRT -1:  ?= 

SQRT 0: 0.0078125 ?= 0

SQRT 0.001: 0.03125 ?= 0.0316227766016838
[download]