Metoda Seidla implementacja

0

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:
user image

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

0

Indeksowanie tablic w C/C++ zaczyna się od 0

0

przy wczytywaniu n zwiększam je o 1 pomijam 0 indeks. Wiem ze to nie jest optymalnie ale Pani profesor sobie tak życzyła.

1 użytkowników online, w tym zalogowanych: 0, gości: 1