% Try PWDM for a circular stadium.
% Plot results of PWDM as a function of k. Illustrate difference between pwdm and pwdm2.

k = 0.2:0.01:6;
L = 1.0
tension=0; %initialize to prevent confusion
Nk = numel(k);
for ii = 1:Nk
   tension(ii) = pwdm2(k(ii),L,12);
endfor
figure(1)
hold off
plot(k,log10(tension),'r')
