python - 如何用python解决具有有限差分的非线性ode二阶？

Question

我目前正在尝试用有限差分和python解决这个非线性颂歌：d/dy((du/dy)^n)=a

带有 n 和 a 常数（实数）。这是一个 bvp，在 y=0 时 u=0，在 y=1 时 u=1

我知道如何找到解析解：它是泊肃叶流，幂律液体限制在 2 个无限板之间。

我找不到用 python 和有限差分来处理这种非线性 ode 的方法。有人能帮助我吗？