 t=linspace(-1,1,20);
 x=t;
 y=2*cosh(t/2);
 % Vectores tangentes unitarios interiores
  t1i=(1)./(sqrt(1+(sinh(t/2)).^2));
  t2i=sinh(t/2)./(sqrt(1+(sinh(t/2)).^2));
  % Vectores tangentes unitarios 
  hold on
  plot(x,y,'LineWidth',2);
  quiver(x,y,t1i,t2i);
 hold off
 % Centrado de la gráfica en el origen de coordenadas
 ax = gca;
 ax.XAxisLocation = 'origin';
 ax.YAxisLocation = 'origin';
 % Etiquetas
 title('Vector tangente unitario')
 legend("Catenaria","Vector tangente unitario")
 axis("equal")
 ax.XAxisLocation = 'origin';
 ax.YAxisLocation = 'origin';
 box on
 grid minor