Given a set of $n$ vectors in a $n$-dimensional vector space $\{u_1, u_2,...,u_n\}$ I am interested in a method to get a $n \times n$ matrix $A$ whose first $n$ powers have each of those vectors as rowsum which means that for $1\leq k\leq n$ the matrix $A$ verifies
$$A^k \begin{align} \begin{bmatrix} 1 \\ 1 \\ \vdots \\ 1 \\ 1 \end{bmatrix} \end{align} = u_k.$$
This isn't always guaranteed. You need additional constraints on what your vectors can be.
Counterexample Consider the vectors $$ u_1 \;\; =\;\; \left [ \begin{array}{c} 1\\1\\ \end{array} \right ] \hspace{3pc} u_2 \;\; =\;\; \left [\begin{array}{c} -1 \\ 1 \\ \end{array} \right ]. $$
Then the we want a $2\times 2$ matrix $A$ that satisfies $Au_1 = u_1$ and $A^2u_1 = u_2$. By the first constraint we want that $$ \left [\begin{array}{cc} a & b \\ c & d \\ \end{array} \right ]\left [ \begin{array}{c} 1\\1\\ \end{array} \right ] \;\; =\;\; \left [ \begin{array}{c} a+b \\ c+d \\ \end{array} \right ] \;\; =\;\; \left [ \begin{array}{c} 1\\1\\ \end{array} \right ]. $$
When we try to find $A^2$ to obtain our second vector we find that \begin{eqnarray*} \left [\begin{array}{cc} a & b \\ c & d \\ \end{array} \right ]\left [\begin{array}{cc} a & b \\ c & d \\ \end{array} \right ]\left [ \begin{array}{c} 1\\1\\ \end{array} \right ] & = & \left [\begin{array}{cc} a^2 + bc & b(a+d) \\ c(a+d) & d^2 + bc \\ \end{array} \right ]\left [ \begin{array}{c} 1\\1\\ \end{array} \right ] \\ & = & \left [ \begin{array}{c} a^2 + b(a+c+d) \\ d^2 + c(a+b+d) \\ \end{array} \right ] \\ & = & \left [ \begin{array}{c} a^2 + ab + b \\ d^2 + cd+c \\ \end{array} \right ] \\ & = & \left [ \begin{array}{c} a+b \\ c+d \\ \end{array} \right ]\;\; =\;\;\left [ \begin{array}{c} 1\\1\\ \end{array} \right ]. \end{eqnarray*}
Therefore no matrix $A$ exists where $Au_1 = u_1$ and $A^2u_1 = u_2$.