  We present the possibility of calculating the quark flavor changing neutral
current decays $B_{s}^0\to \ell \bar \ell$ and $B_{d}^0\to \ell \bar \ell$ for
a large variety of supersymmetric models. For this purpose, the complete
one-loop calculation has been implemented in a generic form in the Mathematica
package SARAH. This information is used by SARAH to generate Fortran source
code for SPheno for a numerical evaluation of these processes in a given model.
We comment also on the possibility to use this setup for non-supersymmetric
models.
