# set terminal pngcairo transparent enhanced fontscale 1.0 size 600, 400 # set output 'fit.5.png' set key bmargin center horizontal Right noreverse enhanced autotitle nobox set title "initial parameters for realistic model function" set xlabel "Temperature T (°C)" set xrange [ * : * ] noreverse writeback set x2range [ * : * ] noreverse writeback set ylabel "Density (g/cm^3)" set yrange [ * : * ] noreverse writeback set y2range [ * : * ] noreverse writeback set zrange [ * : * ] noreverse writeback set cbrange [ * : * ] noreverse writeback set rrange [ * : * ] noreverse writeback set colorbox vertical origin screen 0.9, 0.2 size screen 0.05, 0.6 front noinvert bdefault set fit nologfile brief errorvariables nocovariancevariables errorscaling prescale nowrap v5 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) NO_ANIMATION = 1 myencoding = "utf8" y0 = 1.06865074176486 m = -0.000943519626924529 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 FIT_ERROR = 0 FIT_COV_y0_y0 = 2.20240170189915e-06 FIT_COV_m_y0 = -4.67890413749725e-08 FIT_COV_y0_m = -4.67890413749725e-08 FIT_COV_m_m = 9.9466831564506e-10 ml = -0.0005 mh = -0.0005 dens_Tc = 1.02 Tc = 45.0 g = 1.0 b = 0.01002 ## Last datafile plotted: "lcdemo.dat" plot 'lcdemo.dat' using 1:2:5 with errorbars, density(x) ## Last fit command: "fit l(x) 'lcdemo.dat' using 1:2:5 yerr via y0, m"