gset output "serie_3_plot.ps";
gset term post enhanced;
gset key right bottom;
gset title "Poissonverteilung Asbestfasern";
gset xlabel "S={/Symbol S} _i^n{X_i}";
gset ylabel "Verteilungsfunktion";
k=[0:1:30];
y_15 = poisson_cdf(k,15);
y_poisson = [0,y_15(1:30)];
y_poisson = y_15 .- y_poisson;
plot(k,y_15,['o',';kumulative
Verteilung;'],[0,30],[0.05,0.05],[';5%;'],[7,7],[0,1],[';7;'],k,y_poisson,['*',';Verteilung;']);
printf("y_15:\n",k,y_poisson,'x');
y_15(1:15)
y_10 = poisson_cdf(k,10);
printf("y_10:\n");
y_10(1:15)
printf("Letzter Wert =< 0.05: possion(5,10) = %f\n",y_10(5));
printf("Fehler 2. Art: 1-poisson_cdf(5,10) = %f\n",1-poisson_cdf(5,10));

y_15:
ans =

 Columns 1 through 6:

   3.0590e-07   4.8944e-06   3.9308e-05   2.1138e-04   8.5664e-04   2.7924e-03

 Columns 7 through 12:

   7.6319e-03   1.8002e-02   3.7446e-02   6.9854e-02   1.1846e-01   1.8475e-01

 Columns 13 through 15:

   2.6761e-01   3.6322e-01   4.6565e-01

y_10:
ans =

 Columns 1 through 6:

   4.5400e-05   4.9940e-04   2.7694e-03   1.0336e-02   2.9253e-02   6.7086e-02

 Columns 7 through 12:

   1.3014e-01   2.2022e-01   3.3282e-01   4.5793e-01   5.8304e-01   6.9678e-01

 Columns 13 through 15:

   7.9156e-01   8.6446e-01   9.1654e-01

Letzter Wert =< 0.05: possion(5,10) = 0.029253
Fehler 2. Art: 1-poisson_cdf(5,10) = 0.932914