function y = area(x,epsilon,alpha,k)
y = 4*pi*x.^2./k.^2.*epsilon^4.*alpha.^5./(epsilon.^2-1)./x.^5.*sqrt((epsilon.^2-1)./(alpha^2*epsilon.^2-x.^2));
