Dovresti calcolarti l'epsilon macchina (non ricordo se serve quello assoluto o relativo) e confrontarlo con un valore significativo dell'uguaglianza.
Qualcosa del tipo (siano 'a' e 'b' i due reali da confrontare):

if ( abs((a/b)-1) < eps_macc)
Pero' me lo sono inventato al momento, potrebbe non funzionare!!!

P.S. L'epsilon macchina e' il piu' piccolo numero floating point x tale che:

(x + 1) > 1

Dovrebbe essere pari a 2^(-m), dove m sono i bit dedicati alla mantissa.