  We present a fully automated procedure providing an easy way to perform,
systematically, phenomenological analyses in flavor physics for general BSM
scenarios. This procedure relies on MARTY-1.4, is model independent and
requires as input only the Lagrangian of the theory. Once the Lagrangian has
been defined, tree-level and one-loop Wilson coefficients can be calculated
symbolically by MARTY, from which flavor observables can be computed
numerically by available software programs. We focus in particular on $b\to
s\gamma$ and the recently updated $b\to s\ell^+\ell^-$ observables which are in
tension with the SM, and present a general procedure to extract the relevant
one-loop coefficients, such as $C_7$, $C_7^\prime$, $C_9$ and $C_{10}$.
