function d = RD_2nd_deriv_poly(x,y)
%

[p,s] = polyfit(x,y,3);
d = 6*p(1)*x + 2*p(2)*ones(1,length(x));
