
function integrand=integrand(D,D0,lambda,n,np,xsi)

x=pi*D/lambda;
ro=2*x*(n-1);
beta=atan(np/(n-1));

Qc=2-4*exp(-ro.*tan(beta)).*(cos(beta).*sin(ro-beta)./ro+(cos(beta)./ro).^2.*cos(ro-2*beta))+4*(cos(beta)./ro).^2.*cos(2*beta);

integrand=Qc.*pi.*D.^2/4.*(D/D0).^-xsi;