We will review some recent results concerning the numerical approximation of SDEs with singular drift driven by an additive rotationally invariant alpha-stable process. We will focus on the error on the densities for drifts which can belong to Hölder, Lebesgue or Besov spaces. Some open related problems will be discussed as well.