Hej!
Mam pewien problem. Przyznam od razu, nie jestem zbyt dobra z c++ czy z matmy. Mam znaleźć miejsca zerowe wielomianu (używając metody Newtona). Znalazłam kod, trochę go przerobiłam...Program dobrze oblicza, ale tylko jedno miejsce zerowe. Jak są dwa, nadal wyświetla tylko jedno. Czy mógłby mi ktoś pomóc, przerobić go, aby otrzymać pozostałe?
Kod:
#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);
}
C++ 
