Нахождение числа методом половинного деления

#include <iostream>
using namespace std;

int main()
{


float a=0;
float b=1;
float x=3;
while ((b-a)>=0.001)
{
float c=(a+b)/2.0;

if (c*c<=x)

a=c;
else
b=c;
cout<<c<<endl;
}

}