syms x;
f(x) = tan(x);
f_x = diff(f,x);
disp(f_x);