8 files changed, 111 insertions, 0 deletions
diff --git a/70/CH5/EX5.1.1/5_1_1.sci b/70/CH5/EX5.1.1/5_1_1.sci
new file mode 100755
index 000000000..997e4a46f
--- /dev/null
+++ b/70/CH5/EX5.1.1/5_1_1.sci
@@ -0,0 +1,12 @@
+//page 238
+clear;
+close;
+clc;
+A=[3 0;0 2];
+eig=spec(A);
+[V,Val]=spec(A);
+disp(eig,'Eigen values:')
+x1=V(:,1);
+x2=V(:,2);
+disp(x1,x2,'Eigen vectors:');
+//end
+\ No newline at end of file
diff --git a/70/CH5/EX5.1.2/5_1_2.sci b/70/CH5/EX5.1.2/5_1_2.sci
new file mode 100755
index 000000000..1660fc20f
--- /dev/null
+++ b/70/CH5/EX5.1.2/5_1_2.sci
@@ -0,0 +1,13 @@
+//page 238
+clear;
+close;
+clc;
+disp('The eigen values of a projection matrix are 1 or 0.')
+P=[1/2 1/2;1/2 1/2];
+eig=spec(P);
+[V,Val]=spec(P);
+disp(eig,'Eigen values:')
+x1=V(:,1);
+x2=V(:,2);
+disp(x1,x2,'Eigen vectors:');
+//end
+\ No newline at end of file
diff --git a/70/CH5/EX5.2.1/5_2_1.sci b/70/CH5/EX5.2.1/5_2_1.sci
new file mode 100755
index 000000000..62388df91
--- /dev/null
+++ b/70/CH5/EX5.2.1/5_2_1.sci
@@ -0,0 +1,11 @@
+//page 238
+clear;
+close;
+clc;
+A=[1/2 1/2;1/2 1/2];
+[V,Val]=spec(A);
+disp(Val,'Eigenvalue matrix:');
+disp(V,'S=');
+disp(A*V,'AS=S*eigenvaluematrix')
+disp('Therefore inv(S)*A*S=eigenvalue matrix')
+//end
+\ No newline at end of file
diff --git a/70/CH5/EX5.2.2/5_2_2.sci b/70/CH5/EX5.2.2/5_2_2.sci
new file mode 100755
index 000000000..b7c21513d
--- /dev/null
+++ b/70/CH5/EX5.2.2/5_2_2.sci
@@ -0,0 +1,15 @@
+//page 238
+clear;
+close;
+clc;
+disp('The eigenvalues themselves are not so clear for a rotation.')
+disp('90 degree rotation')
+K=[0 -1;1 0];
+disp(K,'K=')
+eig=spec(K);
+[V,Val]=spec(K);
+disp(eig,'Eigen values:')
+x1=V(:,1);
+x2=V(:,2);
+disp(x1,x2,'Eigen vectors:');
+//end
+\ No newline at end of file
diff --git a/70/CH5/EX5.2.3/5_2_3.sci b/70/CH5/EX5.2.3/5_2_3.sci
new file mode 100755
index 000000000..86315146b
--- /dev/null
+++ b/70/CH5/EX5.2.3/5_2_3.sci
@@ -0,0 +1,15 @@
+//page 249
+clear;
+close;
+clc;
+disp('K is rotation through 90 degree,then K^2 is rotation through 180 degree and inv(k is rotation through -90 degree)')
+K=[0 -1;1 0];
+disp(K,'K=')
+disp(K*K,'K^2=')
+disp(K*K*K,'K^3=')
+disp(K^4,'K^4=')
+[V,D]=spec(K);
+disp('K^4 is a complete rotation through 360 degree.')
+disp(D,'Eigen value matrix,D of K:');
+disp(D^4,'and also D^4=')
+//end
+\ No newline at end of file
diff --git a/70/CH5/EX5.3.1/5_3_1.sci b/70/CH5/EX5.3.1/5_3_1.sci
new file mode 100755
index 000000000..589c51652
--- /dev/null
+++ b/70/CH5/EX5.3.1/5_3_1.sci
@@ -0,0 +1,23 @@
+//page 249
+clear;
+close;
+clc;
+A=[0 4;0 1/2];
+disp(A,'A=');
+eig=spec(A);
+disp(eig,'Eigen values:')
+[v,D]=spec(A);
+u0=[v(:,1)];//Taking u0 as the 1st eigen Vector.
+for k=0:5
+    disp(k,'k=');
+    u=A*u0;
+    disp(u,'U(k+1)(K from 0 to 5)')
+    u0=u;
+end
+u0=[v(:,2)];//Taking u0 as the 2nd eigen vector.
+for k=0:5
+    disp(k,'k=');
+    u=A*u0;
+    disp(u,'U(k+1)=')
+    u0=u;
+end
+\ No newline at end of file
diff --git a/70/CH5/EX5.5.1/5_5_1.sci b/70/CH5/EX5.5.1/5_5_1.sci
new file mode 100755
index 000000000..75c535fec
--- /dev/null
+++ b/70/CH5/EX5.5.1/5_5_1.sci
@@ -0,0 +1,12 @@
+//page282
+clear;
+close;
+clc;
+i=sqrt(-1);
+x=3+4*i;
+disp(x,'x=');
+x_=conj(x);
+disp(x*x_,'xx_=');
+r=sqrt(x*x_);
+disp(r,'r=')
+//end
diff --git a/70/CH5/EX5.5.2/5_5_2.sci b/70/CH5/EX5.5.2/5_5_2.sci
new file mode 100755
index 000000000..2e0b68631
--- /dev/null
+++ b/70/CH5/EX5.5.2/5_5_2.sci
@@ -0,0 +1,10 @@
+//282
+clear;
+close;
+clc;
+i=sqrt(-1);
+x=[1 i]';
+y=[2+1*i 2-4*i]';
+disp(x'*x,'Length of x squared:');
+disp(y'*y,'Length of y squared:');
+//end
+\ No newline at end of file