/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/

/* [wxMaxima: title   start ]
Binomialverteilung
   [wxMaxima: title   end   ] */

/* [wxMaxima: section start ]
Problembeschreibung
   [wxMaxima: section end   ] */

/* [wxMaxima: comment start ]
******************************************
Dokumentnummer: D1070
Quelle: Johann Weilharter,
             Stochastik
******************************************

Von einer Binomialverteliung kennt man die Parameter n und p.
Zu erstellen sind:
a) die Verteilung
b) die kumulierte Verteilung
c) eine Grafik der Verteilung
d) ein definierter Teilbereich der Verteilung
   [wxMaxima: comment end   ] */

/* [wxMaxima: section start ]
Problemlösung
   [wxMaxima: section end   ] */

/* [wxMaxima: comment start ]
EINGABE
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
n:read("Gib den Parameter n ein");
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
p:read("Gib die Wahrscheinlichkeit p ein");
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
VERARBEITUNG
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
W(k):=binomial(n,k)*p**k*(1-p)**(n-k);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
kumW(k):=sum(W(i),i,0,k);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
AUSGABE
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
for j:0 thru n do
    print(j,?truncate(W(j)*1000+0.5)/1000.0,?truncate(kumW(j)*1000+0.5)/1000.0);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
Die Wahrscheinlichkeitsverteilung kann grafisch
dargestellt werden:
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
wxplot2d([W(x)], [x,0,n])$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
Häufig benötigt man nur Teilbereiche der gesamten
Verteilung:
   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
EINGABE
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
a:read("Gib die Untergrenze ein");
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
b:read("Gib die Obergrenze ein");
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
AUSGABE
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
for j:a thru b do
    print(j,floor(W(j)*1000+0.5)/1000.0,floor(kumW(j)*1000+0.5)/1000.0);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
summe:sum(W(k),k,a,b),numer;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
summe:floor(summe*1000+0.5)/1000.0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
print("Summe von",a,"bis",b,"=",summe)$
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$