global a b rho f w mu lambda vs vp ks kp kp2 ks2 b2 epsi
a=2000;                                    %//.............[m]        rayon Omega_a
b=4000;                                    %//.............[m]        rayon Omega_b
rho=1;                                  %//.............[kg/m3]    masse volumique
f=4;                                       %//.............[Hz]       frequence
w=2*pi*f;                                  %//.............[rad/s]    pulsation
mu=4e6;                                    %//.............[Pa]       deuxieme coefficient de lame & shear modulus
lambda = 8e6;                              %//.............[Pa]       premier coefficient de lame
vs=sqrt(mu/rho);                           %//.............[m/s]      celerite de l'onde s
vp=sqrt((lambda+2*mu)/rho);                %//.............[m/s]      celerite de l'onde p
ks=w/vs;                                   %//.............[1/m]      nombre d'onde
kp=w/vp;                                   %//.............[1/m]      nombre d'onde
epsi =1e-6;                                %//.............[]         Condition d'arret
kp2=kp*kp;
ks2=ks*ks;
b2=b*b;
 
