Suppose I have a vector equation:
     
    Y'[t]==rhs[Y[t]] 
and

    Y[0]==ConstantArray[0,n]

where "rhs[Y[t]]" is a black box function which numerically calculates a N-length vector given calculates given a lenght-N vector Y[t].

How can I solve such an equation in "NDSolve" ?