faster version of "fsolve" function

Does anybody know a function in MATLAB that can perform a similar job as "fsolve" does but faster? I have a user-defined function in Simulink and I use "fsolve" to solve a set of nonlinear equations. The number of unknown variables equal to the number of equations. The system I want to solve is a big system (168 variables and equations) and the simulation is too slow takes more than 40 min. for each run with acceptable accuracy. Any idea? I found this link but did not try it yet. http://www.walkingrandomly.com/?p=1488 Thanks