1 files changed, 31 insertions, 0 deletions
diff --git a/62/CH7/EX7.47/ex_7_47.sce b/62/CH7/EX7.47/ex_7_47.sce
new file mode 100755
index 000000000..6aa12c5c0
--- /dev/null
+++ b/62/CH7/EX7.47/ex_7_47.sce
@@ -0,0 +1,31 @@
+clc;
+A=[2 1 0;0 2 1;0 0 2];
+V=[2 0 0;0 2 0; 0 0 2];
+N=[0 1 0;0 0 1;0 0 0];
+z=N*N*N;
+if z==0 then
+    disp("N is a nilpotent of index 3");
+else
+    disp("N is not nilpotent")
+end
+if V*N==N*V then
+    disp("V and N commute")
+else
+    disp("V and N dont commute")
+end
+syms t
+//e^(A*t)=e^(v*t)*e^(N*t)
+e=spec(V);//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;
+eVt=b0*eye(V)+b1*V+b2*V^2;
+disp(eVt,"e^(Vt)")
+e=spec(N);//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;
+eNt=b0*eye(N)+b1*N+b2*N^2;
+disp(eNt,"e^(Nt)")
+eAt=eVt*eNt;
+disp(eAt,"e^(A*t)=")
+\ No newline at end of file