  % Definición de los parámetros
  a=-1;
  b=1;
  h=0.09;
  t=a:h:b;
  % Definición de la curva
  x=t;
  y=cosh(t);
  % Vectores normales unitarios orientación interior
 n1i=sinh(t/2)./(sqrt(1+(sinh(t/2)).^2));
 n2i=-(t./t)./(sqrt(1+(sinh(t/2)).^2));
 % Vectores normales unitarios orientación exterior
 n1e=-sinh(t/2)./(sqrt(1+(sinh(t/2)).^2));
 n2e=(t./t)./(sqrt(1+(sinh(t/2)).^2));
 hold on
 plot(x,y,'LineWidth',2);
 quiver(x,y,n1i,n2i);
 quiver(x,y,n1e,n2e);
 hold off
 % Centrado de la gráfica en el origen de coordenadas
 ax = gca;
 ax.XAxisLocation = 'origin';
 ax.YAxisLocation = 'origin';
 % Etiquetas
 title('Vectores normales')
 legend("Catenaria","Vector normal interior","Vector normal exterior")
 axis("equal")
 ax.XAxisLocation = 'origin';
 ax.YAxisLocation = 'origin';
 box on
 grid minor