Experiment!
d = 1e-3; % plate thickness [m]
h = d/2;
ct = 2.26 * 1e-3 / 1e-6; % transverse wave speed for aluminium[m/s]
omega = linspace(10e6,15e6,10);
k = (0:100)*40;
cla
hold on
for i = 1:length(omega)
w = sin(sqrt((omega(i).^2/ct.^2)-k.^2).*h);
plot(k,w,'color',rand(1,3))
text(k(1),w(1),num2str(i))
end
hold off
