c++ drugi wynik wielomianu

#include <iostream>
#include <math.h>
#include <iomanip>

using namespace std;

int sWielomianu, p, q;
float liczbaStopni[100];


float algHor(int k, float x){ //algorytm Hornera - obliczanie wartosci wielomianu
    if (k == sWielomianu){
        return liczbaStopni[sWielomianu];}
    else{
        return algHor(k + 1,x)* x + liczbaStopni[k];}
}


float algSTs(int j){           //Algorytm Show-Trauba funkcja pomocnicza s(j)
    return (sWielomianu - j) % q;
}

float algSTR(int j){           //Algorytm Show-Trauba funkcja pomocnicza r(j)
    if (j % q == 0)
        return q;
    else
        return 0;
}

float T(int i, int j, float x){ //Algorytm Show-Trauba - glowna funkcja
    if (x == 0)                    
        return liczbaStopni[j];
    else
        if (j == -1)
            return liczbaStopni[sWielomianu - i - 1]* pow(x,algSTs(i + 1));
        else
    if (i == j)
        return liczbaStopni[sWielomianu] * pow(x,algSTs(0));
    else
        return T(i-1, j-1, x)+T(i-1, j, x) * pow(x,algSTR(i - j));
}

float pochodna(int stopien, float punkt){
    if (punkt == 0)
        return T(sWielomianu,stopien,punkt);
    else return T(sWielomianu,stopien,punkt) / pow(punkt,stopien%q);
}

main(){

    int k;
    int iteracja = 1000; // Liczba interacji
    float poczatek, koniec, wynik;
    p = 1;
    q = sWielomianu + 1;

 
    cout << "\nPodaj stopien wielomianu: ";
    cin >> sWielomianu;
        if (sWielomianu > 100){
            cout << "Za duzy stopien" << endl;
        }
        if (sWielomianu < 2){
            cout << "Za maly stopien";
        }

    cout << "\nPodaj teraz wspolczynniki wielomianu (zaczynij od tego z najwieksza potega).\n";
    for(k = sWielomianu; k >= 0; k--){
        cout <<"a" << k << " - ";
        cin >> liczbaStopni[k];
    }

    cout << "\nPodaj poczatek przedzialu: ";
    cin >> poczatek;
    cout << "Podaj koniec przedzialu: ";
    cin >> koniec;


    if (pochodna(2,poczatek)*algHor(0,poczatek)>0){
            wynik = poczatek;
    } else {
        wynik = koniec;
    }

    for (k = 1; k <= iteracja; k++){
        wynik = wynik - (algHor(0,wynik)/pochodna(1,wynik));
        if (algHor(0,wynik)==0){
            break;
        }
    }


    if (algHor(0,wynik) == 0){
        cout << "\nDokladny pierwiastek wynosi: " << setprecision(5) <<  wynik << endl;
    } else{

        cout << "\nPrzyblizony pierwiastek wynosi: "  << setprecision(5) << wynik << endl;
    }

    return(0);
}