1 files changed, 24 insertions, 0 deletions
diff --git a/examples/NRIplot.m b/examples/NRIplot.m
new file mode 100644
index 0000000..3d77fc1
--- /dev/null
+++ b/examples/NRIplot.m
@@ -0,0 +1,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;