Using a new numerical EOS (equation of state) table calculated by Shen et al., we performed numerical simulations of protoneutron star cooling. The EOS is based on the relativistic mean field theory, and the parameters in its Lagrangian have been chosen to reproduce the experimental properties of both stable and unstable nuclei. Furthermore, the numerical table covers such a wide range of thermodynamical quantities (temperature, 0 ∼ 100 MeV; electron fraction, 0 ∼ 0.56; density, 10 5.1 ∼ 10 15.4 g/cc) that cooling simulations even for 50 seconds could be done without troubles. The quasistatic evolution of protoneutron stars was investigated in detail with a numerical code including neutrino transfer (MGFLD scheme). Dependencies both on EOS and on initial models and implications to the SN1987A constraint on neutrino oscillation models are discussed.