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 rate | order 0 | order 1 | order 2 | order 3 | order 4 | order 5 | order 6 |
|---|---|---|---|---|---|---|---|
| 5 | 2.8e-3 | 4.4e-6 | 5.0e-7 | 2.6e-7 | 2.2e-8 | 1.1e-8 | 1.1e-9 |
| 10 | 1.4e-3 | 1.1e-6 | 7.8e-8 | 1.7e-8 | 8.7e-10 | 1.8e-10 | 1.0e-11 |
| 20 | 7.1e-4 | 2.6e-7 | 1.1e-8 | 1.1e-9 | 3.0e-11 | 2.9e-12 | 8.9e-14 |
| 40 | 3.5e-4 | 6.3e-8 | 1.4e-9 | 6.9e-11 | 9.8e-13 | 4.7e-14 | 7.8e-16 |