С++ Задача, Нужны помощь

 double x, y0, y1, eps;
cin >> y1 >> x >> eps;
do {
  y0 = y1;
  y1 = (y1 + x / y1) / 2;
} while (fabs(y0 * y0 - y1 * y1) >= eps);
cout << y1 << ' ' << sqrt(x);