summaryrefslogtreecommitdiff
path: root/62/CH7/EX7.27/ex_7_27.sce
diff options
context:
space:
mode:
Diffstat (limited to '62/CH7/EX7.27/ex_7_27.sce')
-rwxr-xr-x62/CH7/EX7.27/ex_7_27.sce11
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