Obliczanie pierwiastka kwadratowego - metodą Newtona-Raphsona

Specyfikacja problemu algorytmicznego.
Problem: Obliczanie pierwiastka kwadratowego z liczb całkowitych
Dane wejściowe:
pole c C dla pole>0 - liczba podpierwiastkowa
epsilon c R dla epsilon>0 - zadana dokładność
Dane wyjściowe: bok c R dla bok>0 - wartość pierwiastka

Lista kroków
1. Pobierz pole,epsilon
2. Za bok podstaw 1
3. Jeżeli fabs(pole-bok*-bok)<=epsilon to wypisz bok i zakończ
4. W przeciwnym razie wylicz:
bok=(bok-pole/bok)/2
bok=fabs(bok)
i wróć do kroku trzeciego

do góry