In this work, we present a Mathematica package Peng4BSM@LO which calculates the contributions to the Wilson Coe±cients of certain e®ective operators originating from the one-loop penguin Feynman diagrams. Both vector and scalar external legs are considered. The key feature of our package is the ability to ¯nd the corresponding expressions in almost any New Physics model which extends the SM and has no °avor changing neutral current (FCNC) transitions at the tree level.