[richy]

Und schon fertig :

restart;
N_range:=36; # Klassenintervalle
N:=N_range*10; # Verhulst Schritte > Intervalle
N_vorlauf:=10;
N_big:=100; # Anzahl Parameter r Werte
r_min:=3; # r-range
r_max:=3.9;


buch[1]:=17.40:
buch[2]:=9.78:
.........
nor:=0;
for i from 1 to 27 do nor:=nor+buch[i]; od:
for i from 1 to 27 do buch[i]:=buch[i]/nor;od:

# Grosse v-Schleife fuer alle r
################################
for v from 1 to N_big do

# Verhulst Parameter
r[v]:=evalf(r_min+(v-1)/(N_big-1)*(r_max-r_min)):
for i from 1 to N_range do
b[i]:=0; # Anzahl Treffer initialisieren
od:

####### Vorlauf Verhulst
s:=0.1:
for i from 1 to N_vorlauf do s:=r[v]*s*(1-s);od:
####### Verhulst Schleife
for i from 1 to N do
s:=r[v]*s*(1-s);
j:=floor(s*(N_range-1)+1);
b[j]:=b[j]+1;
od:
#b sortieren nach Wert(Bubblesort)
for i from 1 to N_range-1 do
for j from i to N_range do
if (b[i]<b[j]) then
mem:=b[i];b[i]:=b[j];b[j]:=mem;
mem:=a[i];a[i]:=a[j];a[j]:=mem;
fi;
od: od:
Nmax:=27:
#Normieren der Messwerte
s:=0:
for i from 1 to Nmax do s:=s+b[i]; od:
for i from 1 to Nmax do b[i]:=b[i]/s; od:

# Gauss Fehlerintegral
gauss[v]:=0:
for i from 1 to Nmax do

gauss[v]:=gauss[v]+(buch[i]-b[i])**2; # Gaussches Fehlerintegral
od:
od: # Ende grosse Schleife
################################

druck1:=seq([r[v],gauss[v]],v=1..N_big):
plot([[druck1]]);