%求解Lame波對稱部分的基本模態So
format long
syms c;   %定義速度c為符號函數
kh=5      %k為波數，h為總的板厚
cd=5790;cs=3200; %cd為縱波波速，cs為橫波波速
y1=4*tanh(kh/2*(1-(c/cd)^2)^(1/2))*(1-(c/cd)^2)^(1/2)*(1-(c/cs)^2)^(1/2);
y2=(2-(c/cs)^2)^2*tanh(kh/2*(1-(c/cs)^2)^(1/2));
ezplot(y1,[0,3500]); %作圖求解
hold on
grid on
ezplot(y2,[0,3500]);
[xx,yy]=ginput(1);
[x(1),fval,exitflag]=fzero(inline(y1-y2),xx(1),[]);
c=x(1)
y1=4*tanh(kh/2*(1-(c/cd)^2)^(1/2))*(1-(c/cd)^2)^(1/2)*(1-(c/cs)^2)^(1/2);
y2=(2-(c/cs)^2)^2*tanh(kh/2*(1-(c/cs)^2)^(1/2));
e1=y1-y2