global a b rho_a rho_b omega mu_a mu_b lambda_a lambda_b vs vp ks_a ks_b kp_a kp_b
a=2000;                                    %//.............[m]        rayon Omega_a
b=8000;                                    %//.............[m]        rayon Omega_b
f=4;                                       %//.............[Hz]       frequence
omega=2*pi*f;                                  %//.............[rad/s]    pulsation

% domaine_b couronne
rho_b=1;                                  %//.............[kg/m3]    masse volumique
mu_b=4e6;                                    %//.............[Pa]       deuxieme coefficient de lame & shear modulus
lambda_b = 8e6;                              %//.............[Pa]       premier coefficient de lame
vs=sqrt(mu_b/rho_b);                           %//.............[m/s]      celerite de l'onde s
vp=sqrt((lambda_b+2*mu_b)/rho_b);                %//.............[m/s]      celerite de l'onde p
ks_b=omega/vs;                                   %//.............[1/m]      nombre d'onde
kp_b=omega/vp;                                   %//.............[1/m]      nombre d'onde                              

% domaine_a disk
rho_a=2;                                  %//.............[kg/m3]    masse volumique
mu_a=32e6;                                    %//.............[Pa]       deuxieme coefficient de lame & shear modulus
lambda_a = 64e6;                              %//.............[Pa]       premier coefficient de lame
vs_a=sqrt(mu_a/rho_a);                           %//.............[m/s]      celerite de l'onde s
vp_a=sqrt((lambda_a+2*mu_a)/rho_a);                %//.............[m/s]      celerite de l'onde p
ks_a=omega/vs_a;                                   %//.............[1/m]      nombre d'onde
kp_a=omega/vp_a; 
