RPA calculations with Gaussian expansion method