fn_0 = 'f(f(x))'; max_n = 10; %fn = simplify (diff (fn_0)); for n = [6 : max_n] fn(n) = simplify (diff (fn(n - 1))); n end