function E = E_field_LP_Mode(l,gamma_c,gamma_g, r_c, theta,rho)
%   Cette fonction calcule la distribution du champ électrique d'un mode
%   guidé LP_lm.
%   Variables d'entrée :
%     l : valeur du paramèttre l (L) du mode;
%     gamma_c et gamma_g : valeur de ces paramètres définit en cours pour
%     un mode donné
%     r_c : rayon du coeur
%     theta,rho : coordonnées polaires des points où le champ doit être
%     calculé.
%-------------------------------------------------------------------------
f = zeros(size(rho)); %il est essentiel de prédéfinir les dimensions de f
Rho1 = rho <= r_c; % Dans coeur, fonction Bessel_J
Rho2 = rho > r_c;  % Dans gaine, fonction Bessel_K
g = cos(l*theta);
B = besselj(l,gammma_c*r_c)/besselk(l,gamma_g*r_c);
f(Rho1) = besselj(l, gamma_c*rho(Rho1)); %calcule à l’intérieur du coeur
f(Rho2) = B*besselk(l,gamma_g*rho(Rho2)); %calcule à l’extérieur du cœur (i.e. dans gaine)
E = g.*f;
end