fig = scf(1);
fig.figure_size = [630, 891];
clf;

x = [1:10, 10:10:100, 100:100:1000, 1000:1000:10000];
lnx = log(x);
lnxmin = min(lnx);
lnxmax = max(lnx);
nx = size(x, 'c');

F = [0.01:0.005:0.05, 0.03:0.01:0.1, 0.12:0.02:0.4, 0.45:0.05:0.95, 0.96:0.01:0.99];
y = log(-log(1-F));
ymin = min(y);
ymax = max(y);
ny = size(F, 'c');

for i = 1:ny
    yy = y(i); FF = F(i);
    xpoly([lnxmin, lnxmax],[yy, yy]);
    h1 = gce();
    h1.thickness = 0.5;
    xstring(-0.6, yy-0.08, string(FF));
    h2 = gce();
    h2.font_size = 0.3;
end
FF = 1-exp(-1);
yy = log(log(1-FF));
xpoly([lnxmin, lnxmax],[yy, yy])
h1 = gce();
h1.thickness = 0.5;
h1.line_style=3;

for i = 1:nx
    xx = lnx(i);
    xpoly([xx, xx], [ymin, ymax])
    h1 = gce();
    h1.thickness = 0.5;
end
a = gca();
a.data_bounds=[0, ymin ; 10, ymax];