clear all
load cac40.m;
u=diff(cac40(1:2807));
j=length(u);
i=j-10;
q=[u(1:i+1) u(2:i+2) u(3:i+3) u(4:i+4) u(5:i+5) u(6:i+6) u(7:i+7) u(8:i+8) u(9:i+9) u(10:i+10)];
X=((2806.^(-(1/2)))*q); 
V=X'*X;
[c,a2]=eig(V);
b2=diag(a2);
[Descend,Ind]=sort(b2);
s1=Descend;
b=diag(s1);
k=diag(b);
Q=V*c;
g=b*c;
f=(X*c)'*(X*c);
z=X*c;

subplot(3,3,1);
plot(X);
title('X');
axis([2 2806 -5 5]);
grid;

subplot(3,3,2);
plot(X');
title('(transp(X))' );
axis([1 10 -5 5]);
grid;

subplot(3,3,3);
plot(V);
title('V=((Transp(X)).X)');
axis([1 10 -10 1000]);
grid;

subplot(3,3,4);
plot(c);
title('c: Eigenvectors of V');
axis([1 10 -1 1]);
grid;

subplot(3,3,5);
plot(b);
title('b: Eigenvalues of V');
axis([1 10 0 1000]);
grid;

subplot(3,3,6);
plot(Q);
title('V.c');
axis([1 10 -500 500]);
grid;

subplot(3,3,7);
plot(g);
title('b.c');
axis([1 10 -1000 1000]);
grid;

subplot(3,3,8);
plot(f);
title('(Transp(X.c)).(X.c)');
axis([1 10 0 1000]);
grid;

subplot(3,3,9);
plot(z);
title('X.c');
axis([1 2806 -5 5]);
grid;

disp('Press "Enter" to continue')
pause

subplot(3,3,1);
plot(c(:,1),'b')
;title('C(1)');
axis([1 10 -0.5 0.5]);
grid;

subplot(3,3,2);
plot(c(:,2),'b');
title('C(2)');
axis([1 10 -0.5 0.5]);
grid;

subplot(3,3,3);
plot(c(:,3),'b');
title('C(3)');
axis([1 10 -0.5 0.5]);
grid;

subplot(3,3,4);
plot(c(:,4),'b');
title('C(4)');
axis([1 10 -0.5 0.5]);
grid;

subplot(3,3,5);
plot(c(:,5),'b');
title('C(5)');
axis([1 10 -0.5 0.5]);
grid;

subplot(3,3,6);
plot(c(:,6),'b');
title('C(6)');
axis([1 10 -0.5 0.5]);
grid;

subplot(3,3,7);
plot(c(:,7),'b');
title('C(7)');
axis([1 10 -0.5 0.5]);
grid;

subplot(3,3,8);
plot(c(:,8),'b');
title('C(8)');
axis([1 10 -0.5 0.5]);
grid; 

subplot(3,3,9);
plot(c(:,9),'b');
title('C(9)');
axis([1 10 -0.5 0.5]);
grid;

disp('Press "Enter" to continue')
pause

subplot(3,3,1);
plot(z(:,1),'b');
title('X.C (1)');
axis([2 2806 -5 5]);
grid; axis('off');


subplot(3,3,2);
plot(z(:,2),'b');
title('X.C (2)');
axis([2 2806 -5 5]);
grid;
axis('off')

subplot(3,3,3);
plot(z(:,3),'b');
title('X.C (3)');
axis([2 2806 -5 5]);
grid; axis('off')

subplot(3,3,4);
plot(z(:,4),'b');
title('X.C (4)');
axis([2 2806 -5 5]);
grid; axis('off')

subplot(3,3,5);
plot(z(:,5),'b');
title('X.C (5)');
axis([2 2806 -5 5]);
grid; axis('off')

subplot(3,3,6);
plot(z(:,6),'b');
title('X.C (6)');
axis([2 2806 -5 5]);
grid; axis('off')

subplot(3,3,7);
plot(z(:,7),'b');
title('X.C (7)');
axis([2 2806 -5 5]);
grid; axis('off')

subplot(3,3,8);
plot(z(:,8),'b');
title('X.C (8)');
axis([2 2806 -5 5]);
grid; axis('off')

subplot(3,3,9);
plot(z(:,9),'b');
title('X.C (9)');
axis([2 2806 -5 5]);
grid;; axis('off')

disp('Press "Enter" to continue')
pause

subplot(3,3,1);
plot(z(:,1),z(:,2),'b');
title('X.C (1,2)');
grid;

subplot(3,3,2);
plot(z(:,2),z(:,3),'b');
title('X.C (2,3)');grid;

subplot(3,3,3);
plot(z(:,3),z(:,4),'b');
title('X.C (3,4)');grid;

subplot(3,3,4);
plot(z(:,4),z(:,5),'b');
title('X.C (4,5)');grid;

subplot(3,3,5);
plot(z(:,5),z(:,6),'b');
title('X.C (5,6)');grid;

subplot(3,3,6);
plot(z(:,6),z(:,7),'b');
title('X.C (6,7)');grid;

subplot(3,3,7);
plot(z(:,7),z(:,8),'b');
title('X.C (7,8)');grid;

subplot(3,3,8);
plot(z(:,8),z(:,9),'b');
title('X.C (8,9)');grid;

subplot(3,3,9);
plot(z(:,9),z(:,10),'b');
title('X.C (9,10)');grid;


% "Complex and Chaotic Nonlinear Dynamics. 
% Advances in Economics and Finance, 
% Mathematics and Statistics" 
% T.Vialar, Springer 2009.
% Copyright(c).