Starting from a simple model introduced by L. Pietronero et al., we have employed a programming scheme which can greatly simplify the calculation, and have completed the whole numerical evaluation. A method for the exact determination of τs and the specific calculation are given. Extending the results obtained we further discuss the variation of τs due to surface relaxation.