Попытка найти sqrt чисел между 0 и 1 без функции sqrt

Question

Проблема, которую я не могу решить, - это то, как я нахожу sqrt чисел между 0 и 1. После этого я думаю, что все работает. Единственными номерами ввода, которые приводят к сбою этой программы, являются номера между 0 (не включены) и 1 (не включены)

Любая помощь?

#include <iostream> 
#include <cmath> 
#include <cassert> 
#include <stdlib.h> 
#include <cmath> 

using namespace std; 

double squareroot(double x) /* computes the square root of x */ 

{ 
assert(x >= 0); /* make sure x is not negative*/ 
if (x==0) return 0; 

/* the sqrt must be between xhi and xlo */ 
double xhi = x; 
double xlo = 0; 
double guess = x/2; 

/* We stop when guess*guess-x is very small */ 

while (abs(guess*guess-x) > 0.00001) 
    { 
    if (guess*guess > x) xhi = guess; 
    else xlo = guess; 
    guess = (xhi + xlo)/2; 
    } 

return guess; 
} 

/* Test Stub */ 


int main() 
{ 
    double testvalue; 
    cout << "\n Enter a TestValue= " ; 
    cin >> testvalue; 
    cout << endl; 
    double testresult = squareroot(testvalue); 
    cout << "\n Square Root= " << testresult << "\n" ; 
    } 

Answer 1

Ваша программа рушится из-за этой линии:

while (abs(guess*guess-x) > 0.00001)

Если у вас есть число между 0 и 1, ваш алгоритм всегда приведет к ряду больше, чем 0,00001, что приводит к бесконечному циклу , Исправить будет либо изменение вашего алгоритма на другой метод аппроксимации квадратного корня (см. Метод Ньютона, теорема о неподвижной точке, методы полиномиальной аппроксимации и т. Д.).

Подсказка: попробуйте трассировку кода, чтобы увидеть, что происходит, даже если вам нужно записать его на бумаге