h=1/20;                 %paso de muestreo
u=1/3:h:1;              %Intervalo [1/3,1]
v=-1:h:1;               %Intervalo [-1,1]
[U,V]=meshgrid(u,v);    %Matrices de las coordenadas u y v
figure(1)  
xx=U.*V;                 %Paramerización
yy=1/2*(U.^2-V.^2); 
%Tensiones normales
figure(7)
Eu=(-14*U.^2-6*V.^2)./(U.^2+V.^2); % Campo escalar
Ev=(-10*U.^2-2*V.^2)./(U.^2+V.^2); % Campo escalar
subplot(1,2,1);surf(xx,yy,Eu)          % mallado
axis([-1,1,-1,1])     % región
view(2)                % punto de vista
subplot(1,2,2);surf(xx,yy,Ev)          % mallado2
axis([-1,1,-1,1])      % región
view(2)                % punto de vista