void bisekcja(double a, double b, double epsilon) {

    cout <<endl<<endl;
    cout << "Metoda bisekcji"<<endl;
    cout << "a = " <<a<<endl;
    cout << "b = " <<b<<endl;

    double s;
    int i=0;
    do {
        //cout <<fabs(funkcja(s))<<endl;
    i++;
    s = (a+b) / 2;

    if(abs(funkcja(s)) < epsilon)
			break;

    if( funkcja(a) * funkcja(s) < 0 ) {
        b = s;
    } else {
        a = s;
    }

    //cout << fabs( f(s) )<<endl;

    //cout <<fabs(funkcja(s))<<endl;
    } while( (b-a ) >= epsilon);


    cout << "X = "          << s<<endl;
    cout << "Dokladnosc = " << epsilon<<endl;
    cout << "kroki = " << i<<endl;
}