quad8/gd/teoria/02_problema_1_8/gir.m

function anggir = gir(C)
  % Arestes i les seves longituds per la corba
  arestes = C(:, 2:end) - C(:, 1:end-1);
  longs = sqrt(sum(arestes.*arestes));

  % Angles
  prods = sum(arestes(:, 1:end-1).*arestes(:, 2:end));
  cosanggir = prods./(longs(1:end-1).*longs(2:end));
  anggir=acos(cosanggir);