function [Z, ModuleZ, Phi] = rlc(R,L,C, Freq)
%--------------------------------------------------------------------------------
% calcul de l'impédence effective Z du rircuit RLC en parallèle, du module
% de Z et de son argument (angle de phase). Evitez les boucles for, while, ...
%*** paramètres:
% R, L, C       -in-    paramètres des pièces R [Ohm], L [H], C [F]
% Freq(n,1)     -in-    fréquences du courant alterné en [Hz]
% Z(n,1)        -out-   impédence effective Z(i) pour chaque fréquence f=Freq(i)
% ModuleZ(n,1)  -out-   module de Z
% Phi(n,1)      -out-   argument (angle de phase) de Z en degrés, -180<=Z(i)<=180
%--------------------------------------------------------------------------------

%*** nombre des lignes dans le vecteur Freq:
  n = size(Freq,1);

%*** init aléatoire (distribution normale N(0,1)) de la partie réelle et imaginaire:
  Z = randn(n,1) + j*randn(n,1);
%*** init aléatoire (distrubution uniforme dans l'intervalle ouvert (0,1)):
  ModuleZ = rand(n,1);
%*** init à zéro de l'argument de Z (=Phi):
  Phi = zeros(n,1);

%=== A REMPLIR ICI: ...
w = 2*pi*Freq;
wrl = w.*(R*L);
sub = w.*L + j.*((w.^2).*(L*C)-R);
Z = wrl./sub

ModuleZ = abs(Z)

Phi = angle(Z).*(180/pi)

%=== FIN DE SOLUTION
  return;       % n'est pas obligatoire ici, la fonction peut se terminer avec le end final
end