function r=testrouchecapelli(A,b)
  % Controlla se il rango della matrice associata al sistema
  % omogeneo e' uguale al rango di quella associata al sistema completo
  %   A*x = b

  % Esempio:
  %   A=rand(3,2)*rand(2,3); b=rand(3,1); testrouchecapelli(A,b)
  %     ===> 0

  % Esempio:
  %   A=rand(3,3); b=A*rand(3,1); testrouchecapelli(A,b)
  %     ===> 0

  r=(rank(A) == rank([A,b]));