set term postscript eps enhanced color set out 'hysteresis.eps' set xlabel 'Applied field (kA/m)' set ylabel 'M / Ms' versor_x = 1/sqrt(2) versor_y = 1/sqrt(2) versor_z = 0.0 scalar_prod(x1,x2,x3) = x1*versor_x + x2*versor_y + x3*versor_z set mxtics 5 # minor tics and grid set ytics 1 set mytics 5 set grid xtics ytics mxtics mytics lt -1 lw 0.5, lt 0 plot [-1050:1050] [-1.2:1.2] \ 'plot.dat' u (scalar_prod($1,$2,$3)/1000):(scalar_prod($4,$5,$6)) t 'Stoner-Wohlfarth' w lp 4