% demolinkuv.m   30.9.03
% Shear

A=[1 3;0 1]

neliox=[0 2 2 0 0] 
nelioy=[0 0 2 2 0]

plot(neliox,nelioy)

nelio=[neliox;nelioy]

kuva=A*nelio

plot(kuva(1,:),kuva(2,:))

     %% Parannellaan:

fill(neliox,nelioy,'b')
     axis([-0.5 2.5 -0.5 2.5])
shg

     figure(2)
fill(kuva(1,:),kuva(2,:),'y')
tol=0.5
     a=min(kuva(1,:))-tol;b=max(kuva(1,:))+tol
     c=min(kuva(2,:))-tol;d=max(kuva(2,:))+tol
     axis([a b c d])
shg
