function dfunc = fdiff(func,n) if(n>1) for I = 1:n, syms x; symf = feval(deriv,x); deriv = diff(symf,x); end end func = inline(char(deriv));