Problem z double (niedokładny wynik)

#include<iostream>
#include <iomanip>
#include <cmath>
long double poPrzecinku( long double liczba, long poPrzecinku )
{
    long long iLiczba = 10;
    iLiczba = static_cast < long long >( pow( iLiczba, poPrzecinku ) );
   
    return round( liczba * iLiczba ) / iLiczba;
}

using namespace std;
int main()
{
    double a = 0, b = 5, Pole = 0;
    cout << "b= " << b << endl;
    cout << "a= " << a << endl;
    cout << "Pole= " << Pole << endl;
    cout << "Operacja pierwsza.....\n";
    a = b / sqrt( 2 );
    a = poPrzecinku( a, 5 );
    cout << "b= " << b << endl;
    cout << "sqrt(2)= " << sqrt( 2 ) << endl;
    cout << "a= " << a << endl;
    cout << "Druga operacja...\n";
    a = pow( a, 2 );
    cout << "a= " << a << endl;
    cout << "Pole= " << Pole << endl;
    cout << "Pole tego tr˘jkĄta wynosi: " << setprecision( 8 ) << a;
    system( "pause>>NULL" );
    return( 0 );
}