#include <stdio.h>
#include <math.h>

int
main ()
{
  double x1, x2, p, q;
  double x1esatto;

  x1esatto = 1/3.0;
  printf ("x1 esatto = %20.18lg\n", x1esatto);
  x2 = 1e10;

  p = (x1esatto + x2)/2.0;
  q = x1esatto*x2;

  printf ("Algoritmo 1:\n");
  x1 = p - sqrt(p*p -q);
  printf ("x1 = %20.18lg\n", x1);
  printf ("errore: %lg\n", x1esatto - x1);

  printf ("Algoritmo 2:\n");
  x1 = q/(p + sqrt(p*p -q));
  printf ("x1 = %20.18lg\n", x1);
  printf ("errore: %lg\n", x1esatto - x1);
}