sub min { eigenstates( any(@_) <= all(@_) ) }

sub max { eigenstates( any(@_) >= all(@_) ) } 

##</code><code>##

sub max {
  my $max = shift;
  $max = $max > $_ ? $max : $_ for @_;
  return $max
}