Use this smooth two-dimensional rotational velocity field. Compare Euler and fixed-step Heun at equal NFE, then inspect an adaptive Euler–Heun pair that uses their discrepancy to select steps.