#!/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 "%.3f"
set format y "%.3f"
set format x2 "%.3f"
set format y2 "%.3f"
set format z "%.3f"
set angles radians
set nogrid
set key title ""
set key left 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 ticscale 1 0.5
set mxtics default
set mytics default
set mx2tics default
set my2tics default
set xtics border mirror norotate autofreq 
set ytics border mirror norotate 0.002
set ztics border nomirror norotate autofreq 
set nox2tics
set noy2tics
set title "C_{44} for bcc W at a = 5.93 Bohr" 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 "x^2 (= e_6^2)" 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 "E (Ry)" 0.000000,0.000000  ""
set y2label "" 0.000000,0.000000  ""
set yrange [ 0 : 0.008 ] noreverse nowriteback
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(e0,c2,c4,c6,x) = e0 + 0.5*v*x*(c2+x*(c4+x*c6))
a = 5.93
v = 104.2639285
e0 = 0.000538217706293788
cbohr = 0.0129143246974814
c4 = 0.156781800910638
c6 = -8.31741418396513
c44 = 189.976385257225
plot "SKENG" using 1:4 t "TB Points" w p,e(e0,cbohr,c4,c6,x) t "Fit" w l 3
## fit e(e0,cbohr,c4,c6,x) "SKENG" using 1:4 via e0,cbohr,c4,c6
#    EOF