Based on the fully relativistic multiconfiguration Dirac-Hartree-Fock(MCDHF)method and the corresponding pro-gram package GRASP2018,a new program for calculating the polarizabilities is developed.As the first application,the static electric-dipole polarizabilities of the ground state 2s21 S0 and excited state 2s2p 3P0 of beryllium are calculated.By means of these polarizabilities,the blackbody radiation(BBR)shift of the 2s2p 3P0→2s2 1S0 clock transition is deter-mined.The present results agree very well with other available theoretical results.