# set terminal pngcairo  transparent enhanced fontscale 1.0 size 600, 400 
# set output 'fit.5.png'
set bar 1.000000 front
set style circle radius graph 0.02, first 0.00000, 0.00000 
set style ellipse size graph 0.05, 0.03, first 0.00000 angle 0 units xy
set key bmargin center horizontal Right noreverse enhanced autotitle nobox
set style textbox transparent margins  1.0,  1.0 border
unset logscale
unset paxis 1 tics
unset paxis 2 tics
unset paxis 3 tics
unset paxis 4 tics
unset paxis 5 tics
unset paxis 6 tics
unset paxis 7 tics
set title "initial parameters for realistic model function" 
set xlabel "Temperature T (°C)" 
set ylabel "Density (g/cm^3)" 
set paxis 1 range [ * : * ] noreverse nowriteback
set paxis 2 range [ * : * ] noreverse nowriteback
set paxis 3 range [ * : * ] noreverse nowriteback
set paxis 4 range [ * : * ] noreverse nowriteback
set paxis 5 range [ * : * ] noreverse nowriteback
set paxis 6 range [ * : * ] noreverse nowriteback
set paxis 7 range [ * : * ] noreverse nowriteback
set colorbox vertical origin screen 0.9, 0.2, 0 size screen 0.05, 0.6, 0 front  noinvert bdefault
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)
myencoding = "utf8"
y0 = 1.06865074176486
m = -0.000943519626924529
GPFUN_l = "l(x) = y0 + m*x"
x = 0.0
FIT_CONVERGED = 1
FIT_NDF = 35
FIT_STDFIT = 0.148686593680091
FIT_WSSR = 0.773769609906597
FIT_P = 1.0
FIT_NITER = 6
y0_err = 0.00148404909012443
m_err = 3.15383626024729e-05
ml = -0.0005
mh = -0.0005
dens_Tc = 1.02
Tc = 45.0
g = 1.0
b = 0.01002
GPFUN_high = "high(x) = mh*(x-Tc) + dens_Tc"
GPFUN_lowlin = "lowlin(x)  = ml*(x-Tc) + dens_Tc"
GPFUN_curve = "curve(x) = b*tanh(g*(Tc-x))"
GPFUN_density = "density(x) = x < Tc ? curve(x)+lowlin(x) : high(x)"
## Last datafile plotted: "lcdemo.dat"
plot 'lcdemo.dat' using 1:2:5 with errorbars, density(x)
## fit l(x) 'lcdemo.dat' using 1:2:5 yerr via y0, m