#include <stdio.h>
#include <math.h>
double fh(double h)
{
double kao,kat,kath,y,c,kw,hta,ha,a,oh;
kw=pow(10,-14)*pow(10,7);
kao=0.059;
kat=0.000064;
kath=0;
c=0.1;
kao=kao*pow(10,7);
kat=kat*pow(10,7);
kath=kath*pow(10,7);
c=c*pow(10,7);
y=h*h+kao*h+kao*kat;
hta=c*h*h/y;
ha=c*kao*h/y;
a=c*kao*kat/y;
oh=kw/h;
double fh;
fh=h-ha-2*a-oh;
return fh;
}
main()
{
double x,e,b,r,dlt;
e=0;
b=10000000;
dlt=1;
while(fabs(b-e)>dlt)
{
r=(e+b)/2;
if(fh(r)==0)
{break;}
else
{
if(fh(r)*fh(b)<0)
{
e=r;
}
else
{
b=r;
}
}
}
printf ("%f",r);
}