function [Xxx,Xzz]=hex_polarisabilities(V,eps,w)

Lz=0.5*((1-3/w)/(1+3/w)+1);
Lx=0.25*((1-0.5*w^0.9)/(1+0.5*w^0.9)+1);
Xxx=V/(4*pi)*(eps-1)/(Lx*(eps-1)+1);
Xzz=V/(4*pi)*(eps-1)/(Lz*(eps-1)+1);
