Boczek11

Próbowałem przedstawić metodę eliminacji Gaussa w mathematice, lecz coś mi się "psuje" tzn. wynik dla x3 jest poprawny natomiast błędne wyniki wychodzą w przypadku x1 jaki i x2. Czy ktoś mógłby mi pomóc?? Gauss[a1_, b1_] := Module[{n, a, b, p, x}, a = a1; b = b1; n = Length[b1]; x = Table[0, {i, 1, n}]; Do[Do[p = a[[i, j]]/a[[j, j]]; b[[i]] = b[[i]] - p*b[[j]]; a[[i]] = a[[i]] - p*a[[j]], {i, j + 1, n}], {j, 1, n - 1}]; x[[n]] = b[[n]]/a[[n, n]]; Do[x[[n - i]] = (b[[n - i]] - Do[x[[n - k]]*a[[n - i, n - k]]{k, 0, i - 1, 1}])/a[[n - i, n - i]], {i, 1, n - 1}]; Return[x]; ]Gauss[{{1/4, 1/5, 1/6}, {1/3, 1/4, 1/5}, {1/2, 1, 2}}, {9, 8, 8}]

Dziękuje bardzo za pomoc.

Witam, Mój problem wygląda następująco. Musze napisać program na rozwiązywanie układu liniowego metodą eliminacji Gaussa w języku C++. Znam algorytm(korzystałem z książki autorstwa D.Kincaid oraz W. Cheney "Analiza numeryczna"),lecz przekształcenie go w program jest dla mnie czarną magią .Bardzo proszę o jakąkolwiek pomoc.