diff --git a/algorithms/arithmetic_analysis/fixed_point_iter.m b/algorithms/arithmetic_analysis/fixed_point_iter.m
new file mode 100644
index 0000000..1a36423
--- /dev/null
+++ b/algorithms/arithmetic_analysis/fixed_point_iter.m
@@ -0,0 +1,29 @@
+function fpi(f,a, b, x0, N)
+  
+% generate the cobweb diagram associated with
+% the iteration: x_{n+1}=f(x_n).
+% N is the number of iterates, and
+% (a,b) is the interval
+% x0 is  the initial point.
+% plot y=f(x)
+
+x=a:.01:b;
+y=f(x);
+figure
+hold on;
+title('fixed point iteration')
+plot(x,y,'k'); % plot the function
+plot(x,x,'r'); % plot straight line
+x(1)=x0; % plot orbit starting at x0
+for i=1:N
+    x(i+1)=f(x(i));
+    if i==1
+        line([x(i),x(i)],[0,x(i+1)]);
+    else 
+        line([x(i),x(i)],[x(i),x(i+1)]);
+    end
+    pause
+    line([x(i),x(i+1)],[x(i+1),x(i+1)]);
+    pause
+end
+hold off;