function d = RD_2nd_deriv_quad(y,h)
%

num_pts = length(y);
d = zeros(1,num_pts);
for p = 2:num_pts-1
    d(p) = (y(p-1) - 2*y(p) + y(p+1))/h^2;
end
