f_euler = dmodel.f(x, u, k=0)