並前往https://www.mathworks.com/matlabcentral/answers/106715-how-to-add-a-constraint-condition-to-fsolve
