Witam

Znalazłem na angielskiej wikipedii numeryczny wzór dla metody Jacobbiego więc go zaimplementowałem ( oczywiście brak sprawdzania dzielenia przez zero, ale pomińmy ten fakt na razie ) Oto kod:

double B[17]; double X[17];double w;
double XX[17];

cout<<"Wprowadz dane"<<endl;

for(i=0;i<17;i++)
{
    cin>>w;
    B[i]=w;
}

for(i=0;i<17;i++) {X[i]=0;XX[i]=0;}
for(int k=0;k<5;k++)
{
for(i=0;i<17;i++)
{
    double s=0;
    for(j=0;j<17;j++)
    {
        if(j!=i) s+=tab[i][j]*XX[i];
    }
    X[i]=(B[i]-s)/(tab[i][i]);
    XX[i]=X[i];
}
}

Dane wejściowe mam w pliku, Matlab zwraca mi dobre wyniki natomiast program dla małych układów ok, dla mojego jakieś bzdury..

Proszę o pomoc, pozdrawiam