risingfactorial(m_,n_):=if(m<0,if(m+n>0,0, (-1)^n*(-m)!/(-m-n)!), gamma(m+n)/gamma(m) ); 
risingfactorial(m_,1):=m;