function vizpdf(F,P)

Ef = F*P;
del = abs(F-Ef);

del = round(del*100)/100;

D = unique(del);
pD = zeros(size(D));
for i=1:length(D)
  jj = find(del==D(i));
  pD(i) = sum(P(jj));
end
Ef
D
pD
plot(D,pD,'.');
return