Newton- Raphson

#include<iostream>
#include<cstdlib>
#include<iomanip>
#include<cmath>
#include<conio.h>
using namespace std;

double oblicz(double p, double epsilon, long int ilosc_iteracji)
{
     double a = p;
     int i=0;

       while(abs(a-p/a)>=epsilon && i<ilosc_iteracji)
           {
               a=(a+p/a)/2;
               i++;
           }
       return a;
}

int main()
{   char warunek;
    do
    {
        double p;
        cout<<"PROGRAM LICZACY PIERWIASTEK KWADRATOWY Z LICZBY NIEUJEMNEJ \nALGORYTMEM NEWTONA-RAPHSONA (metoda HERONA)\n"<<endl;
        cout<<"Podaj liczbe nieujemna do wyliczenia pierwiasteka kwadratowego: ";
        cin>>p;

        double epsilon;
        cout<<"Podaj dokladnosc obliczen(epsilon): ";
        cin>>epsilon;

        long int ilosc_iteracji;
        cout<<"Podaj maksymalna ilosc teracji: ";
        cin>>ilosc_iteracji;

        cout<<endl<<endl;
        cout<<"Pierwiastek z "<<p<<" wynosi ";
        cout<<fixed<<setprecision(6)<<oblicz(p, epsilon, ilosc_iteracji)<<endl<<endl;

        cout<<"Czy chcesz wyliczyc kolejny pierwiastek? t/n"<<endl;
        warunek=getch();
        system("cls");
    }
    while(warunek=='t');

    cout<<"KONIEC PROGRAMU"<<endl;
    system("pause");
    return 0;
}