#!/usr/bin/env perl

my $Nz = 4;

sub getE
{
  my $dir = shift;
  my $sum = 0;
  for ($i = 1; $i <= $Nz; $i++) {
    $sum += `gettotalenergy $dir/$i/log`;
  }
  my $E = $sum/$Nz;
  return $E;
}

my $E1 = getE('../02_td_kondo');
my $E0 = getE('../02_td_0');

my $Eimp = $E1-$E0;

print "Impurity binding energy: $Eimp\n";

