diff options
Diffstat (limited to '62/CH7/EX7.27/ex_7_27.sce')
-rwxr-xr-x | 62/CH7/EX7.27/ex_7_27.sce | 11 |
1 files changed, 11 insertions, 0 deletions
diff --git a/62/CH7/EX7.27/ex_7_27.sce b/62/CH7/EX7.27/ex_7_27.sce new file mode 100755 index 000000000..92878c268 --- /dev/null +++ b/62/CH7/EX7.27/ex_7_27.sce @@ -0,0 +1,11 @@ +clc;
+syms n;
+A=[2 0 0;0 -2 1;0 4 1];
+eig=spec(A);
+m=poly(eig(1:2),"A",["roots"]);
+disp(m,"minimal polynomial is")
+e=spec(A);
+b1=(e(2)^n-e(1)^n)/(e(2)-e(1));
+b0=e(2)^n-b1*e(2);
+An=b0*eye(A)+b1*A;
+disp(An,"A^n")
\ No newline at end of file |