h=1/10;
x=0:h:10;
 y=-1:h:1;
 [xx,yy]=meshgrid(x,y);
 %Tensiones normales en las direcciones de los ejes coordenados 
 tni=xx./40;
 tnj=(3.*xx)./40;
 tnk=xx./40;
  %Representación tensiones normales
  figure
 subplot(1,3,1)
 surf(xx,yy,tni)
 axis image
 view(2)
 colorbar
 xlabel('x')
 ylabel('y')
 zlabel('z')
 title('Tensión normal dirección eje i')
 subplot(1,3,2)
 surf(xx,yy,tnj)
 axis image
 view(2)
 colorbar
 xlabel('x')
 ylabel('y')
 zlabel('z')
 title('Tensión normal dirección eje j')
 subplot(1,3,3)
 surf(xx,yy,tnk)
 axis image
 view(2)
 colorbar
 xlabel('x')
 ylabel('y')
 zlabel('z')
 title('Tensión normal dirección eje k')
