No BSD License
function f=diffAAW(x) t=0:0.01:0.2; y=feval(@eqAAW,t); z=diff(y)./diff(t); f=interp1(t(1:end-1),z,x); end
Contact us at files@mathworks.com