Newton- Raphson ver2

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

void oblicz(void)
{
    double p;
    double epsilon;
    int ilosc_iteracji;

    cout<<"Podaj liczbe nieujemna do wyliczenia pierwiasteka kwadratowego: ";
    cin>>p;

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

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

    double a = p;
    int i=0;

    while(abs(a-(p/a))>=epsilon && i<ilosc_iteracji)
       {
           cout<<"Krok "<<i+1<<endl;
            double roznica= abs(a-(p/a));
            a=(a+(p/a))/2;
            double b=p/a;
            cout<<"obecnie a wynosi: "<<fixed<<setprecision(2)<<a<<endl;
            cout<<"obecnie b wynosi: "<<fixed<<setprecision(2)<<b<<endl;

            cout<<"roznica miedzy bokami wynosi: "<<fixed<<setprecision(6)<<roznica<<endl<<endl;

           i++;
       }

    cout<<"\n\nPierwiastek z liczby "<<fixed<<setprecision(2)<<p<<" wynosi ";
    cout<<fixed<<setprecision(4)<<a<<endl<<endl<<endl;
}

int main()
{
    cout<<"PROGRAM LICZACY PIERWIASTEK KWADRATOWY Z LICZBY NIEUJEMNEJ \nALGORYTMEM NEWTONA-RAPHSONA (metoda HERONA)\n"<<endl;

    oblicz();

    system("pause");
    return 0;
}