Before revision:  To solve this system, we first fix the optimal exercise boundary for each regime resulting in multi-variable fixed domains.
After revision:  To solve this system, we first employ the logarithmic transformation to map the free boundary for each regime resulting in multi-variable fixed domains.