#!/usr/bin/gnuplot -persist
#
#    
#    	G N U P L O T
#    	Linux version 3.7
#    	patchlevel 0
#    	last modified Thu Jan 14 19:34:53 BST 1999
#    
#    	Copyright(C) 1986 - 1993, 1998, 1999
#    	Thomas Williams, Colin Kelley and many others
#    
#    	Type `help` to access the on-line reference manual
#    	The gnuplot FAQ is available from
#    		<http://www.uni-karlsruhe.de/~ig25/gnuplot-faq/>
#    
#    	Send comments and requests for help to <info-gnuplot@dartmouth.edu>
#    	Send bugs, suggestions and mods to <bug-gnuplot@dartmouth.edu>
#    
# set terminal x11 0
# set output
set noclip points
set clip one
set noclip two
set bar 1.000000
set border 31 lt -1 lw 1.000
set xdata
set ydata
set zdata
set x2data
set y2data
set boxwidth
set dummy x,y
set format x "%g"
set format y "%.2f"
set format x2 "%g"
set format y2 "%g"
set format z "%g"
set angles radians
set nogrid
set key title ""
set key right top Right noreverse box linetype -2 linewidth 1.000 samplen 4 spacing 1 width 0
set nolabel
set noarrow
set nolinestyle
set nologscale
set offsets 0, 0, 0, 0
set pointsize 3
set encoding default
set nopolar
set noparametric
set view 60, 30, 1, 1
set samples 100, 100
set isosamples 10, 10
set surface
set nocontour
set clabel '%8.3g'
set mapping cartesian
set nohidden3d
set cntrparam order 4
set cntrparam linear
set cntrparam levels auto 5
set cntrparam points 5
set size ratio 0 1,1
set origin 0,0
set data style points
set function style lines
set xzeroaxis lt -2 lw 1.000
set x2zeroaxis lt -2 lw 1.000
set yzeroaxis lt -2 lw 1.000
set y2zeroaxis lt -2 lw 1.000
set tics in
set ticslevel 0.5
set tics scale 1
set mxtics default
set mytics default
set mx2tics default
set my2tics default
set xtics border mirror norotate autofreq 
set ytics border mirror norotate 0.05
set ztics border nomirror norotate autofreq 
set nox2tics
set noy2tics
set title "Face-Centered Cubic Tight-Binding Platinum" 0.000000,0.000000
set timestamp "" bottom norotate 0.000000,0.000000
set rrange [ * : * ] noreverse nowriteback  # (currently [-0:10] )
set trange [ * : * ] noreverse nowriteback  # (currently [-5:5] )
set urange [ * : * ] noreverse nowriteback  # (currently [-5:5] )
set vrange [ * : * ] noreverse nowriteback  # (currently [-5:5] )
set xlabel "Volume/Atom (Bohr^3)" 0.000000,0.000000
#set x2label "" 0.000000,0.000000  ""
set timefmt "%d/%m/%y\n%H:%M"
set xrange [ * : * ] noreverse nowriteback  # (currently [-10:10] )
#set x2range [ * : * ] noreverse nowriteback  # (currently [-10:10] )
set ylabel "Energy/Atom (Ry)" 0.000000,0.000000
#set y2label "" 0.000000,0.000000  ""
set yrange [ * : * ] noreverse nowriteback  # (currently [-10:10] )
#set y2range [ * : * ] noreverse nowriteback  # (currently [-10:10] )
set zlabel "" 0.000000,0.000000
set zrange [ * : * ] noreverse nowriteback  # (currently [-10:10] )
set zero 1e-08
set lmargin -1
set bmargin -1
set rmargin -1
set tmargin -1
set locale "C"
min(a,b) = a<b ? a:b
max(a,b) = a>b ? a:b
step(a) = a>0 ? 1:0
e(eo,vo,ko,kop,v) = eo + 1.125*ko*vo*t(vo,v)*t(vo,v)* (1.0 + 0.5*(kop-4.0)*t(vo,v))
t(vo,v) = (vo/v)**.666666666666666667 - 1.0
vo = 100
eo = -0.001
ko = 0.01
kop = 4
fit e(eo,vo,ko,kop,x) "SKENG" using 2:4 via eo,vo,ko,kop
ao = (4*vo)**.3333333333333333333333
kogpa = 14710.5164*ko
plot "SKENG" using 2:4 t "TB Points" w p 1 1,e(eo,vo,ko,kop,x) t "Birch Fit" w l 3
print "Equilibrium volume = ",vo," Bohr^3"
print "Equilibrium lattice constant = ",ao," Bohr"
print "Equilibrium energy = ",eo," Ry"
print "Equilibrium bulk modulus = ",kogpa," GPa"
save "ptplot.gnu"
pause -1 "Hit Enter to Exit"
#    EOF