Is it possible to use any custom optimization method from standard MATLAB's toolkits? It seems not to be so difficult in comparison with lsqnonlin().