MSSM Higgs production cross section in bottom-quark annihilation is evaluated at next-to-next-to-leading order (NNLO) in QCD. Scale dependence for both the factorization and renormalization scales for center of mass energies 7, 8 TeV was found. With the help of computer program HIGLU the neutral MSSM Higgs production cross section via gluon fusion at LHC with center of mass energies 7, 8, 14 TeV including next-to-next-to-leading order QCD corrections is presented. The result for neutral Higgs boson mass is in accordance with the last experimental data received at the LHC.