A variational method for calculating dynamic polarizabilities