summaryrefslogtreecommitdiff
path: root/50/CH6/EX6.3/ex_6_3.sce
diff options
context:
space:
mode:
Diffstat (limited to '50/CH6/EX6.3/ex_6_3.sce')
-rwxr-xr-x50/CH6/EX6.3/ex_6_3.sce17
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