function dfunc = fdiff(func)
syms x;
symf = feval(func,x);
deriv = diff(symf,x);
dfunc = inline(char(deriv));