normalized propagation constant vs v anyone please can u help me to correct the program

clc;

clear all; % Given Parameters

lambda = 1550 * 10^-9; % in meter

n_2 = 1.444; % refractive index of the cladding

a = 8 * 10^-6 ; % Core radius % given core diameter = 16 micro meter

% n_eff = 1.451190;

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);

%b = f*1e-3./(n_1-n_2);

plot (V,b);

xlabel('frequence normalise V')

ylabel('constante de propagation b')

%xlim([0 12])

%ylim([0 1])

Answer 1

[Jack Goodhew]

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])