常微分方程式 python フィッティング

pythonでdy/dx=-ky^2の式でフィッティングしたいのですがどうやってやればいいですか。
一次式の場合は、下記のコードで実行できたのですが常微分方程式のdy/dx=-ky^2でフィッティングしたいです。
よろしくお願いします。
x = np.array([20, 30, 40, 50, 60])
y = np.array([30, 38, 41, 49, 62])

p = np.polyfit(x, y, 1)

f = np.poly1d(p)

plt.scatter(x, y)
plt.plot(x, f(x))