How to Define a Complete Bell Polynomial in Mathematica?

  • Context: Mathematica 
  • Thread starter Thread starter anthony2005
  • Start date Start date
  • Tags Tags
    Mathematica Variable
Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
2 replies · 3K views
anthony2005
Messages
24
Reaction score
0
Hello,
how do you define a function or make a list with n elements, where n is any? More precisely there is a function in mathematica, BellY

[itex]BellY[n,k,\{x_{1},...x_{n-k+1}\}][/itex]

which gives the partial Bell polynomial. I would like to define in mathematica the complete Bell polynomial defined as

[itex]CBellY[n,\{x_{1},...x_{n}\}]=\sum_{k=1}^{n}BellY[n,k,\{x_{1},...x_{n-k+1}\}][/itex]

How can I do that?
Thank you.
 
Physics news on Phys.org
Temporarily I have Mathematica forget it knows BellY so you can see the details.

In[1]:= CBellY[n_,v_]:=Sum[BellY[n,k,Take[v,n-k+1]],{k,1,n}]

In[2]:= CBellY[3,{1,2,3}]

Out[2]= BellY[3,1,{1,2,3}]+BellY[3,2,{1,2}]+BellY[3,3,{1}]
 
Great, thank you very much!