clear; clf;

// Initialisation

resultat=zeros(12, 12, 12);

// Calcul : toutes les valeurs sont calculées
// mais toutes ne sont pas affichées

for imax = 4:2:12
   for jmax = imax:2:12
      for i = 1:imax
         for j = 1:jmax
            if j > i then resultat(imax, jmax, j) = resultat(imax, jmax, j) + 1;
            else resultat(imax, jmax, i) = resultat(imax, jmax, i) + 1; end;
         end
      end
   end
end

// Transformation dénombrement -> probabilité
// d4

somme44 = sum(resultat(4, 4, :));
for i = 1:4
   y44(i) = 100*resultat(4, 4, i)/somme44;
end

somme46 = sum(resultat(4, 6, :));
for i = 1:6
   y46(i) = 100*resultat(4, 6, i)/somme46;
end

somme48 = sum(resultat(4, 8, :));
for i = 1:8
   y48(i) = 100*resultat(4, 8, i)/somme48;
end

somme410 = sum(resultat(4, 10, :));
for i = 1:10
   y410(i) = 100*resultat(4, 10, i)/somme410;
end

somme412 = sum(resultat(4, 12, :));
for i = 1:12
   y412(i) = 100*resultat(4, 12 , i)/somme412;
end

// d6

somme66 = sum(resultat(6, 6, :));
for i = 1:6
   y66(i) = 100*resultat(6, 6, i)/somme66;
end

somme68 = sum(resultat(6, 8, :));
for i = 1:8
   y68(i) = 100*resultat(6, 8, i)/somme68;
end

somme610 = sum(resultat(6, 10, :));
for i = 1:10
   y610(i) = 100*resultat(6, 10, i)/somme610;
end

somme612 = sum(resultat(6, 12, :));
for i = 1:12
   y612(i) = 100*resultat(6, 12, i)/somme612;
end

// d8

somme88 = sum(resultat(8, 8, :));
for i = 1:8
   y88(i) = 100*resultat(8, 8, i)/somme88;
end

somme810 = sum(resultat(8, 10, :));
for i = 1:10
   y810(i) = 100*resultat(8, 10, i)/somme810;
end

somme812 = sum(resultat(8, 12, :));
for i = 1:12
   y812(i) = 100*resultat(8, 12, i)/somme812;
end

// d10

somme1010 = sum(resultat(10, 10, :));
for i = 1:10
   y1010(i) = 100*resultat(10, 10, i)/somme1010;
end

somme1012 = sum(resultat(10, 12, :));
for i = 1:12
   y1012(i) = 100*resultat(10, 12, i)/somme1012;
end

// d12

somme1212 = sum(resultat(12, 12, :));
for i = 1:12
   y1212(i) = 100*resultat(12, 12, i)/somme1212;
end

// tracé
// seules certaines courbes sont tracées

subplot(2, 2, 1)

plot2d(1:4, y44, rect = [1, 0, 12, 45], syle = 1)
plot2d(1:6, y46, rect = [1, 0, 12, 45], style = 2)
plot2d(1:8, y48, rect = [1, 0, 12, 45], style = 3)
plot2d(1:10, y410, rect = [1, 0, 12, 45], style = 4)
plot2d(1:12, y412, rect = [1, 0, 12, 45], style = 5)

xtitle("d4", "rés.", "%")

subplot(2, 2, 3)

plot2d(1:8, y48, rect = [1, 0, 12, 25], style = 1)
plot2d(1:8, y68, rect = [1, 0, 12, 25], style = 2)
plot2d(1:8, y88, rect = [1, 0, 12, 25], style = 3)
plot2d(1:10, y810, rect = [1, 0, 12, 25], style = 4)
plot2d(1:12, y812, rect = [1, 0, 12, 25], style = 5)

xtitle("d8", "rés.", "%")

subplot(2, 2, 2)

plot2d(1:12, [y412, y612, y812, y1012, y1212], rect = [1, 0, 12, 16])

xtitle("d12", "rés.", "%")

subplot(2, 2, 4)

xstring(0.5, 0.8, "Légende")

xsegs([0, 0, 0, 0.6, 0.6 ; 0.1, 0.1, 0.1, 0.7, 0.7], ...
   [0.6, 0.4, 0.2, 0.6, 0.4 ; 0.6, 0.4, 0.2, 0.6, 0.4], ...
   1:5)
xstring(0.15, 0.55, "d4")
xstring(0.15, 0.35, "d6")
xstring(0.15, 0.15, "d8")
xstring(0.75, 0.55, "d10")
xstring(0.75, 0.35, "d12")