clc; clear; close all;
f = @(x,y) x*y*sin(x)*cos(y)

h=0.1
xvalues = 0:h:7;
yvalues = 0:h:5;
for i=1:length(xvalues)
    for j=1:length(yvalues)
        fvalues(i,j) = f(xvalues(i),yvalues(j));
    end
end
%mesh(fvalues)
mesh(yvalues, xvalues, fvalues)

volum = 0;
for i = 1:length(xvalues)
    areal = integrer(yvalues, fvalues(i,:));
    volum = volum + areal * h;
end
title(volum)