function[xm,fv]=POS(fitness,N,c1,c2,w,M,D)