Pomoc w implementacji LU

0

http://metodynumeryczne.pl/rozklad_lu.htm

Jest tu stronka i fajnie, ale mam już problem na początku implementacji pierwszych pętli, może nie byłoby to dziwne gdyby nie to, że nie wydaje mi się zbyt praktyczne bo chyba według powyższej strony komputer sam będzie wiedział ile wynosi x jak mu zapisze

2x=4;.....

A to tak prosto na pewno nie będzie...

Moje pierwsze kroki:

for(int i=0;i<W;i++) //Zwiększenie wiersza
{
for(int j=0;j<K;j++) //Zwiększenie kolumny
{
for(int k=0;k<K-2;k++) //Pętla licząca sumę
{
if(i==k) pomoc=1;
else pomoc=0;
suma=suma+pomoc*tabu[k][j];
}
tabu[i][j]=tab[i][j]-suma;
}
}
0

Za wiele to i tak mi nie daje bo na kartce rozwiąże to bez żadnego problemu ;-) Chodzi mi o to, że na moje te wzory nie mają zastosowania wprost przy implementacji bo komputer sam nie rozwiąże równania:

x=1-2x...

Trzeba to poprzekształcać możesz mi pomóc? Najgorzej to powyznaczać te elementy... :(

0

Popatrz np. Dla pierwszej pętli mamy:

u11=1-(1u11+0u12+0*u13)

Dla macierzy:

123
456
789

Zrobiłem

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