#!/bin/bash
gnuplot --persist <<EOF
set termoption enh
set title "Single impurity Anderson model - renormalization flow diagram"
set xlabel 'NRG step number N (even)'
set ylabel 'Energy/{/Symbol w}_N'
set yrange [0:3.3]
plot 'e-S.dat' w lp title 'Q=0, S=0', \
 'e-D.dat' w lp title 'Q=1, S=1/2', \
 'e-D2.dat' w lp title 'Q=-1, S=1/2', \
 'e-T.dat' w lp title 'Q=0, S=1'
EOF