Vasicek with jumps at a switching intensity

Upward jumps in the intensity that arrive more often in one regime.

The intensity follows $dx = \kappa(\theta_y - x)\,dt + \sigma_y\,dW + dJ$, where $J$ jumps at rate $\ell_y$ by exponential amounts of mean $m$. Jumps change only the time function:

$$g_i(t) = -\kappa\theta_i B + \tfrac12\sigma_i^2 B^2 + \ell_i\Big(\frac{1}{1 + m B} - 1\Big), \qquad B = \frac{1 - e^{-\kappa t}}{\kappa}.$$

Here $\kappa = 2$, $\theta = (0.05, 0.02)$, $\sigma = (0.02, 0.01)$, $\ell = (3, 0.2)$ and $m = 0.03$. The table gives the error in $a_1(3)$ after each order.

switching rateorder 0order 1order 2order 3order 4order 5order 6
52.8e-34.4e-65.0e-72.6e-72.2e-81.1e-81.1e-9
101.4e-31.1e-67.8e-81.7e-88.7e-101.8e-101.0e-11
207.1e-42.6e-71.1e-81.1e-93.0e-112.9e-128.9e-14
403.5e-46.3e-81.4e-96.9e-119.8e-134.7e-147.8e-16