/* Quadratwurzel ohne Fehlerausgabe fuer negative Radikanden       */
#include <stdio.h>

#define  DELTA  1.0e-16

double wurzel(double);   /* Funktionsprototyp */

int main() {
  int i;
  for (i=0; i < 10; i++)
    printf( "%d  %20.16f\n", i, wurzel(i) );
}

double betrag(double x) {
  return x < 0 ? -x : x;
}

double wurzel(double n) {
  double  x = n, alt_x = 0;
  if (n <= 0)
    return 0;

  while ( betrag(x - alt_x) > DELTA ) {
    alt_x = x;
    x = (n/x + x)/2;
  }
  return x;
}