No solution found.
fsolve stopped because the problem appears regular as measured by the gradient,
but the vector of function values is not near zero as measured by the
default value of the function tolerance.
<stopping criteria details>
fsolve stopped because the sum of squared function values, r, has gradient with
relative norm 7.341307e-07; this is less than options.OptimalityTolerance = 1.000000e-06.
However, r = 7.954251e-01, exceeds sqrt(options.FunctionTolerance) = 1.000000e-03.
Optimization Metric Options
norm(grad r) = 7.34e-07 OptimalityTolerance = 1e-06 (default)
r = 7.95e-01 sqrt(FunctionTolerance) = 1.0e-03 (default)