clear all
close all
clf
% Morisonanvendelse
h_lambda=0:0.01:0.5;
H_D=2*pi* (sinh(2*pi.*h_lambda).^2)./(sinh(4*pi.*h_lambda)./4+pi.*h_lambda);

a=(H_D(3)-H_D(2))/(h_lambda(3)-h_lambda(2));


figure(1)
plot(h_lambda,H_D,...
    [0 h_lambda(end)], [4*pi 4*pi],'--',...
    h_lambda,a*h_lambda,'-.')
ht=title('Fig IV.5 HdynJorde');
hx=xlabel('h / \lambda');
hy=ylabel('H / D');
set([hx hy gca],'fontsize',14)
set(ht,'fontsize',12)