x=1/2+(0:1); % abscisas de los puntos medios [X,Y,Z]=meshgrid(x); % malla de puntos medios f=X+Y-Z; % valores de la función sobre la malla sum(f(:)) % suma de Riemann (el volumen de cada subcaja es 1)Esta suma es una aproximación de la integral triple de $f(x,y,z)$ sobre la caja. El valor de esa integral también puede obtenerse en el ordenador, escribiendo
syms x y z % se declaran estas variables simbólicas int(int(int(x+y-z,z,0,2),y,0,2),x,0,2) % cálculo de la integral
for n=2:10 % ciclo para dibujar las secciones hold on fill3([0 0 1 1], (1/n)*ones(1,4),[0 1 1 0],'g') % cada sección es un cuadrado alpha(.2) % transparencia end axis([0 1 0 1 0 1]) % tamaño de los ejes box on % incluye caja xlabel('x');ylabel('y');zlabel('z'); % etiquetas de los ejes hold off
[x,y,z] = meshgrid(0:.1:2); v = -x.*y.*z; xslice = 2; yslice = (0:10)/5; zslice = 0; slice(x,y,z,v,xslice,yslice,zslice) colormap jet shading interp xlabel('x');ylabel('y');zlabel('z'); box on view([-120,32])
figure(1) subplot(1,3,1) fill3([1 0 0],[0 2 0],[0 0 0],'g', 'EdgeColor','m') axis equal xlabel('x');ylabel('y');zlabel('z'); view([141,32]) axis([0 1.2 0 2.2 0 3.2]) subplot(1,3,2) fill3([1 0 0],[0 0 0],[0 0 3],'r', 'EdgeColor','m') axis equal xlabel('x');ylabel('y');zlabel('z'); view([141,32]) axis([0 1.2 0 2.2 0 3.2]) subplot(1,3,3) fill3([0 0 0],[0 2 0],[0 0 3],'y', 'EdgeColor','m') axis equal xlabel('x');ylabel('y');zlabel('z'); view([141,32]) axis([0 1.2 0 2.2 0 3.2]) figure(2) fill3([1 0 0],[0 2 0],[0 0 3],'b', 'EdgeColor','m') axis equal xlabel('x');ylabel('y');zlabel('z'); alpha(.6) view([141,32])