Matvec Users’ Guide

50 CHAPTER 8. 2D/3D OBJECT-ORIENTED PLOTTING
8.3.3 Example 3
The following Matvec script draw the graph of normal pdf function for two sets of parameters.
D1 = StatDist("Normal",0,1);
D2 = StatDist("Normal",2,0.5);
x = -4:4:0.1;
xy1 = [x; D1.pdf(x)]’;
xy2 = [x; D2.pdf(x)]’;
xy1.save("Normal.dat1",1);
xy2.save("Normal.dat2",1);
curve1 = " ’Normal.dat1’ title ’N(0,1)’ with lines 1,";
curve2 = " ’Normal.dat2’ title ’N(2,0.5)’ with lines 4";
curve = curve1 + curve2;
P = Plotter();
P.set("nozeroaxis");
P.set("size 4/5.0, 2.5/3.0")
P.set("title ’Normal Probability Density Function (pdf)’ ");
P.set("xlabel ’x’");
P.set("ylabel ’pdf(x)’");
P.set("key -2.0,0.5");
P.plot(curve);
The graph is shown on page 90.
8.3.4 Example 4
The following Matvec script draw the graph of the pdf of b(25, 0.15).
D = StatDist("Binomial",25,0.15);
x = -0.5:19.5;
xy = [x; D.pdf([0:20])]’;
xy.save("Binomial.dat",1);
P = Plotter();
P.set("title ’Binomial Probability Density Function (pdf)’ ");
P.set("nozeroaxis");
P.set("format y ’%3.2f’ ");
P.set("xtics 0, 1.0, 13");
P.set("xrange [-0.5 : 12.5]");
P.set("size 4/5.0, 2.5/3.0")
P.set("xlabel ’k’");
P.set("ylabel ’pdf(k)’");
P.plot(" ’Binomial.dat’ title ’b(25,0.15)’ with steps");
The graph is shown on page 103.