Calculation of molecular polarizabilities for large systems within the random phase approximation