Solving simultaneous equations involving complex numbers

[Helge Gabert]

In two dimensions:

fsolve([eq1, eq2],[var1, var2],[startvalue_for_var1, startvalue_for_var2])

eq1 could be x^y-y^2=9, var1 would be x, var2 would then be y, startvalue_for_var1 could be 0.1.

eq2 could be . . . you choose.

Just add more equations, variables, starting values for higher dimensions.

Starting values are important, because generally, your solution will not be unique.