Octave / SLAU

f=fopen('SLAU1.txt','rt+');
d = dlmread(f)
n=d(1,1)
a=d(2 : n+1, 1 : n)
b=d(2 : n+1, n + 1)
if (det(a)~=0)
  x=a\b
  disp('Проверка')
  s=a*x;
  if(s==b)
  disp('Верно')
  else
  disp('Неверно')
  endif
else
  disp('Нет решения')
endif
fclose(f);