#!/bin/gnuplot # Gnuplot 5 script using terminal qt (export is done via this terminal) # Plots number of particles times IPR versus frequency # for systems in dimension 4 at Δϕ/ϕ_J = 0.030 set term qt 0 persist enhanced font 'serif,14' f_1k = "IPR_1k.dat" f_4k = "IPR_4k.dat" f_8k = "IPR_8k.dat" f_16k = "IPR_16k.dat" nb = 4 set format "%g" # Coloring function (shading red to blue) color(i) = 65536*int(255*(i-1)/(nb-1)) + 255-int(255*(i-1)/(nb-1)) set logscale xy set xrange [1e-2:] set yrange [:1000] set key title "N" at 0.08,250 set xlabel "ω" set ylabel "N×Y(ω)" set label 1 sprintf("(a) 4d - Δϕ/ϕ_J = 0.030") at 0.013,600 plot f_1k using 1:($2*1024) lc rgb color(1) w l t "1024", \ f_4k using 1:($2*4096) lc rgb color(2) w l t "4096", \ f_8k using 1:($2*8192) lc rgb color(3) w l t "8192", \ f_16k using 1:($2*16384) lc rgb color(4) w l t "16384"