2 files changed, 152 insertions, 16 deletions
diff --git a/math2373/lab7/lab7.txt b/math2373/lab7/lab7.txt
index e69de29..446d49a 100644
--- a/math2373/lab7/lab7.txt
+++ 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
index 26ce350..3559149 100644
--- a/math2373/lab7/p100.m
+++ b/math2373/lab7/p100.m
@@ -1,19 +1,20 @@
-clear d;
-A=[0 1/3 1/2 1/2 0;
-   1/4 0 0 0 0;
-   1/4 1/3 0 0 0;
-   1/4 1/3 1/2 0 1;
-   1/2 0 0 1/2 0];
+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];
-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:500
-     %p1=A*p1;
-     %p2=A*p2;
+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;
+     p4=A*p4;
+     p5=A*p5;
 end