We investigate theB0→J/ψηK0andB+→J/ψηK+decay by using the Dalitz plot analysis. As we know there are tree, penguin, emission, and emission-annihilation diagrams for these decay modes in the factorization approach. The transition matrix element is factorized into aB→ηKform factor multiplied byJ/ψdecay constant and also aB→Kform factor multiplied byJ/ψηdecay constant. According to QCD factorization approach and using the Dalitz plot analysis, we calculate the branching ratios of theB0→J/ψηK0andB+→J/ψηK+three-body decay in view of theη-η'mixing and obtain the value of the(9.22-1.47+2.67)×10-5, while the experimental results of them are(8±4)×10-5and(10.8±3.3)×10-5, respectively. In this research we also analyze theB(Bs)→ηcπK*decay which is similar to the previous decay, but there is no experimental data for the last decay. Since for calculations of theB(Bs)→ηcπK*decay we use assumptions of theB→J/ψηKdecay, we hope that if this decay will be measured by the LHCb in the future, the experimental results will be in agreement with our calculations.