syms x a b n f(x) = nthroot(a*x-b,n); F(x) = int(f, x); disp(F);