Probably, you can do `Add` for `p` and `q` separately, then multiply to `sin(x**2)` or `cos(x**2)`.  No need to use subs or \_eval_nseries, right?