  A new accurate method for reconstructing the arrival direction of an
extensive air shower (EAS) is described. Compared to existing methods, it is
not subject to minimization of a function and, therefore, is fast and stable.
This method also does not need to know detailed curvature or thickness
structure of an EAS. It can have angular resolution of about 1 degree for a
typical surface array in central regions. Also, it has better angular
resolution than other methods in the marginal area of arrays.
