% beta2 Shear wave velocity of layer (i+1) [m/s] % rho1 Mass density of layer i [kg/m^3] % rho2 Mass density of layer (i+1) [kg/m^3] p1 = C_beta*X(2)+s*S_beta*X(3 ...