A brief review is presented for the pre-history and discovery of fullerenes. Single-site potentials with parameters proposed by Girifalco was used to describe the interactions of the fullerene molecules C70 and C96. We present theoretical model for calculating the thermodynamic properties of liquid for both C70 and C96 by means of an improved equation of state (EOS), in which the particles are interacting via pair wise interaction composed of suitable linear combination of three Yukawa functions (3YK). The proposed equation of state provides a powerful mathematical formalism for the Helmholtz free energy and the pressure within the series mean-spherical approximation (SMSA) which are the basic ingredients to compute the liquid-vapour coexistence curve of C70 and C96 as well as the other thermodynamic properties for the bulk liquid and the vapour phases. The comparisons with Gibbs ensemble Monte Carlo (GEMC) simulations and the self-consistent Ornstein-Zernike approximation (SCOZA) were carried out. The estimated critical parameters for both C70 and C96 are TC=2176K, ρC=0.44nm-3, PC=51.64 bars and TC=2477K, ρC=0.32nm-3, PC=44.28 bars respectively. It is to be noted that the obtained results of the thermodynamic properties along the binodal curves of C70 and C 96 are exhibit interesting features.