// Scilab script (matrix function)
NN=4;
NN1=2*NN+1;
X=kron([-NN:NN],ones(1,NN1)');
Y=X';
// draw domain
for k=1:NN1
plot(X(k,:),Y(k,:));
plot(X(:,k),Y(:,k));
end
// draw square
a1=[0 0 -1 -1];
b1=[0 1 1 0];
xfpoly(a1,b1);
square(-NN,-NN,NN,NN);
//
// transformation matrix
A=[1 0;1 2];
// transform
for k=1:NN1
for m=1:NN1
  z(1)=X(k,m);
  z(2)=Y(k,m);
  w=A*z;
  X(k,m)=w(1);
  Y(k,m)=w(2);  
end
end
// draw range
set("current_figure",2);
for k=1:NN1
plot(X(k,:),Y(k,:),'r');
plot(X(:,k),Y(:,k),'r');
end
square(-NN,-NN,NN,NN);
//
// draw pararallogram
a2=[0 1 -1 -2];
b2=[0 2 1 -1];
xfpoly(a2,b2);
//