#include <stdio.h>

void main()
{float x, y, z1, z2;
 x = 1.0;
 y = 1.0e-20;
 rd_near();
 z1 = x - y;
 z2 = y - x;
 z1 = z1 - x;
 z2 = z2 + x;
 printf("near, z1 = %17.10e, z2 = %17.10e \n",z1, z2);
 rd_minf();
 z1 = x - y;
 z2 = y - x;
 z1 = z1 - x;
 z2 = z2 + x;
 printf("minf, z1 = %17.10e, z2 = %17.10e \n",z1, z2);
 rd_pinf();
 z1 = x - y;
 z2 = y - x;
 z1 = z1 - x;
 z2 = z2 + x;
 printf("pinf, z1 = %17.10e, z2 = %17.10e \n",z1, z2);
 rd_zero();
 z1 = x - y;
 z2 = y - x;
 z1 = z1 - x;
 z2 = z2 + x;
 printf("zero, z1 = %17.10e, z2 = %17.10e \n",z1, z2);
};