The determination of a large class of regular and singular Euclidean self-dual Yang-Mills field configurations is reduced to the solution of a set of linear algebraic equations. The matrix of the coefficients is a polynomial functions ofx and the rules for its construction are elementary.