Inside a xyz fsolve function, I used a good initial point that needed to be used but now I am struck here getting this
fsolve completed because the vector of function values at the initial point
is near zero as measured by the value of the function tolerance, and
the problem appears regular as measured by the gradient.
<stopping criteria details>
Equation solved. The final point is the initial point.
The sum of squared function values, r = 1.067653e-22, is less than sqrt(options.FunctionTolerance) = 1.000000e-03.
The relative norm of the gradient of r, 1.038534e-11, is less than options.OptimalityTolerance = 1.000000e-06.
Now, the query is that I want the fsolve function to solve all steps irrespective of how little the value we get and not stop the function and give a constant output. So, what change in that function shoul be implemented keeping the above code output lines in mind.