/* [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: D1068
Quelle: Johann Weilharter,
             Stochastik
******************************************

Für gegebene Werte von n und p sind
a) die Dichtefunktion der Binomialverteilung
b) die Binomialverteilung
c) der Erwartungswert
d) die Varianz und
e) die Streuung der Binomialverteilung
zu bestimmen.
   [wxMaxima: comment end   ] */

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

/* [wxMaxima: comment start ]
UNTERPROGRAMM LADEN
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
load(distrib);
/* [wxMaxima: input   end   ] */

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

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

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

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

/* [wxMaxima: input   start ] */
Dichte(x):=pdf_binomial(x,n,p);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Verteilung(x):=cdf_binomial(x,n,p);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
D:makelist(Dichte(x),x,0,n);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
D:D,numer;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
V:makelist(Verteilung(x),x,0,n);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
V:V,numer;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Erwartungswert:mean_binomial(n,p);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Erwartungswert:Erwartungswert,numer;
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: input   start ] */
Varianz:var_binomial(n,p);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Varianz:Varianz,numer;
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: input   start ] */
Streuung:std_binomial(n,p);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Streuung:Streuung,numer;
/* [wxMaxima: input   end   ] */

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

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

/* [wxMaxima: input   start ] */
print("Dichtefunktion =",D)$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
print("Verteilung =",V)$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
print("Erwartungswert =",Erwartungswert)$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
print("Varianz =",Varianz)$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
print("Streuung =",Streuung)$
/* [wxMaxima: input   end   ] */

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