While different methods, very involved from computational point of view, as Hartree-Fock [1]},[2]}, RPA [3]}, Density Functional Theory [4]},[5]},etc. have been used to study the optical spectrum in fullerene cluster, we intend in this paper to provide a study from the perspective of a semi-classical method, namely Time Dependent Thomas Fermi (TDTF). The method, or other variations of the same approximation, has been used before in connection with nuclei[6]}, atoms [7]}, clusters [8]}, strong laser interaction in clusters [9]}, thin metal films [10]} and even transport processes in semiconductors [11]}. Therefore, it is eligible to consider such a method, which, to the best of our knowledge, has never been simulated before in its fully nonlinear form on \(C_{60}\)  fullerene.
