d[1,m_]/;m≥1:=0;
d[_,0]=1;
d[n_,m_]:= d[n,m]=(n-1)d[n,m-1]+d[n-1,m];

Table[d[m,n],{n,10},{m,10}] // MatrixForm