...but unsure on the syntax for a problem like this.
The syntax does not change. It is simply that your equation expressions will be functions of a vector variable z instead of a scalar.
F = @(z)[z(1)+z(2); 2*Z(1)+z(2)^2];
[x,y] = fsolve(F,[z1,z2]);
