One last question, if you do not mind. Is there any other “recommended” method to handle this kind of problems, other than the Julia code snippet I posted? I was able to find an example in the documentation of the package that had an external function, but it did not have other parameters on which the solution depended, so I was forced to take a different approach.