ВМЛА лаб1 на дом оставленное

function X = gauss_simple(A,B)
  #Проверка на совпадение размерности
  if (rows(B)!=rows(A) || columns(B)!=1 || rows(A)!=columns(A))
    disp("\nWrong input\n");
  else
   Ag = [A B]; #делаем расширенную матрицу
   n = rows(Ag);
   m = columns(Ag);
   #приведем к треугольному виду:
   for k = 2:n
    for i = k:n
      temp = Ag(i,k-1);
      for j = 1:m
        Ag(i,j)=Ag(i,j)-Ag(k-1,j)*temp/Ag(k-1,k-1);
      endfor
     endfor
   endfor
   #раскручиваем матрицу:
   X = zeros(n,1);
   for i = 1:n
     X(n-i+1)=Ag(n-i+1,m);#записали самую праву строчку (идем снизу вверх)
     for j = m-i+1:m-1 #идем по одной строчке, вычитая слева направо
       X(n-i+1)=X(n-i+1)-Ag(n-i+1,j)*X(j);
     endfor
     X(n-i+1)=X(n-i+1)/Ag(n-i+1,m-i);
   endfor
  endif
endfunction