In supersymmetry with large {\tan \beta} the decays B⁰(\bar {B}⁰) \rightarrow l⁺l^- are dominated by the scalar and pseudoscalar Higgs penguin diagrams leading to strong enhancement of leptonic decay rates with potentially large CP asymmetries in the \tau ⁺\tau ^- decay modes measurable in BELLE or BaBar experiments. The TAUOLA  \tau -lepton decay library supplemented by its universal interface can efficiently be used to search for B⁰(\bar {B}⁰)\rightarrow \tau ⁺\tau ^- decays, and to investigate how the CP asymmetry is reflected in realistic experimental observables.

PACS numbers: 13.25.Hw, 14.60.Fg, 11.30.Er