We present numerical tests and predictions of the KrkNLO method for matching of NLO QCD corrections to hard processes with LO parton shower Monte Carlo generators (NLO+PS). This method was described in detail in our previous publications, where it was also compared with other NLO+PS matching approaches (MC@NLO and POWHEG) as well as fixed-order NLO and NNLO calculations. Here we concentrate on presenting some numerical results (cross sections and distributions) for $Z/gamma^*$ (Drell-Yan) and Higgs-boson production processes at the LHC. The Drell--Yan process is used mainly to validate the KrkNLO implementation in the Herwig 7 program with respect to the previous implementation in Sherpa. We also show predictions for this process with the new, complete, MC-scheme parton distribution functions and compare them with our previously published results. Then, we present the first results of the KrkNLO method for Higgs production in gluon-gluon fusion at the LHC and compare them with MC@NLO and POWHEG predictions from Herwig 7 fixed-order results from HNNLO and a resummed calculation from HqT, as well as with experimental data from the ATLAS collaboration.