How to solve 'The Levenberg-Marquardt algorithm does not handle bound constraints' error

Question

I am trying to use LM Algorithm to optimize using lsqnonlin instead of interior-point to allow for gradient free optimization. However I get the following error

Error using lsqncommon (line 64)
The Levenberg-Marquardt algorithm does not handle bound constraints and the trust-region-reflective algorithm requires at least as
many equations as variables; aborting.

Even though it is solved here as follows 1, This says as LM algorithm can be modified to solve the problem in hand with boundaries2. Is it possible to solve this error by modifying algorithm using lsqnonlin function?