How to determine the coefficients of a polynom

split_s(input::Num, v::Num) = 
  vcat([Symbolics.coeff(input, v^i) for i in Symbolics.degree(input, v):-1:1],
  substitute(input, v => 0))

gives:

julia> split_s(input,s)
  4U*a
 2 + Q
   J*Q
1 Like