%----------------------------------------------
%   輸入參數in為波長（單位：m）
%   輸出out為對應的電導率
%----------------------------------------------
function out=getcon(in)
nn=in*1e6;
coon=(-6.21198035452752*nn.^5+ 36.9470450716619*nn.^4-83.1702487377279*nn.^3+...
   94.5079035255191*nn.^2-48.61825472792*nn+9.64235361750187);
eee=(-108.909148353719*nn.^5+ 671.044293885632*nn.^4-1547.0900731642*nn.^3+ ...
   1674.14670244375*nn.^2-773.494216929895*nn+130.284211701206);

e0=(1e-9)/(36*pi);
w=2*pi*3e8./(in);     % 光頻率

out=(coon+i*eee)*e0.*w;