clear;
clf;
STB = [4 2 1 0 8 1 1 1];
A = STB'*STB;
B = [0 1 0 0 0 0 0 0;
     0 0 1 0 0 0 0 0;
     0 0 0 1 0 0 0 0;
     0 0 0 0 1 0 0 0;
     0 0 0 0 0 1 0 0;
     0 0 0 0 0 0 1 0;
     0 0 0 0 0 0 0 1;
     4 2 1 0 8 1 1 1];
 k = 0;
 I = 1;
 for I = 1:1000
     C = A + (10*sin(k*pi/100))*B;
     lambdaC = eig(C);
        for J = 1:8
            x(J,I)=lambdaC(J);
        end
     k = k + 10;
 end
 plot(x,'o'), grid on, xlabel('Sumbu Nyata'),ylabel('Sumbu Khayal')