0001
0002 eps_infty = 1.54;
0003 lambda_p = 177.5;
0004 mu_p = 14500.0;
0005 A1=1.27;
0006 lambda1=470.0;
0007 mu_p1=1900.0;
0008 A2=1.1;
0009 lambda2=325.0;
0010 mu_p2=1060.0;
0011 phi=-pi/4;
0012
0013 epsAu = eps_infty * (1 - 1 ./ (lambda_p.^2 *( (1./lambda).^2 + i./(mu_p.*lambda))))...
0014 + A1 / lambda1 *(exp(i*phi)./(1/lambda1-1./lambda-i/mu_p1)+exp(-i*phi)./(1/lambda1+1./lambda+i/mu_p1))...
0015 + A2 / lambda2 *(exp(i*phi)./(1/lambda2-1./lambda-i/mu_p2)+exp(-i*phi)./(1/lambda2+1./lambda+i/mu_p2));
0016