function makeGraph3()
d = [10,30,100,300];
g     = 9.8;
N     = 100;
lam0  = 1;
lamN  = 1000;
lam   = lam0*exp((log(lamN/lam0)/N)*(0:N))';
k     = 2*pi./lam;
omega = NaN(N+1,length(d));

for n = 1:length(d)
    omega(:,n) = sqrt(g*k.*tanh(k*d(n)));
end
T     = 2*pi./omega;
T_inf = 2*pi./sqrt(g*k);

plot(lam,T_inf,'k')
hold on
plot(lam,T)
hold off

texText('$\lambda$ [m]','x')
texText('$T$ [s]','y')
ylim([0,60])
texText('$d=10$m','g',400,55)
texText('$d=30$m','g',600,45)
texText('$d=100$m','g',730,34)
texText('$d=300$m','g',800,26)
texText('$d=\infty$','g',850,21)
graphToPDF('period.pdf',5,4)
end
