1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
% Illustrates the use of NewtonsRingsI by plotting Newton's rings
thetaoffset=20*pi/180;
rseed=[0:0.001:sqrt(0.1)].^2;
thetaseed=[0:1:360]*pi/180;
[theta,r]=meshgrid(thetaseed,rseed);
[x,y]=pol2cart(theta,r);
subplot(111);
h=pcolor(x,y,NewtonRingsI([0.01,100,14e3,-.1e3,1,1,1.024e-6,thetaoffset],r)); hold on
%h_contours=contour(x,y,sqrt(x.^2+y.^2),[0.02,0.02],'r');
colormap('bone');
set(h,'EdgeColor','none');
set(h,'FaceColor','interp');
set(gca,'Visible','off');
set(gcf,'Color','black');
%figtext(4.3,-0.7,['(',num2str(p),',',num2str(m),')'],8);
axis square
hold off
shg;
|
