Is there an algorithm to find a family of polynomials with given finite set of zeros? Question: Is there an algorithm for finding a collection $S$ of multivariate polynomials which has a pre-specified finite set of zeros in $\mathbb{C}^n$, for $n>1$? 
E.g., how to find a set of bivariate polynomials whose only common zeros are $z_1 =(1,2)$ and $z_2 = (3,4)$, or how to find a set of trivariate polynomials whose only common zero is $z_1 =(1,1,1)$.
A pointer to a reference would suffice for an answer since this question is almost assuredly incredibly basic and not really worth anyone's time to answer in-depth.
Attempt: It's obvious how to do this for $n=1$, i.e. polynomials with coefficients in $\mathbb{C}$ (i.e. just use $S=\{(x-z_1)(x-z_2)\dots(x-z_m)\}$). But the analogous trick doesn't work in higher dimensions, even if the set of zeros consists of only one point. For example, $(x-1)(y-1)(z-1)$ has more zeros than just the point $z_1=(1,1,1)$ -- any point for which any of the $x,y,z$ coordinates is equal to $1$ is also a zero, whereas we are only interested in the point for which all three are equal to 1.
My solutions for these problems are increasingly ad hoc in higher dimensions and I feel like I might be wasting my time by not using a general procedure. For example, I am able to check my answers using Bezout's theorem and the properties of intersection multiplicity, but it would be easier not needing to check any answers because I knew I was using a provably correct algorithm.
Note: Actually I am interested in $\mathbb{P}^n(\mathbb{C})$, but it isn't difficult to use a solution for the problem for $\mathbb{C}^n$ to solve it for $\mathbb{P}^n(\mathbb{C})$, so for the sake of simplicity I am only asking about the problem for $\mathbb{C}^n$.

Let your desired set of zeros be $\{z_1, \ldots, z_m\}$.
Applying a linear transformation, we may assume the first coordinates $(z_j)_1$ are all distinct.
We can construct interpolation polynomials $P_k$ so that $(z_j)_k = P_k((z_j)_1)$ for $k = 2 \ldots n$.  Then take your set of polynomials to be 
$\prod_{j=1}^m (Z_1 - (z_j)_1)$ and $Z_k - P_k(Z_1)$ for $k = 2 \ldots n$.
Thus in your first example ($(1,2)$ and $(3,4)$) the polynomials could be $(x-1)(x-3)$ and $y - x - 1$.
