The coupled system, where one is a degenerate parabolic equation and the other has not a diffusion term arises in the modeling of European options with liquidity shocks. Two implicit-explicit (IMEX) schemes that preserve the positivity of the differential problem solution are constructed and analyzed. Numerical experiments confirm the theoretical results and illustrate the high accuracy and efficiency of the schemes in combination with Richardson extrapolation