MathematicaのFindFitで非線形関数を任意のモデル関数にあてはめその係数を求めたい

Question

前提・実現したいこと

MathematicaのFindFitを使って非線形関数を任意のモデル関数にあてはめその係数を求めたいです。

発生している問題・エラーメッセージ

Power: Indeterminate expression 0.^0. encountered.
FindFit: The Jacobian is not a matrix of real numbers at {a,b,c}={1., 1., 1.,}.

該当のソースコード

data = {{0,125},{0.018,124.92},{0.022,124.89},{0.026,124.84},{0.032,124.77},{0.038,124.68},{0.049,124.48},{0.062,124.15},{0.085,123.42},{0.099,122.85},{0.113,122.21},{0.143,120.51},{0.197,116.53},{0.229,113.47},{0.262,110.07},{0.301,105.46},{0.351,98.84},{0.400,91.82},{0.453,83.69},{0.510,74.66},{0.562,66.46},{0.616,58.12},{0.653,52.63},{0.708,44.70},{0.763,37.35},{0.815,30.75},{0.860,25.42},{0.912,19.69},{0.961,14.76},{1.010,10.28},{1.049,7.38},{1.095,5.04},{1.139,3.57},{1.184,2.21},{1.159,1.31}}

model = 1 + a / (1 + (x/b)^c)
fit = FindFit[data, model, {a,b,c},x]

{a→1., b→1., c→1.}

試したこと

FindFitの例題をWOLFRAMの公式サイトより調べ実行しました。
model関数の方に問題があるように感じましたが、問題を特定することができませんでした。