Drude模型计算介电常数matlab

clear
clc
lamda=1:0.001:5;
c=310^8;
q=1.60218910^(-19);%单个电子的电荷量
ee=5.5;%背景折射率
me1=9.1095310^(-31);%单个电子的质量
me=0.21me1;
n=3.710^26;%载流子浓度
u=35910^(4);%载流子迁移率
e0=8.8510^(-12);%真空的介电常数
gama=q./(ume);%阻尼
wp=((nq^2)./(me*e0))(0.5);%等离子体频率
w=(2pic)./(lamda10^(-6));
jiedian=ee-((wp)^2)./(w.(gama1i+w));
plot(lamda,real(jiedian))
figure
plot(lamda,imag(jiedian))