CALCUL_RAD_PATR

#include <cstdio>
#include <cmath>
#define eps 0.00001
using namespace std;
float radical(float a)
{
/// avand valoarea egala cu precizia dorita in calculul radicalului
 float xn1=1, xn=(1+xn1)/2;
 do { xn1=xn;
 xn = (xn1+a/xn1)/2;
 } while (fabs(xn-xn1)>eps);
 /// fabs(x) = |x|, unde x este numar real
 return xn;
/// rezultatul este ultimul termen calculat al sirului
}

int main()
{
 float a;
 printf("Dati a:"); scanf("%f",&a);
 printf("\nsqrt (%7.3f)=%10.6f",a,sqrt(a));
 /// se compara cele doua functii
 printf("\nradical(%7.3f)=%10.6f",a,radical(a));
 printf("\nCinci zecimale exacte!");
}