A number of new results on sufficient conditions for the solvability and numerical algorithms of the following general algebraic inverse eigenvalue problem are obtained: Given n+1 real n×n matrices A = (aij), Ak = (a(k)ij))(k = 1, 2 , n) and n distinct real numbers λ1, λ2 , λn, find n real numbers c1, c2,..., cn such that the matrix A(c) = A + n∑k=1 ckAk has eigenvalues λ1, λ2,..., λn.