j'ai simplement recherche les quadruplets (a,b,c,d) verifiant par exemple 20<=a<=b<=c<=d<=70
Ici avec une fonction en Xcas et en français.
fonction Resolution()
local L,a,b,c,d,D,m,s;
L:=[];
pour a de 20 jusque 70 faire
pour b de a jusque 70 faire
pour c de b jusque 70 faire
pour d de c jusque 70 faire
D:=[30.0,35,32,31,18,25,29,50,36,44,32,26,20,24,46,45,a,b,c,d];
m:=moyenne(D);
s:=ecart_type_population(D);
si m==33.05 et abs(s-8.941182072)<=10^-9 alors
L.append([a,b,c,d]);
fsi
fpour
fpour
fpour
fpour
retourne L
ffonction:;