#!/bin/bash
gnuplot --persist <<EOF
set termoption enh
set title "SIAM"
set logscale x
set format x '10^{%L}'
set xlabel 'Temperature'
set xrange [1e-8:1e-1]
set ylabel 'C(T)/k_B'
plot 'td-C.dat' with lp title 'impurity heat capacity'
EOF
