Figures 3 and 4 show similar results as Figs. 1 and 2, but for the strength functions in . As for the case, the IS and IV strength functions from the Arnoldi iteration agree very well with the strength functions of Ref. [17]. The convergence of strength functions is as fast as for ; after iterations the smoothed strength functions change only by about . We thus have to make only iterations to calculate reasonably accurate strength functions for the RPA problem whose dimension is . The large double spikes observed in Fig. 4 below 10MeV are due to the lowest RPA phonons, which by the smoothing procedure acquire 100-keV widths and move slightly down in excitation energy.