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

figure(1);

set(gcf, 'PaperUnits', 'centimeters');
xSize = 20; 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;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(1, 2, 1);
hold on;
t = title(['Plot der Funktion ', '$f_{1}(x_{1}, x_{2}) = \frac{x_{1} x_{2}^{3}}{(x_{1}^{2} + x_{2}^{2})^{2}}$'], 'interpreter', 'latex');
set(t, 'FontSize', 16);
Z = X.*Y.^3./((X.^2 + Y.^2).^2);
mesh(X, Y, Z);
view([125, 56]);

xlabel('x_{1}');
ylabel('x_{2}');
zlabel('f_{1}(x_{1}, x_{2})');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(1, 2, 2);
hold on;
t = title(['Plot der Funktion ', '$f_{2}(x_{1}, x_{2}) = \frac{x_{1}^{3} x_{2}^{2}}{(x_{1}^{2} + x_{2}^{2})^{2}}$'], 'interpreter', 'latex');
set(t, 'FontSize', 16);
Z = X.^3.*Y.^2./((X.^2 + Y.^2).^2);
mesh(X, Y, Z);
view([125, 56]);

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