# this works for Octave only

a = randn(30,1); 
b = randn(30,1)+2;

x=(-3:0.1:5); 
pa = normal_pdf(x,0,1); 
pb=normal_pdf(x,2,1);

plot(a,zeros(1,30),'1*',b,zeros(1,30),'3x',x,pa,'1-',x,pb,'3-');