f_euler_derivative, _ = dmodel.f_jacobian(x, u)