正規分布ｘ - 経済統計練習帳

var point=new Array();
 
function sample(){
 
for(s=1;s<101;s++){    
point[s]=0;
}
 
for(s=1;s<10001;s++){
x1=f();
for(sx=1;sx<101;sx++){    
if(x1>0.06*sx-3)y1=sx;    
}
point[y1]=point[y1]+1;
}
 
 paint();
 
 
 
}
 
function paint(){
 
 var sub=new svgpaint();
 
str="<svg width=\"500\" height=\"500\">";
 
for(s=1;s<101;s++){    
str=str+sub.rect(5*s,500-point[s],5,point[s],"blue");
}
str=str+"</svg>";
 
$("#memo").html(str);   
 
}
 
 
 
 
function f(){
 
var pi=Math.PI;
 
var x=Math.random();
var y=Math.random();
 
var z1=Math.cos(2*pi*y);
var z2=Math.sin(2*pi*y);
 
var k1=-2*Math.log(x);    
 
var z1=Math.sqrt(k1)*z1;
var z2=Math.sqrt(k1)*z2;   
 
 
return z1;    
}