#!/bin/sh
gnuplot --persist <<EOF
set termoption enh
set title "Two-channel Kondo model - renormalization flow diagram"
set xlabel 'NRG step number N (even)'
set yrange [0:3.2]
set ylabel 'Energy/{/Symbol w}_N'
plot 'e-D1.dat' w lp title 'Q=0, S=1/2', \
 'e-D2.dat' w lp title 'Q=2, S=1/2', \
 'e-D3.dat' w lp title 'Q=-2, S=1/2', \
 'e-S1.dat' w lp title 'Q=1, S=0', \
 'e-S2.dat' w lp title 'Q=-1, S=0', \
 'e-T1.dat' w lp title 'Q=1, S=1', \
 'e-T2.dat' w lp title 'Q=-1, S=1'
EOF
