diff options
Diffstat (limited to '62/CH7/EX7.44/ex_7_44.sce')
-rwxr-xr-x | 62/CH7/EX7.44/ex_7_44.sce | 16 |
1 files changed, 16 insertions, 0 deletions
diff --git a/62/CH7/EX7.44/ex_7_44.sce b/62/CH7/EX7.44/ex_7_44.sce new file mode 100755 index 000000000..e42f88c09 --- /dev/null +++ b/62/CH7/EX7.44/ex_7_44.sce @@ -0,0 +1,16 @@ +clc;
+syms t
+A=[0 -2 1;0 0 3;0 0 0];
+z=A*A*A;
+if z==0 then
+ disp("A is a nilpotent of index 3");
+else
+ disp("A is not nilpotent")
+end
+//we know %e^(at)=b0*I+b1*A+b2*A^2;
+e=spec(A);//eigen values
+b2=t^2*%e^(e(3)*t);
+b1=t*%e^(e(2)*t)-b2*2*e(2);
+b0=%e^(e(1)*t)-b1*e(1)-b2*e(1)^2;
+eAt=b0*eye(A)+b1*A+b2*A^2;
+disp(eAt,"e^(At)")
\ No newline at end of file |