  In this paper, it is proved that the higher dimensional Hardy operator is
bounded from Hardy space to Lebesgue space. The endpoint estimate for the
commutator generated by Hardy operator and (central) BMO function is also
discussed.
