For the preliminary estimation of the vertical hull natural frequency, the Schlick's or Schlick-type formulae have been traditional ones and are still in common use today. Some investigators have made their efforts, based on statistical data of ships' system parameters, to extend the applicability of Schlick-type formulae to higher modes, or to utilize the Rayleigh method. For instance, the work done by Dinsenbacher et al.[5] belongs to the former and that of Nagamoto et al.[6] to the latter. In a part of his previous paper[7], the author, investigating the case of a cargo ship of medium size, suggested that provided statistically simplified curves such as trapezoid of system parameter distributions are available in hands, direct utlization of an ordinary computer program can be also an another convenient method by which we can obtain both natural frequencies and normal mode shapes. In this paper, to confirm the feasibility of the above suggestion, four high speed boats are investigated. The system parameters of them are originally given in [5]. The computer program used here is one confiled based on a calculation method derived from Myklestal-Prohl modeling of hull, transfer matrix formulation and an extended Gumbel's initial value method for solving frequency equation. The results of the investigation show that the direct calculation based on statistically oriented and reasonably assumed trapezoidal mean curves of system parameter distributions can give us natural frequencies within about 5% deviation up to several-noded modes and normal mode shapes serviceable at least up to 4- or 5-noded modes in comparision with those based on actual distributions of system parameters. For this simplified method the actual data required for input are only of ship length, displacement, total added mass, bending and shear rigidity at amidship. They are available at the early stage of design. By this method we can also easily trace variations of vibration characteristics in the course of ship design cycles.