Listing 4

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

void main()
{
 double x,y ;
 
 for ( x = 1.0; x <= 25.0; (x = x + 1.0) )
  {
   y = sqrt(x) * sqrt(x);
   if (y != x)
     printf("%5.1f %21.17f %10.2e\n", x, y, x-y);
  }
}