[X,Y]=meshgrid(-1:1); surf(X,Y,2-3*Y) shading interp hold on quiver3(0,0,2,0,3,1) hold off alpha(.7) axis equal xlabel('x');ylabel('y');zlabel('z')
[X,Y]=meshgrid(0:.1:2,-1:1); surf(X,Y,sqrt(4-X.^2)) shading interp hold on quiver3(1,0,sqrt(3),1/sqrt(3),0,1) hold off alpha(.7) axis equal xlabel('x');ylabel('y');zlabel('z')
[X,Y]=meshgrid(-1:.1:1,0:.1:2); surf(X,Y,1+X.^2+Y.^2) shading interp hold on quiver3(0,1,2,0,-2,1) hold off alpha(.7) axis equal xlabel('x');ylabel('y');zlabel('z')
[X,Y]=meshgrid(-.7:.1:.7,-1.7:.1:0-.3); surf(X,Y,sqrt(4-(X.^2+Y.^2))) shading interp hold on quiver3(0,-1,sqrt(3),0,-1/sqrt(3),1) hold off alpha(.7) axis equal xlabel('x');ylabel('y');zlabel('z')
[X,Y]=meshgrid(2:.1:4,-1:.1:1); surf(X,Y,2-sqrt(X.^2+Y.^2)) shading interp hold on quiver3(3,0,-1,1,0,1) hold off alpha(.7) axis equal xlabel('x');ylabel('y');zlabel('z')
[X,Y]=meshgrid(1:3,1:.1:3); surf(X,Y,Y.^2) shading interp hold on quiver3(2,2,4,0,-4,1) hold off alpha(.7) axis equal xlabel('x');ylabel('y');zlabel('z')
Halla el vector normal en los puntos $P_{ij}$, evalúa los cuatro productos $|{\bf N}|\mbox{área}(R_{ij})$ y súmalos. Esa suma es una aproximación (poco precisa, pues sólo tomamos cuatro subrectángulos de $R$) del área de $S$.
Supón ahora que $g(x,y,z)=x^2+z^2$ diera en cada punto el valor de la densidad de masa superficial. Calcula una aproximación de la masa de $S$ sumando los cuatro valores $$g(P_{ij})|{\bf N}|\mbox{área}(R_{ij})$$
[X,Y]=meshgrid(-2:.1:2,-1:.1:1); Z=1-2*X+3*Y; AU=ones(size(X)); subplot(1,2,1) surf(X,Y,Z,'EdgeColor','none') colormap gray shading interp subplot(1,2,2) surf(X,Y,2*(5-4*X+5*Y)) shading interp xlabel('x');ylabel('y');zlabel('z') subplot(1,2,1) hold on [X,Y]=meshgrid(-2:1:2,-1:1:1); AU=ones(size(X)); Z=1-2*X+3*Y; quiver3(X,Y,Z,2*AU,-3*AU,AU); quiver3(X,Y,Z,2*Z,-2*ones(size(X)),-2*Y,'r');%rotacional view([-30,58]) alpha(.3) axis equal xlabel('x');ylabel('y');zlabel('z') hold off