format longR=27.9749e-3;fun=@(r,alpha)r.*(R-r.*cos(alpha)./(R^2+R^2-2*R*r.*cos(alpha)).^(3/2));a=quad2d(fun,16.129e-3,16.51e-3,0,2*pi)format short
