function maxInd= maxInd(prob,NP)
%return NP most likely particles
[n,m]=size(prob);
prob=reshape(prob,n*m,1);
x=[(1:n*m)' -prob];
x=sortrows(x,2);
maxInd=x(1:NP,1);
