Usage is as fsolve. The difference is that it can compute solutions in C^m for m-variables. Input can be both anonymous functions and symbolic equations. Subsequent calls to this function in a matlab session will not cause the recomputation of the internal matlab functions unless the input equations have change.
Corrected Line82: as options=optimset(options, ...