解释numpy中Euler方法的代码

Question

我是一个python新手，尝试用它来数值求解微分方程。我是从这门课上学的：http://nbviewer.jupyter.org/github/numerical-mooc/numerical-mooc/blob/master/lessons/01_phugoid/01_02_Phugoid_Oscillation.ipynb

我一直在努力理解这些代码的作用：

for n in range(1,N):

    u = u + dt*numpy.array([u[1], g*(1-u[0]/zt)])

    z[n] = u[0]

如果有人能解释一下这是怎么回事，我会很感激的。到目前为止，一切都是有意义的，我知道这应该只是欧拉的方法，但我不明白它是如何工作的。在