function [neighboor] = nearest(v, data)
// be sure v is a row vector
  [m,n]=size(v);
  if n<m then v=v'; end
  sub_data = data(:, 2:$);
  [m,n] = size(sub_data);
  I = ones(m, 1);
  M = I*v;
  norme = sum((sub_data-M).^2, 'c');
  [val,i]=min(norme);
  neighboor = i;
endfunction


function [p,erreur] = classe(test, data)
  [m,n]=size(test);
  erreur=[];
  status = 0
  for i=1:m
    neighboor=nearest(test(i,2:$), data);
    if (data(neighboor,1)~=test(i,1)) then
      status=status+1;
      erreur=[erreur;i,neighboor];
    end
    if modulo(10*i,m)==0 then
      printf("%d \% des images traitees\n", 100*i/m);
    end
  end
  p = status/m;
endfunction