diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /62/CH7/EX7.44 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '62/CH7/EX7.44')
| -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 |