Tag: runge-kutta

def second(y, z, y0, z0, t0, tf, dt): N = (tf-t0)/dt y.append(y0) for i in range(int(N)): t1 = t0 + dt y1 = y0 + z0*dt zk1 = z(t0, y0, z0) zk2 = z(t0+1/2*dt, y0+1/2*dt, z0+1/2*zk1*…