# set terminal pngcairo transparent enhanced font "arial,10" fontscale 1.0 size 500, 350 # set output 'fit.35.png' set ticslevel 0.2 set xlabel "x" set ylabel "y" set zlabel "z" set zrange [ -3.00000 : 3.00000 ] noreverse nowriteback 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) h(x,y) = sqrt(r*r - (abs(x-x0))**2.2 - (abs(y-y0))**1.8) + z0 phi(x) = (x - phi0)/360.0*2.0*pi main(x) = c11*sin(phi(x))**2 + c33*cos(phi(x))**2 + c44 mixed(x) = sqrt( ((c11-c44)*sin(phi(x))**2 +(c44-c33)*cos(phi(x))**2)**2 +(2.0*(c13+c44)*sin(phi(x))*cos(phi(x)))**2 ) vlong(x) = sqrt(1.0/2.0/rho*1e9*(main(x) + mixed(x))) vtrans(x) = sqrt(1.0/2.0/rho*1e9*(main(x) - mixed(x))) f(tc)= (tanh(Y(tc)) + abs(cos(2.*tb)) * tanh(abs(Y(tc)*cos(2.*tb)))) / (Y(tc)*(1.+(cos(2.*tb))**2)) W(x) = 1./(sqrt(2.*pi)*eta) * exp( -1. * x**2 / (2.*eta**2) ) Y(tc) = tc/sin(tb) * Fhkl * r0liV Q(tc) = (r0*Fhkl/V)**2 * (lambda**3/sin(2.*tb)) * P * f(tc) a(x) = W(x) * Q(tc) / mu R(x) = sinh(A*a(x)) * exp(-1.*A*(1.+a(x))) y0 = 0.000600937793093794 m = -0.000943519626921127 GPFUN_l = "l(x) = y0 + m*x" FIT_CONVERGED = 1 FIT_NDF = 117 FIT_STDFIT = 0.100111507833005 FIT_WSSR = 1.17261073806994 ml = -0.00100003942457045 mh = -0.000831266464088135 dens_Tc = 1.02497044549157 Tc = 46.0899236389968 g = 3.85585054766573 b = 0.00153901241268336 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)" FIT_LIMIT = 1e-10 r = 1.00101991715444 x0 = -0.000309703201613002 z0 = 0.00253574954583606 GPFUN_h = "h(x,y) = sqrt(r*r - (abs(x-x0))**2.2 - (abs(y-y0))**1.8) + z0" FIT_MAXITER = 0 rho = 1000.0 phi0 = -0.161896293627274 GPFUN_phi = "phi(x)\t = (x - phi0)/360.0*2.0*pi" c11 = 5.3401483346306 c33 = 12.4010700704255 c44 = 1 GPFUN_main = "main(x) = c11*sin(phi(x))**2 + c33*cos(phi(x))**2 + c44" c13 = 4 GPFUN_mixed = "mixed(x) = sqrt( ((c11-c44)*sin(phi(x))**2\t\t\t\t +(c44-c33)*cos(phi(x))**2)**2 +(2.0*(c13+c44)*sin(phi(x))*cos(phi(x)))**2 )" GPFUN_vlong = "vlong(x) = sqrt(1.0/2.0/rho*1e9*(main(x) + mixed(x)))" GPFUN_vtrans = "vtrans(x) = sqrt(1.0/2.0/rho*1e9*(main(x) - mixed(x)))" GPFUN_f = "f(tc)= (tanh(Y(tc)) + abs(cos(2.*tb)) * tanh(abs(Y(tc)*cos(2.*tb)))) / (Y(tc)*(1.+(cos(2.*tb))**2))" mu = 0.113046900551349 t0 = 0.18 tb = 0.199278608299778 A = 0.020759275611633 P = 0.924693446208538 Fhkl = 3.42318325539711 r0 = 2.81794092e-13 lambda = 7.09338062818239e-09 V = 1.62253546981499e-23 r0liV = 123.194394853936 eta = 0.000100781677728629 GPFUN_W = "W(x) = 1./(sqrt(2.*pi)*eta) * exp( -1. * x**2 / (2.*eta**2) )" GPFUN_Y = "Y(tc) = tc/sin(tb) * Fhkl * r0liV" GPFUN_Q = "Q(tc) = (r0*Fhkl/V)**2 * (lambda**3/sin(2.*tb)) * P * f(tc)" tc = 0.0020212816909931 GPFUN_a = "a(x) = W(x) * Q(tc) / mu" GPFUN_R = "R(x) = sinh(A*a(x)) * exp(-1.*A*(1.+a(x)))" a0 = 3.11977682286917 a1 = 0.511296458281725 a2 = 0.546666545458654 splot a0/(1+a1*x**2+a2*y**2), 'fit3.dat' in 3 u 1:2:4 ## fit a0/(1+a1*x**2+a2*y**2) 'fit3.dat' in 3 u 1:2:4:(1) via a0,a1,a2