The bubble model by Keller and Prosperetti is adapted to solve the nonlinear oscillation of a gas bubble. This formulation leads to accurate results since it introduces the energy equation instead of the polytropic assumption for the bubble interior. The numerical method used in this study is stable enough to handle large amplitude of bubble oscillation. The numerical results show some interesting nonlinear phenomena fur the bubble oscillator. The excitation changes the natural frequency of the bubble and makes some harmonic resonances at $f/f_0=1/2, 1/3$ and so on. The natural frequency of a bubble oscillator decreases compared with the linear case result, which means that the nonlinear bubble oscillation system is a "softening"system. In addition, the frequency response curve jumps up or down at a certain frequency. It is also found that there exist multi-valued regions in the frequency response curve depending on the initial conditions of bubble. The dependency of the bubble motion on the initial condition can generate extremely large pressure and temperature which might be the cause of the acoustic cavitation and the sonoluminescence.inescence.