  We investigate the contributions of neutral Higgs bosons to nonleptonic
transition $b \to q \bar{s} s, q=d, s$ under the supersymmetric context. Their
effects to decay width and CP violation in corresponding exclusive decays are
explored. The anomalous dimension matrices of the operators which have to be
incorporated to include the contributions of neutral Higgs bosons are given. We
find that when tan$\beta$ is large (say, 50) and neutral Higgs bosons are not
too heavy (say, 100 GeV), contributions of neutral Higgs penguin can dominate
electroweak penguin contributions, and for some processes, they can greatly
modify both decay width and CP asymmetry.
