Błędy w programie wyliczającym pole kwadratu z trójkąta

Witam,
mój program ma za zadanie obliczyć maksymalne pole kwadratu, jaki można wyciąć z trójkąta równobocznego o wpisanym polu. Ma zaokrąglać je do dwóch miejsc po przecinku. niestety wyniki wychodzą inne, niż w zadaniu. dla 10 powinno być 4.97, dla 59 29.35, a dla 69 34.32

#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
int main(){
long double P, a, b, P2, h;
cin>>P;
b=sqrt(4+P/sqrt(3));
h=b*sqrt(3)/2;
a=b*h/(b+h);
P2=pow(b,2);
cout<<fixed<<setprecision(2)<<P2;
return 0;
}

masz złe wzory