x = 100:400 ;
y = 600:1500 ;
[X, Y] = meshgrid(x, y);
Z = 3 * X + 2 * Y;
subplot(1, 2, 1);
surf(X, Y, Z, 'EdgeColor', 'none');
fontSize = 15;
xlabel('X', 'FontSize', fontSize);
ylabel('Y', 'FontSize', fontSize);
zlabel('Z', 'FontSize', fontSize);
title('Z = 3 * X + 2 * Y', 'FontSize', fontSize);
colormap(jet(256));
colorbar;
P = X ./ Z;
subplot(1, 2, 2);
imshow(P, [], 'XData', x, 'YData', y);
xlabel('X', 'FontSize', fontSize);
ylabel('Y', 'FontSize', fontSize);
zlabel('Z', 'FontSize', fontSize);
axis('on', 'xy');
title('P = X ./ Z', 'FontSize', fontSize);
colorbar;
impixelinfo;
g = gcf;
g.WindowState = 'maximized'
