# set terminal png transparent nocrop enhanced font arial 8 size 420,320 
# set output 'fit.15.png'
set key right below Right noreverse noinvert enhanced box linetype -2 linewidth 1.000 samplen 4 spacing 1 width 0 height 0 autotitles
set title "fitted to realistic model function" 0.000000,0.000000  font ""
set xlabel "Temperature T  [deg Cels.]" 0.000000,0.000000  font ""
set ylabel "Density [g/cm3]" 0.000000,0.000000  font ""
l(x) = y0 + m*x
high(x) = mh*(x-Tc) + dens_Tc
lowlin(x)  = ml*(x-Tc) + dens_Tc
curve(x) = b*tanh(g*(Tc-x))
density(x) = x < Tc ? curve(x)+lowlin(x) : high(x)
y0 = 1.06865074176476
m = -0.000943519626922319
ml = -0.00100003889377453
mh = -0.00083126648503697
dens_Tc = 1.02497044445617
Tc = 46.0899249399587
g = 3.85580365358305
b = 0.00153901455759968
plot 'lcdemo.dat', density(x)
## fit density(x) 'lcdemo.dat' via 'start.par'