function linsol(A,b)
  % Risolve il sistema lineare
  %   A x = b

  if not(testrouchecapelli(A,b))
    disp 'Il sistema lineare non ha soluzione';
    return
  end

  if not(testrankmax(A))
    m=size(A,2);
    C=[A,b];
    C=rankmax(C);
    A=C(1:i,m);
    b=C(1:i,m+1);
  else

  mostralinsol(A,b);