n = 100;
X = rand(1, n, "norm"); // tirage aléatoire de n valeurs
// avec une loi normale centrée réduite

p = 0.1:0.1:0.9;
q = 1 - p;

qexp = quantile(X, p); // déciles expérimentaux

for i = 1:size(p, "*") // déciles théorique pour une loi normale centrée réduite
    qth(i) = cdfnor("X", 0, 1, p(i), q(i));
end

scf(0);
clf();
qqplot(qth', qexp); // diagramme quantile-quantile
xtitle("diagramme quantiles-quantiles", "quantiles théoriques",...
    "quantiles expérimentaux");