在numpy中解释欧拉方法的代码

Question

我是一个蟒蛇新手，试图学习用数字来解决差异化的问题。我是从这门课程中学到的：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]

如果有人可以解释这是什么我会很感激。达到这一点的一切都是有道理的，我知道这应该只是欧拉的方法，但我只是不明白它是如何工作的。