set terminal png                   #出力をpngファイルに指定
set output "fig_test7.png"         #pngファイル名指定
#
set title "Probability oval"
set size square
set parametric        #媒介変数使用
set dummy t           #媒介変数tを明示的に指定
set trange [0:2*pi]   #tの範囲指定
set samples 360       #tのサンプリング数を360に設定
#
xm1=0.0;sx1=1.0;ym1=0.0;sy1=1.0;r1=0.7;p1=0.6
f1(t)=xm1+sx1*sqrt(-2*(1-r1*r1)*log(1-p1)/(1-2*r1*sin(t)*cos(t)))*cos(t)
g1(t)=ym1+sy1*sqrt(-2*(1-r1*r1)*log(1-p1)/(1-2*r1*sin(t)*cos(t)))*sin(t)
#
xm2=1.0;sx2=1.0;ym2=1.0;sy2=1.0;r2=0.7;p2=0.6
f2(t)=xm2+sx2*sqrt(-2*(1-r2*r2)*log(1-p2)/(1-2*r2*sin(t)*cos(t)))*cos(t)
g2(t)=ym2+sy2*sqrt(-2*(1-r2*r2)*log(1-p2)/(1-2*r2*sin(t)*cos(t)))*sin(t)
#
xm3=2.0;sx3=1.0;ym3=2.0;sy3=1.0;r3=0.7;p3=0.6
f3(t)=xm3+sx3*sqrt(-2*(1-r3*r3)*log(1-p3)/(1-2*r3*sin(t)*cos(t)))*cos(t)
g3(t)=ym3+sy3*sqrt(-2*(1-r3*r3)*log(1-p3)/(1-2*r3*sin(t)*cos(t)))*sin(t)
#
plot \
f1(t),g1(t), \
f2(t),g2(t), \
f3(t),g3(t)
#