I would lie to display the coefficients before the expansion terms of the cos one by one.
syms x y
f = expand(sin(x + 3*y));
coeffs_cosy = collect(f, cos(y))
How can the terms before the cos^n be displayed one by one?
This is short but for more complicated expressions, several dimensions the taylor series output text that cannot be displayed.