% show soln
global A C lnFinv finv k n result
n=256; k=4;
load lnFinv
finv=exp(lnFinv);
disp("showing result"); fflush(1);
result=zeros(n);
for i=1:n
    kq=k *finv(i);
    j=1; while ((finv(j)<kq)&&(j<256)) ; j=j+1; end%while
    k=[i j]
    result(j,i)=255;
end%for
imshow(result);

