with(stats):
with(plots):
with(statplots):
with(describe): Warning, these names have been redefined: anova, describe, fit, importdata, random, statevalf, statplots, transform
Warning, these names have been redefined: boxplot, histogram, scatterplot, xscale, xshift, xyexchange, xzexchange, yscale, yshift, yzexchange, zscale, zshift
Warning, these names have been redefined: coefficientofvariation, count, countmissing, covariance, decile, geometricmean, harmonicmean, kurtosis, linearcorrelation, mean, meandeviation, median, mode, moment, percentile, quadraticmean, quantile, quartile, range, skewness, standarddeviation, sumdata, variance
#niggli:
count:=proc(ll, x)
local i,r;
r:=0:
for i from 0 to nops(ll) do
if (op(i,ll) = x) then r:=r+1: end if:
end do:
return r:
end proc:
Total:=proc(ll)
local i,r,x;
r:=0:
for i from 1 to nops(ll) do
#x:=op(
r:=r+op(i,ll):
end do:
return r:
end proc:
roll:=rand(1..6):
kmax:=10:
w:=[seq( [seq(roll(),i=1..12) ], k=1..kmax)]; NiM+SSJ3RzYiNyw3LiIiJiIiJEYoIiInRihGKCIiIkYrIiIjRisiIiVGLTcuRixGKEYpRi1GK0YrRilGKkYqRitGLEYqNy5GKkYtRipGLEYoRixGKUYtRihGKEYoRik3LkYqRilGKkYqRihGK0YqRihGK0YqRilGKjcuRi1GK0YrRilGLUYpRipGLUYpRitGLUYoNy5GKUYpRilGLEYoRilGKkYsRipGLUYrRig3LkYtRixGK0YsRi1GK0YoRihGK0YpRihGKzcuRipGKUYoRixGLUYtRilGLEYsRixGLEYoNy5GKEYrRilGKUYtRixGKkYoRitGKkYpRik3LkYpRitGLUYqRipGKUYrRitGKkYsRitGLA== dAugenzahl:=[ seq( Total( op(i,w) ), i=1..kmax ) ];
d6er:=[ seq( count(op(i,w),6), i=1..kmax ) ];
sAugenzahl:=evalf( [mean(dAugenzahl), standarddeviation(dAugenzahl)] );
s6er:=evalf( [mean(d6er), standarddeviation(d6er)] );
pa:=listplot(dAugenzahl):
ps:=listplot(d6er):
pm:=plot(mean(dAugenzahl),x=1..10):
display([pa,ps,pm]); NiM+SStkQXVnZW56YWhsRzYiNywiI1UiI1MiI10iI2EiI1IiI1YiI01GKEYnIiNP NiM+SSVkNmVyRzYiNywiIiIiIiQiIiMiIidGJ0YpIiIhRidGKUYo NiM+SStzQXVnZW56YWhsRzYiNyQkIiNVIiIhJCIrcTRrNGQhIio= NiM+SSVzNmVyRzYiNyQkIisrKysrQCEiKiQiKyVRdHpkIkYp LSUlUExPVEc2KC0lJ0NVUlZFU0c2IzcsNyQkIiIiIiIhJCIjVUYsNyQkIiIjRiwkIiNTRiw3JCQiIiRGLCQiI11GLDckJCIiJUYsJCIjYUYsNyQkIiImRiwkIiNSRiw3JCQiIidGLCQiI1ZGLDckJCIiKEYsJCIjTUYsNyQkIiIpRixGMjckJCIiKkYsRi03JCQiIzVGLCQiI09GLC1GJjYjNyw3JEYqRio3JEYwRjU3JEY1RjA3JEY6RkQ3JEY/Rio3JEZERjA3JEZJJEYsRiw3JEZORio3JEZRRjA3JEZURjUtRiY2JDdTRik3JCQiMyMqKioqKlwoUXVoPiIhIzxGLTckJCIzLytdN3RZJ29PIkZmb0YtNyQkIjMxKytEXlAjKWU6RmZvRi03JCQiMy0rK3ZlXzBfPEZmb0YtNyQkIjMhKioqXGlxJG9WJT5GZm9GLTckJCIzKioqKlw3VnNtQUBGZm9GLTckJCIzKikqKlw3KVImRzJCRmZvRi03JCQiMy0rXTdleEApXCNGZm9GLTckJCIzJSkqKlxpJnpQJilvI0Zmb0YtNyQkIjNBKytdWmBJJSlHRmZvRi03JCQiMyMqKioqXFBeUG4wJEZmb0YtNyQkIjNHKytdI0hiM0QkRmZvRi03JCQiMyUpKioqKlxQLHhYTUZmb0YtNyQkIjNFKytdMm5nTE9GZm9GLTckJCIzJikqKlw3My49L1FGZm9GLTckJCIzLysrXTwpM3ErJUZmb0YtNyQkIjNbKysrcTYkKXlURmZvRi03JCQiM1ErXTc4OXF5VkZmb0YtNyQkIjN1KioqKipcVz9jYiVGZm9GLTckJCIzVStdN2onRyhcWkZmb0YtNyQkIjNYK11QKmVsWCRcRmZvRi03JCQiM3MrK0RKTlVGXkZmb0YtNyQkIjMyK11pRHRfL2BGZm9GLTckJCIzWyoqKlxQcGRiXCZGZm9GLTckJCIzIyoqKlw3ZVgpUnAmRmZvRi03JCQiM0srXShvczpuJ2VGZm9GLTckJCIzdSoqKlw3N3FLMCdGZm9GLTckJCIzWSoqKioqXDEqKmZDJ0Zmb0YtNyQkIjN0KioqXGl0WVhWJ0Zmb0YtNyQkIjMnKioqXDcuaihwaCdGZm9GLTckJCIzbCoqKlxQQkwmPm9GZm9GLTckJCIzUCoqKioqXGtSOisoRmZvRi03JCQiM2UrK11QLihlPihGZm9GLTckJCIzZSoqXDdHRyc+UChGZm9GLTckJCIzXysrXUsleVdjKEZmb0YtNyQkIjNTKipcNzgxaVh4RmZvRi03JCQiM24qKlxpJlFtXCR6RmZvRi03JCQiMzErK10oWyczPyIpRmZvRi03JCQiM3EpKlw3eSsqUUopRmZvRi03JCQiMywsKytxZmErJilGZm9GLTckJCIzPSoqKlwoeSZHOXApRmZvRi03JCQiM18rXTckZUkyKSkpRmZvRi03JCQiM2EqKioqKlxZelkwKkZmb0YtNyQkIjNRKioqXFBeV1NEKkZmb0YtNyQkIjMqMysrKyoqZUJWKkZmb0YtNyQkIjNMKipcNzglekNpKkZmb0YtNyQkIjN2KipcKG8iKltXISkqRmZvRi03JEZURi0tJSZDT0xPUkc2JiUkUkdCRyRGVSEiIiRGLEZneEZoeC0lK0FYRVNMQUJFTFNHNiRRIng2IlEhRl15LSUlRk9OVEc2JCUqSEVMVkVUSUNBR0ZVLSUlVklFV0c2JDtGZnhGVDskISQzIiEiIyQiJTNiRmp5