% plot the certainty function from Finv
% use: certaintyfromFinv(Finv);

function plotcertaintyfromFinv(Finv);

gset terminal x11
Uncertainty=diff(Finv);
Certainty=1../Uncertainty;
Certainty=interp_odd_to_even(Certainty);

Certainty=real(Certainty);

%certainty=exp(Certainty); too sharp a peak
%plot(exp(Certainty./10)) % seems good
certainty=exp(Certainty./max(Certainty/exp(1))); % normalize on interval 0:1
certainty=certainty/max(certainty);
axis([0 255 0 max(certainty)*1.05]);
gset nokey
xlabel('pixel integer')
ylabel('relative certainty')
plot(certainty);

gset terminal postscript eps 35 # show it like it is
system('rm deleteme_pork.eps');
gset output 'deleteme_pork.eps';
replot
gset terminal x11

axis;
keyboard