syms x; f(x) = x^x; f_x = diff(f,x); f_xx = diff(f_x,x); f_xxx = diff(f_xx,x); disp("f'(x):"); disp(simplify(f_x,100)); disp("f''(x):"); disp(simplify(f_xx,100)); disp("f'''(x):"); disp(simplify(f_xxx,100));