set term gif set xrange [5.0:6.0] set nokey f0(x,m,s,N)=exp(-N*(x-m)*(x-m)/(2.0*s))*sqrt(N/(2*pi*s)) A=5.0 B=6.0 m=2.0*(B-A)/3.0+A s=(B-A)*(B-A)/18 f(x)=2.0*(x-A)/((B-A)*(B-A)) set pointsize 0.5 set title "基礎となる密度関数" set yrange [0.0:3.0] set output "data/dat-0.gif" plot 'data/dat-0' , f(x) set yrange [0.0:15.0] set output "data/dat-N1.gif" set title "中心極限定理のシミュレーション (N=1)" #plot 'data/dat-N1' plot 'data/dat-N1', f0(x,m,s,1) set output "data/dat-N2.gif" set title "中心極限定理のシミュレーション (N=2)" #plot 'data/dat-N2' plot 'data/dat-N2', f0(x,m,s,2) set output "data/dat-N3.gif" set title "中心極限定理のシミュレーション (N=3)" #plot 'data/dat-N3' plot 'data/dat-N3', f0(x,m,s,3) set output "data/dat-N5.gif" set title "中心極限定理のシミュレーション (N=5)" #plot 'data/dat-N5' plot 'data/dat-N5', f0(x,m,s,5) set output "data/dat-N10.gif" set title "中心極限定理のシミュレーション (N=10)" #plot 'data/dat-N10' plot 'data/dat-N10', f0(x,m,s,10) set output "data/dat-N20.gif" set title "中心極限定理のシミュレーション (N=20)" #plot 'data/dat-N20' plot 'data/dat-N20', f0(x,m,s,20) set output "data/dat-N50.gif" set title "中心極限定理のシミュレーション (N=50)" #plot 'data/dat-N50' plot 'data/dat-N50', f0(x,m,s,50)