function [res] = projection(M, vectpr, k)
  c=M*vectpr;
  c(:,k+1:784)=0;
  res=c*vectpr';
endfunction