%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                            Michael Pokojovy                             % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[X, Y] = meshgrid(-2:0.1:2);

figure(1);

set(gcf, 'PaperUnits', 'centimeters');
xSize = 14; ySize = 10;
xLeft = (21 - xSize)/2; yTop = (30 - ySize)/2;
set(gcf,'PaperPosition', [xLeft yTop xSize ySize]);
set(gcf,'Position',[0 0 xSize*50 ySize*50]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
hold on;
t = title('Das Vektorfeld $x \mapsto (-2x_{1} + x_{2}, 3x_{2})^{T}$ im Rechteck $(0, 3) \times (0, 2)$', 'interpreter', 'latex');
set(t, 'FontSize', 16);

M = 15;
N = 15;

[X, Y] = meshgrid(linspace(0, 3, M), linspace(0, 2, N));
U = -2*X + Y;
V =  3*Y;

axis([-0.2 3.2 -0.2 2.2]);

grid on;
colormap gray;

quiver(X, Y, U, V);

plot([0 3 3 0 0], [0 0 2 2 0], 'r');

xlabel('x_{1}');
ylabel('x_{2}');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%