расчет погрешности

//
#include "stdafx.h"
#include <stdlib.h>
#include <conio.h>
#include <stdio.h>
#include <iostream>
#include <math.h>
using namespace std;
 
static union
{
    float drob;
    unsigned int cel;
} st_un1;
 
void perevod(float chislo)
{
    unsigned int b, e, a, one;
    int poryadok;   
    st_un1.drob = chislo;
    b = st_un1.cel;
    poryadok = b;
    a = b;
    // Находим порядок
    a <<= 1;
    a >>= 24;
    printf("\n \n");
    poryadok = a;
    poryadok = poryadok - 127; // Перевод порядка в дестичный вид
    cout << "Порядок: " << poryadok << endl;
    // Находим мантиссу
    b <<= 9;
    b >>= 9;
    printf("\n");
    b = b | 040000000; // Добавление 1 перед мантиссой
    // Выделение целой части
    a = b;
    a >>= 23 - poryadok;
    cout << "Целая часть: " << a << endl;
    // Выделение дробной части мантиссы
    one = b;
    one <<= 9 + poryadok;
    one >>= 9 + poryadok;
    printf("\n");
    st_un1.cel = one;
    //cout << st_un1.drob << endl;
    printf("\nОбратный перевод числа с погрешностью:\n ");
    cout << a << ".";
    //цикл, переводящий в десятичный код
    for (int i = 0; i < 24; i++)
    {
        one = one * 10;
        e = one;
        //выделение дробной части
        e <<= 9 + poryadok;
        e >>= 9 + poryadok;
        //выделение целой части
        one >>= 23 - poryadok;
        cout << one;
        one = e;
    }
}
 
 
 
int main()
{
    
    float chislo;
    setlocale(LC_ALL, "russian");
    cout << "Введите число: ";
    cin >> chislo;
    
    perevod(chislo);
    cout << endl;
    getch();
    cout << endl;
}