def forward_euler(equation: ExplicitNavierStokesODE, time_step: float) -> TimeStepFn:
    return jax.named_call(navier_stokes_rk(ButcherTableau(a=[], b=[1]), equation, time_step), name='forward_euler')
def forward_euler(equation: ExplicitNavierStokesODE, time_step: float) -> TimeStepFn:
    return jax.named_call(navier_stokes_rk(ButcherTableau(a=[], b=[1]), equation, time_step), name='forward_euler')