function [v,u]=HartreeEnergy()

[v,u]=Constant.(UnitSystem()).HartreeEnergy;

end