set term post enh color solid
  set output 'meandensity.ps'
  set encoding iso_8859_1
  set colors classic

# set term pngcairo
# set output 'meandensity.png'
# set encoding utf8
# set minussign
# set colors classic

# kernel = 'cubic'
# kernel = 'quintic'
  kernel = 'wendc2'
# kernel = 'wendc2b'
# kernel = 'wendc6'
# kernel = 'wendc8'

  set title kernel.' kernel'
  set log x
  set xrange [9:400]
# set yrange [0.98:1.04]
# set yrange [0:5]
  set yrange [0.96:1.1]
  set ylabel 'rho_{sim} / rho_{ini}'
  set xlabel 'number of neighbors'
  plot \
    1 w l lt 7 lc rgb '#dddddd' notitle, \
    kernel.'-glass_10x10x10'  u 1:5 w l lt 5 notitle, \
    kernel.'-grid_10x10x10'   u 1:5 w l lt 1 notitle, \
    kernel.'-packed_12x12x12' u 1:5 w l lt 3 notitle, \
    kernel.'-glass_10x10x10'  u 1:5:(sqrt($4)/$2) w errorbars lt 5 pt 7 t 'glass', \
    kernel.'-grid_10x10x10'   u 1:5:(sqrt($4)/$2) w errorbars lt 1 pt 7 t 'cubic grid', \
    kernel.'-packed_12x12x12' u 1:5:(sqrt($4)/$2) w errorbars lt 3 pt 7 t 'close packing'

  ics = 'packed_12x12x12'
# ics = 'grid_10x10x10'
# ics = 'glass_10x10x10'

  set title ics noenh
  plot \
    1 w l lt 7 lc rgb '#dddddd' notitle, \
    'cubic-'  .ics u 1:5 w l lt 1 notitle, \
    'quintic-'.ics u 1:5 w l lt 3 notitle, \
    'wendc2-' .ics u 1:5 w l lt 4 notitle, \
    'wendc6-' .ics u 1:5 w l lt 5 notitle, \
    'wendc8-' .ics u 1:5 w l lt 7 notitle, \
    'cubic-'  .ics u 1:5:(sqrt($4)/$2) w errorbars lt 1 pt 7 t 'cubic kernel', \
    'quintic-'.ics u 1:5:(sqrt($4)/$2) w errorbars lt 3 pt 7 t 'quintic kernel', \
    'wendc2-' .ics u 1:5:(sqrt($4)/$2) w errorbars lt 4 pt 7 t 'wendland c2 kernel', \
    'wendc6-' .ics u 1:5:(sqrt($4)/$2) w errorbars lt 5 pt 7 t 'wendland c6 kernel', \
    'wendc8-' .ics u 1:5:(sqrt($4)/$2) w errorbars lt 7 pt 7 t 'wendland c8 kernel'