%function f = sismo_Rayleigh_disp(k,sigma,g,alphas,betas)
function f = sismo_Rayleigh_disp(k,sigma,alphas,betas)
% This function is zero for Rayleigh waves (without oceans)

R=0; %g/(2*betas);
S=0; %g/(2*alphas);


m=R+sqrt(k.^2-sigma.^2./alphas.^2);
n=S+sqrt(k.^2-sigma.^2./betas.^2);

f=(4.*betas.^4.*k.^2.*(m.*n)-(sigma.^2-2.*k.^2.*betas.^2).^2);
end

