#!/bin/sh
python <<EOF
# -*- coding: utf-8 -*-

import matplotlib as mpl
mpl.use('agg')
import matplotlib.pyplot as plt
import numpy as np

f = plt.figure()

plt.title("Kondo model", fontsize=12)
plt.xlabel(r'\$T/D\$', fontsize=12)
plt.ylabel(r'\$S(T)\ [k_B]\$', fontsize=12)
plt.xlim(1e-8,1e-1)
plt.xscale('log')
plt.ylim(bottom=0)
plt.tick_params(labelsize=12)

plt.tight_layout()

x1,y1 = np.loadtxt('td-S.dat', unpack=True)
l1, = plt.plot(x1,y1, color='black', label='Impurity entropy')

plt.legend()
plt.show()
plt.savefig('entropy.pdf')

EOF