Modelling Cavitation Bubble Radius Oscillation using Equation Based Modelling

Hello, I am trying to model the radius, R(t), oscillation of a cavitation bubble in liquid (compressible) subjected to alternating pressure at ultrasonic frequency. I would like to model the radius oscillation using the Gilmore model equation as shown in the attached .jpeg file.
The link to the published paper is as follows: link.springer.com/content/pdf/10.1007%2Fs11661-012-1188-3

I have tried to model equation (2) through the ODE module and by specifying the derivative of H by d(H,R), and managed to get the model running. However, the result seems to be very different from the results for the change in the bubble radius, R(t), presented in the paper.
The followings are the parameters, variables and also the ODE function I specified for the model based on the values considered in the paper:

Parameters
f 20000[Hz] Frequency of ultrasound
omega 2*pi*f
acoustic_pressure 0.6e6[Pa]

Variables
B 3000[atm]
n 7
rho 970[kg/m^3]
k_exponent 1.6 For air
viscosity 2.66e-3[Pa*s] SCN
sigma 3.85e-2[N/m] SCN
P_static 1[atm]
R_initial 100e-6[m]
P_ext acoustic_pressure*cos(omega*t)
P_infinity P_static+P_ext
Press (2*sigma/R_initial+P_static)*(R_initial/R)^(3*k_exponent)-2*sigma/R-4*Rt*viscosity/R
C sqrt(n/rho*(Press+B)*((Press+B)/(P_static+B))^(-1/n))
H n/(n-1)*((P_static+B)^(1/n))/rho*((Press+B)^(n-1/n)-(P_infinity+B)^(n-1/n))

ODE Function (R)
(1-Rt/C)*R*Rtt+3/2*(Rt)^2*(1-Rt/(3*C))-(1+Rt/C)*H-R*Rt/C*d(H,R)*(1-Rt/C)

Kind regards,

Charles