package Stats::Fermi::FermiDiracStats;
sub new {
        my ($class) = @_;
        my $self = {
                k => 1.380649 * pow(10,-23), 
        };
        $class = ref($class) || $class;
        bless $self, $class;
}
sub funcall {
        my ($self, $e, $mu, $T) = shift;
        return (1 / (exp(($e-$mu)/($self->{k}*$T)) + 1));
}
1;