// L'exécution de ce script nécessite le module CASCI

// déciles

Y = 0.1:0.1:0.9;
X = idfnormal(Y, 0, 1);
n = size(Y, "*");

// lois normales

x = -2:0.1:2

y1 = pdfnormal(x, -1, 1);
y2 = pdfnormal(x, 0, 1);
y3 = pdfnormal(x, 0, 0.5);
y4 = pdfnormal(x, 0, 2);

scf(0)
clf;

// tracé des positions

for i = 1:n
    xpoly([X(i), X(i)], [0, 0.41], "lines");
    h = gce();
    h.line_style=2;
end
plot2d(x', [y1', y2'], style=[1 2])

scf(1)
clf;

// tracé de spositions

for i = 1:n
    xpoly([X(i), X(i)], [0, 0.81], "lines");
    h = gce();
    h.line_style=2;
end
plot2d(x', [y3', y2', y4'], style=[1 2 3])