To je programiranje. To je cist numericki problem i neko treba da ga isprogramira. Ta dva idu zajedno ko hleb i voda.
Ovako postupak ti je sledeci.
X1=b1/a11
x2=(b2-a21*x1)/a22
evo ti algoritam, pa ga ti prebaci u koji god oces jezik
Ovako postupak ti je sledeci.
X1=b1/a11
x2=(b2-a21*x1)/a22
evo ti algoritam, pa ga ti prebaci u koji god oces jezik
Kod:
function [x, flag] = gauss(a, b)
% ulazna matrica
% ulazni vektor kolona
% x rezultat
% flag indikator greske
% 0 - sve je u redu
% 1 - greska
[n, m] = size(a);
% n je broj vrsta
% m je broj kolona
% ovde bi trebalo proveriti da li su n i m isti
nB = length(b);
% ovde bi trebalo proveriti da je nB isto sto i n
% svodjenje na trougaonu matricu
for k=1:n
for i=k+1:n
m = a(i,k)/a(k,k);
for j=k:n
a(i,j) = a(i,j)-m*a(k,j);
end
b(i) = b(i)-m*b(k);
end
end