function yd=deriv3p(x,y)
% Derivada de y respecto de x, usando fórmula de 3 puntos

n=length(x);

h = x(2)-x(1);

if n<3
  fprintf("La fórmula para 3 puntos es...para 3 puntos\n");
  return
endif

yd(1) = (-3*y(1)+4*y(2)-y(3))/(2*h);
yd(n) = (3*y(n)-4*y(n-1)+y(n-2))/(2*h);
for k=2:n-1
  yd(k) = (y(k+1)-y(k-1))/(2*h);
endfor

endfunction
