Summary: The Steklov eigenvalue problem of Helmholtz equation is considered. The Steklov eigenvalue problem is reduced to a new variational formula on the boundary of a given domain, in which the selfadjoint property of the original differential operator is kept and the calculating of hyper-singular integral is avoided. A numerical example showing the efficiency of this method and an optimal error estimate are given.