% spcplot.m

load spc20040915.txt
dtr = pi/180;
tha = 0:360;
thar = tha*dtr;
xtha = sin(thar);
ytha = cos(thar);
a = -0.112:0.0035:0.112;
b = 1:65;
b = b./b;
kx = b'*a;
ky = kx';
for (i=0:700)
 clf
 hold on
 axis([-0.141 0.141 -0.112 0.112])
 ib=1+65*i;
 contour(kx,ky,spc20040915(ib:ib+64,1:65),'k')

 for wl=100:100:300
    r = 2*pi/wl;
    plot(r*xtha,r*ytha,'k')
 end
 plot([0 0],0.112*[-1 1],'k',0.112*[-1 1],[0 0],'k')
 pause(0.1)
end
