Algebraic reduction of Feynman diagrams to scalar integrals: A Mathematica implementation of LERG-I