## *** 描画の設定 ***
set terminal windows enhanced
set xlabel "Volume (Bohr^3 / Pd_{1-x}Rh_x)"
set ylabel "Energy (Ry)"
set format y "%8.3f"

## *** グラフのタイトル ***
set title "fcc PdRh_{100}"

## *** Birch-Murnghan の状態方程式 ***
vm(x)=9*v0*b/16.0*(dbdp*(((v0/x)**(2/3.0)-1)**3)+(((v0/x)**(2.0/3.0)-1)**2.0)*(6-4*((v0/x)**(2.0/3.0))))+e0
v0=115 ; e0=-10000 ; dbdp=4 ; b=0.01 ;
fit [95:110] vm(x) "./PdRh_100.txt" u (($1**3)/4.0):2 via v0,e0,dbdp,b

## *** グラフの描画 ***
plot './PdRh_100.txt' u (($1**3)/4.0):2 w p notitle,\
      vm(x) notitle

## *** フィッティングデータの書き出し ***
save variables "fitting_PdRh_100.txt"

## *** 画像の書き出し ***
set terminal pngcairo color enhanced size 520,390
set output "PdRh_100.png"
rep