function sanityCheck(ka)
kd = pi;
xi = 0:0.05:2*pi;
s = sin(xi);
c = cos(xi);
zeta = -ka*s+kd;
cz = cosh(zeta);
sz = sinh(zeta);
cd = cosh(kd);
sd = sinh(kd);
LHS_D = cz.*s/cd + 0.5*ka*(sz.^2+s.^2)/(sd*cd);
RHS_D = s;
LHS_K = sz.*c/sd - ka*c.*cz.*s/sd;
RHS_K = c;
plot(xi,LHS_D,'r',xi,RHS_D,':r')
hold on
plot(xi,LHS_K,'b',xi,RHS_K,':b')
hold off
leg = {'LHS D','RHS D','LHS K','RHS K'};
legend(leg,'location','southwest')
texText('$x/\lambda$','x')
texText(sprintf('$ka = %0.3f$',ka),'t')
filename = sprintf('comp%03d',1000*ka);
graphToPDF(filename,4,3)
end