clc
clear
format long
KJ=2.3;%1.4:0.2:3;
lamda=0.5:0.1:3.5;
f=0.1:0.1:0.6;
f_eff=0.625;
n=length(f);
N_si=[];
N_he=[];
N_cl=[];
for i=1:n
    d=KJ.*f(i);
    N0=SiN(lamda);
    N1=HE11_EffN(lamda,KJ,d,f_eff);
    N2=CladN(lamda,KJ,d)
    N_si=[N_si,N0];
    N_he=[N_he,N1];
    N_cl=[N_cl,N2];
    i
end
save HE_f.mat N_si N_he N_cl KJ lamda f f_eff