I try to translate a Maple function computing OEIS A144150 as given in [1]:
g := proc(p) local b;
b := proc(n) option remember; local k;
`if`(n=0,1,add(binomial(n-1,k-1)*p(k)*b(n-k),k=1..n))
end end:
A144150 := (n, k) -> (g@@k)(1)(n);
seq(print(seq(A144150(n,k),k=0..8)),n=0..5);
I am somewhat clueless how to achieve this with Sage. Google found [2] which highlights the problem but does not give an answer.
[1] http://oeis.org/A144150 [2] http://groups.google.com/group/sage-support/browse_thread/thread/9cda87e8889b6ead/
