%% H1Doku10R   Harj.1 teht. 10
% Piirretaan 'Strangin talo' ja sen kuvia 2x2-matriisien antamissa lineaarikuvauksissa. 

% Vrt. Gilbert Strang: Linear Algebra (kansikuva)
% Kts. my?s MIT Open courseware: Lectures by G. Strang

% Ohje: Aja skripti osissa, joita rajaavat %%-merkit:

%       Siirra osoitin ao. osaan ja CTR-ENTER  (CMD-ENTER Mac:ssa)

% Loupuksi voit ajaa koko skriptin painamalla keskella olevaa vihreaa
% nuolta.

% Voit modifioida skripti? esim. poistamalla hold on-komennon ja
% komentamalla figure, jolloin avautuu uusi grafiikkaikkuna, tms.

close all
T=[0 0 -1 6 13 12 12 3 3 6 6 0;0 9 8 15 8 9 0 0 5 5 0 0];
plot(T(1,:),T(2,:))
shg
hold on
%% A_1
A = [1 0 ; 0 2];
S = A*T;
plot(S(1,:),S(2,:),'r')
shg
%% A _2
A = [1 0 ; 0 -1];
S = A*T;
plot(S(1,:),S(2,:),'g')
shg
%% A_3
A = [0 1; 1 0];
S = A*T;
plot(S(1,:),S(2,:),'c')
shg
%% A_4

A = [1 3 ; 0 1];
S = A*T;
plot(S(1,:),S(2,:),'m')

shg

%% Vaihtoehto: Jaetaan grafiikkaruutu osiin.

figure    % Uusi grafiikkaikkuna

subplot(2,3,1)
plot(T(1,:),T(2,:))
ylim([0 30])  % Sama skaala kuin seur. kuvassa, jotta vertailu helpompi.

subplot(2,3,2)
A = [1 0 ; 0 2];
S = A*T;
plot(S(1,:),S(2,:))

subplot(2,3,3)
A = [1 0 ; 0 -1];
S = A*T;
plot(S(1,:),S(2,:))

subplot(2,3,4)
A = [0 1; 1 0];
S = A*T;
plot(S(1,:),S(2,:),'c')

subplot(2,3,5)
A = [1 3 ; 0 1];
S = A*T;
plot(S(1,:),S(2,:),'m')