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

k = 9.93610952421767:0.00000000000000005:9.93610952421770;
tension=0; %initialize to prevent confusion
tension2=0; %initialize to prevent confusion
Nk = numel(k);
for ii = 1:Nk
   tension(ii) = pwdm(k(ii),0,12);
   tension2(ii) = pwdm2(k(ii),0,12);
endfor
figure(2)
hold off
plot((k-9.9361095242)*1e12,log10(tension),'r')
hold on
plot((k-9.9361095242)*1e12,log10(tension2),'b')
