csv = csvread("punktman_n.csv");

stopnie = 90;
radiany = stopnie * pi / 180;
sx = 1;
sy = 0.5;
tx = 1;
ty = -0.5;

MacierzObrot = [cos(radiany) -sin(radiany) 0; sin(radiany) cos(radiany) 0; 0 0 1];
MacierzSkala = [sx 0 0; 0 sy 0; 0 0 1];
MacierzPrzesun = [1 0 tx; 0 1 ty; 0 0 1];
moc = MacierzPrzesun*MacierzSkala*MacierzObrot;

MatrixXY1 = [csv(:,1)';csv(:,2)';ones(1,length(csv(:,1)))];

for i = 1:3:length(MatrixXY1(1,:))
   MatrixXY1(:,i:i+2) = moc * MatrixXY1(:,i:i+2);
endfor

subplot (1, 2, 1)
plot(csv(:,1),csv(:,2),'sk');
subplot (1, 2, 2)
plot(MatrixXY1(1,:),MatrixXY1(2,:),'sk');