#5109: [with patch, positive review] add support for Bell polynomials in Sage
---------------------------+------------------------------------------------
Reporter: mhansen | Owner: mhansen
Type: enhancement | Status: new
Priority: minor | Milestone: sage-3.3
Component: combinatorics | Resolution:
Keywords: |
---------------------------+------------------------------------------------
Comment (by wdj):
One more test (also positive):
{{{
sage: n=6
sage: add([bell_polynomial(n,i)((1,)*(n-i+1)) for i in range(1,n+1)]) ==
bell_number(n)
True
sage: n = 7
sage: add([bell_polynomial(n,i)((1,)*(n-i+1)) for i in range(1,n+1)]) ==
bell_number(n)
True
sage: n = 8
sage: add([bell_polynomial(n,i)((1,)*(n-i+1)) for i in range(1,n+1)]) ==
bell_number(n)
True
sage: n = 20
sage: add([bell_polynomial(n,i)((1,)*(n-i+1)) for i in range(1,n+1)]) ==
bell_number(n)
True
sage: bell_number(n)
51724158235372
}}}
Returns these pretty fast too!
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/5109#comment:3>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en
-~----------~----~----~----~------~----~------~--~---
