function p = probdev(P,F,m,n,delta)
Ef = F*P;
del = abs(F-Ef);

ii = find(del>=delta*n);
p = sum(P(ii));
return