Siento mucho decirte que si te quieres complicar vas a tener que trabajar a nivel de bit con el estándar IE754, no es complicado separas parte entera mantisa y exponente y sacas los datos.

NOB2014, puedes aumentar la precisión cambiando en donde dice FLT_EPSILON por la constante DBL_EPSILON

En este caso cambias del épsilon de la máquina para punto flotante al épsilon para precisión doble. Creo que al manejar más de nueve decimales estás considerando una precisión doble, por eso requieres DBL_EPSILON

Saludos, Yoel

Por fin he solucionado lo que querías amigo y sí hay una forma estándar y sencilla de hacerlo te paso mi código y espero que te guste sigue las instrucciones, te pido por favor que me envies un mensaje si te funciono es para cuestiones estadísticas que llevo

#include <stdlib.h>
#include <stdio.h>
 
int main()
{
   div_t d;
   int num, denom;
   double aux;
   printf( "Escriba el numerador y el denominador (separados por un espacio):\n" );
   scanf( "%d %d", &num, &denom );
   d = div( num, denom ); 
    aux=d.rem/(double)denom;
   printf( "ldiv( %d, %d ) : cociente = %d, resto = %f\n", num, denom, d.quot, aux );
 
   return 0;
}

Espero que me contestes

---

Esta solución permite 30 decimales !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

 
 
#include <stdlib.h>
#include <stdio.h>
 
int main()
{
   div_t d;
   int num, denom;
   double aux;
   printf( "Escriba el numerador y el denominador (separados por un espacio):\n" );
   scanf( "%d %d", &num, &denom );
   d = div( num, denom ); 
    aux=d.rem/(double)denom;
   printf( "ldiv( %d, %d ) : cociente = %d, resto = %.30lf\n", num, denom, d.quot, aux );
 
   return 0;
}
 
 
 
 
 
 

Así es jeje no me fije, pero ya con una variable más grande tiene más precisión, un long double, me parece, de cualquier forma creo que mi solución es la que parece más eficaz no sé qué le parezca al dueño del post.