We present a first-principles lattice QCD calculation of the local form factors describing the decay. We focus on the region of large di-muon invariant masses , where the contributions from the four-quarks operators in the effective weak Hamiltonian (which are neglected at present) are expected to be small. We use our results for the form factors to determine the branching fraction for , which has been recently measured by the LHCb in the region .