Neliöjuuren laskenta

Ompas se yksinkertainen juttu. Vaikka arvaus olisikin aika kaukana oikeasta, lopputulos löytyy ihmeen nopeasti. Periaate on tuhansia vuosia vanha mutta soveltuu tietokoneaikaan ihmeen hyvin.

#include <stdio.h> //Neliöjuuren laskeminen iteraatiolla
double luku=0.2,juuri=1, edel=0; //luku=juurrettava luku, juuri=arvaus
int lkm=0;
int main ()
{
while (edel!=juuri){
edel=juuri;
juuri=(juuri luku/juuri)/2;
lkm ;
printf("%i b=%.26f\n",lkm,juuri);
}
printf("Luku=%g juuri=%g\n",luku,juuri);
return 0;
}
Käännetään c-kielessä näin:
gcc -o k nj.c&&./k