ConstStr k=t^b*a*0.473/(2.081-1.608*t),
         dk=((((b*(t^(b-1)))*a)*0.473)*(2.081-1.608*t)-(((t^b)*a)*0.473)*(-1.608))/((2.081-1.608*t)^2);
Parameters a,b;
Variable t,r;
ODEFunction r''=-r'*(1+10*dk/k)/10;
Data;
0.228398899        8.32
0.242900416        13.47
0.255672483        15.64
0.267082228        17.69
0.277390305        18.64
0.286788667        18.96
0.295422884        19.35
0.303406043        20.17
0.310827786        21.69