In the original paper defining the Seiberg-Witten map, I have been confused about the following step in their derivation. Using the gauge transformation constraint, they write \begin{align*} A'_i (A+ \delta_\lambda A) - A'_i(A) - \partial_i \lambda' - i[\lambda', A_i] - i[\lambda, A'_i]= -\frac{1}{2}\theta^{kl} \{\partial_k \lambda, \partial_l A_i\} \end{align*} where $A_i,A'_i, \lambda,\lambda'$ are matrices, and the brackets and braces correspond to commutator and anti-commutator, resp. They then proceed to solve this equation for $\lambda'$ and $A'_i$, to obtain \begin{align*} A'_i(A) =& -\frac{1}{4}\theta^{kl}\{A_k, \partial_l A_i + F_{li}\} \\ \lambda'(\lambda, A) =& \frac{1}{4}\theta^{ij}\{\partial_i \lambda, A_j\} \end{align*} where $F_{li}$ is the standard field strength. This corresponds to Eq. 3.4 and Eq. 3.5 in the original paper.
How did they obtain this result? In particular, how did they get a solution for two variables from a single equation?