diff options
Diffstat (limited to '50/CH6/EX6.3/ex_6_3.sce')
-rwxr-xr-x | 50/CH6/EX6.3/ex_6_3.sce | 17 |
1 files changed, 17 insertions, 0 deletions
diff --git a/50/CH6/EX6.3/ex_6_3.sce b/50/CH6/EX6.3/ex_6_3.sce new file mode 100755 index 000000000..ea77769e4 --- /dev/null +++ b/50/CH6/EX6.3/ex_6_3.sce @@ -0,0 +1,17 @@ +// example 6.3
+// solution to the given IVP
+
+disp('du/dt= A*u');
+// u=[u1 u2]';
+A=[-3 4 ;-2 3]; // given
+B=[1 0;0 1]; // identity matrix;
+
+
+
+
+[x,lam] = geigenvectors(A,B);
+
+// hence;
+disp('u=c1*%e^t*x(:,1)+c2*%e^-t*x(:,2)');
+disp('u1=c1*%e^t+c2*%e^-t*2')
+disp('u2=c1*%e^t+c2*%e^-t')
\ No newline at end of file |