set colors classic
  set view equal xyz
  set ticslevel 0
  set palette model HSV functions 5./6.*gray, 1, 1

  set xlabel 'x'
  set ylabel 'y'
  set zlabel 'z'
  set cblabel 'rho/rhomean'
  splot \
    'density.txt' u 6:7:8:5 w p pt 5 lc palette

  pause -1 'Press enter to continue: '