function dy=difftripts(x,y)
% 插值型三点公式求给定点上的数值导数，函数解析式未知
% 输入的x必须是一个等距序列

%  NUMERICAL METHODS: Matlab Programs
% (c) 2007 by Xie Liling
%  Complementary Software to accompany the textbook:
%  Information and Computing Science: A Laboratory Course

h=x(2)-x(1);
n=length(x);
dy=zeros(1,n);
% 左端点利用前向公式求导数
dy(1)=(-3*y(1)+4*y(2)-y(3))/(2*h);
% 右端点利用后向公式求导数
dy(n)=(y(n-2)-4*y(n-1)+3*y(n))/(2*h);
% 其他点利用中心公式求导数
for i=2:n-1
    dy(i)=(-y(i-1)+y(i+1))/(2*h);
end