Witam,
Mam pewien problem z implementacją wzoru obliczania układu równań za pomocą metody Seidla.
for(int i=1; i<n ; i++){
if(i==1){
for(int j=2 ; j<n ; j++){
x_nowy[1] = (alfa[1][j]*x_stary[j])+beta[1];
}
}else{
for(int j=1, k=i+1 ; j < i-1, k<n ; j++, k++){
x_nowy[i] = (alfa[i][j]*x_nowy[j]) + (alfa[i][k]*x_stary[k]) + beta[i];
}
}
}
Wzór:
Czy ktoś mógłby podpowiedzieć czy podany wzór jest dobrze zaimplementowany? Ew. poprawić czy naprowadzić. Wyniki wychodzą mi w granicy od 0.9 do 1.5 a powinny być bardzo bliskie 1