clc;
clear all;
% Given Parameters
lambda = 1550*10^-9; % in meters
n_2 = 1.444; % refractive index of the cladding
a = 8*10^-6; % Core radius (given core diameter = 16 micro meters)
n_eff = [1.451190 1.453455 1.454277 1.456903 1.458948]; % Effective refractive index for LP_01 mode
n_1 = 1.474; % refractive index of the core
%f=[12.7 10.8 8.3 7.5 5.3 3.7 1.9]; % delta n
% Calculation
k = (2*pi)/lambda;
%V=[2:2:10];
%V = linspace(0,10,5);
V=linspace(0,10,5);
% b =((n_eff).^2 -(n_2)^2)./((V *lambda)./(2*pi*a)).^2;
b = (n_eff-n_2)/(n_1-n_2);
normalized_b = b./max(b);
%b = f*1e-3./(n_1-n_2);
plot (V,normalized_b);
xlabel('frequence normalise V')
ylabel('constante de propagation b')
%xlim([0 12])
%ylim([0 1])
