3 files changed, 174 insertions, 0 deletions
diff --git a/math2373/lab7/lab7.txt b/math2373/lab7/lab7.txt
new file mode 100644
index 0000000..446d49a
--- /dev/null
+++ b/math2373/lab7/lab7.txt
@@ -0,0 +1,135 @@
+p100OG; p2
+
+p2 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+d=200; p100OG; p2
+
+p2 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+d=500; p100OG; p2
+
+p2 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+d=100; p100OG; p3
+
+p3 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+d=100; p100OG; p4
+
+p4 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+d=100; p100OG; p5
+
+p5 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+p100OG; p5
+
+p5 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+clear p5
+p100OG; p5
+
+p5 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+p100OG; p5
+
+p5 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+p100OG; p5
+
+p5 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+p100OG; p5
+
+p5 =
+
+    0.2264
+    0.1132
+    0.0377
+    0.3396
+    0.2830
+
+diary off
+I=eye(5);R=rref(AI);
+{Undefined function or variable 'AI'.
+} 
+I=eye(5);R=rref(A-I);
+{Matrix dimensions must agree.
+} 
+I=eye(6);R=rref(A-I);
+v=-R(:,5);
+v=-R(:,6);
+v=-R(:,6); v(6)=1;
+p=v/sum(v);
+p
+
+p =
+
+    0.0851
+    0.1702
+    0.1064
+    0.1560
+    0.2695
+    0.2128
+
+diary off
diff --git a/math2373/lab7/p100.m b/math2373/lab7/p100.m
new file mode 100644
index 0000000..3559149
--- /dev/null
+++ b/math2373/lab7/p100.m
@@ -0,0 +1,20 @@
+A=[ 0   1/2 0   0 0   0;
+    0   0   1/3 0 1/2 0;
+    1/4 1/2 0   0 0   0;
+    1/4 0   0   0 1/2 0;
+    1/4 0   1/3 0 0   1;
+    1/4 0   1/3 1 0   0; ];
+%Matrix A from lab
+p1=[1;0;0;0;0;0];
+p2=[0;1;0;0;0;0];
+p3=[0;0;1;0;0;0];
+p4=[0;0;0;1;0;0];
+p5=[0;0;0;0;1;0];
+p6=[0;0;0;0;0;1];
+for n=1:100
+     p1=A*p1;
+     p2=A*p2;
+     p3=A*p3;
+     p4=A*p4;
+     p5=A*p5;
+end
diff --git a/math2373/lab7/p100OG.m b/math2373/lab7/p100OG.m
new file mode 100644
index 0000000..87fcc0f
--- /dev/null
+++ b/math2373/lab7/p100OG.m
@@ -0,0 +1,19 @@
+A=[0   1/3 1/2 1/2 0;
+   1/2 0   0   0   0;
+   0   1/3 0   0   0;
+   0   1/3 1/2 0   1;
+   1/2 0   0   1/2 0];
+%Matrix A from lab
+p1=[1;0;0;0;0];
+p2=[0;1;0;0;0];
+p3=[0;0;1;0;0];
+p4=[0;0;0;1;0];
+p5=[0;0;0;0;1];
+for n=1:600
+     p1=A*p1;
+     p2=A*p2;
+     p3=A*p3;
+     p4=A*p4;
+     p5=A*p5;
+end
+