clc
clear
format long
KJ=0.4;
lamda=0.5:0.1:2.5;
f=0.79;
d=KJ.*f;
f_eff=1;
N0=SiN(lamda);
N1=HE11_EffN(lamda,KJ,d,f_eff);
N2=CladN(lamda,KJ,d)
save HE_0.7.mat N0 N1 N2 KJ lamda f f_eff
plot(lamda,N0,'-b',lamda,N1,':',lamda,N2,'-b')