clc
close all
clear all
axis([-20 20 -20 20])
axis square
axis off
set(gcf,'Color','white')
[x,y]=meshgrid(-19:19);
x=x;
y=y;
[m,n]=size(x);
A=[19 1;0 1];
x1=zeros(m,n);
y1=zeros(m,n);
for i=1:m
    for j=1:n
        temp=A*[x(i,j);y(i,j)];
        x1(i,j)=temp(1);
        y1(i,j)=temp(2);
    end
end
line(x1,y1,...
    'linestyle','none',...
    'marker','s',...
    'markersize',7,...
    'color','b',...
    'markerfacecolor','b')



B=[19 9;0 1];
x2=zeros(m,n);
y2=zeros(m,n);
for i=1:m
    for j=1:n
        temp=B*[x(i,j);y(i,j)];
        x2(i,j)=temp(1);
        y2(i,j)=temp(2);
    end
end
line(x2,y2,...
    'linestyle','none',...
    'marker','s',...
    'markersize',7,...
    'color','b',...
    'markerfacecolor','b')



C=[0 1;-1 0];
x3=zeros(m,n);
y3=zeros(m,n);
for i=1:m
    for j=1:n
        temp=C*A*[x(i,j);y(i,j)];
        x3(i,j)=temp(1);
        y3(i,j)=temp(2);
    end
end
line(x3,y3,...
    'linestyle','none',...
    'marker','s',...
    'markersize',7,...
    'color','b',...
    'markerfacecolor','b')
x4=zeros(m,n);
y4=zeros(m,n);
for i=1:m
    for j=1:n
        temp=C*B*[x(i,j);y(i,j)];
        x4(i,j)=temp(1);
        y4(i,j)=temp(2);
    end
end
line(x4,y4,...
    'linestyle','none',...
    'marker','s',...
    'markersize',7,...
    'color','b',...
    'markerfacecolor','b')



D=[0 1;1 0];
x5=zeros(m,n);
y5=zeros(m,n);
for i=1:m
    for j=1:n
        temp=D*A*[x(i,j);y(i,j)];
        x5(i,j)=temp(1);
        y5(i,j)=temp(2);
    end
end
line(x5,y5,...
    'linestyle','none',...
    'marker','s',...
    'markersize',7,...
    'color','b',...
    'markerfacecolor','b')
x6=zeros(m,n);
y6=zeros(m,n);
for i=1:m
    for j=1:n
        temp=D*B*[x(i,j);y(i,j)];
        x6(i,j)=temp(1);
        y6(i,j)=temp(2);
    end
end
line(x6,y6,...
    'linestyle','none',...
    'marker','s',...
    'markersize',7,...
    'color','b',...
    'markerfacecolor','b')



E=[1 0;0 -1];
x7=zeros(m,n);
y7=zeros(m,n);
for i=1:m
    for j=1:n
        temp=E*A*[x(i,j);y(i,j)];
        x7(i,j)=temp(1);
        y7(i,j)=temp(2);
    end
end
line(x7,y7,...
    'linestyle','none',...
    'marker','s',...
    'markersize',7,...
    'color','b',...
    'markerfacecolor','b')
x8=zeros(m,n);
y8=zeros(m,n);
for i=1:m
    for j=1:n
        temp=E*B*[x(i,j);y(i,j)];
        x8(i,j)=temp(1);
        y8(i,j)=temp(2);
    end
end
line(x8,y8,...
    'linestyle','none',...
    'marker','s',...
    'markersize',7,...
    'color','b',...
    'markerfacecolor','b')