В общем, у меня граничные условия задаются гиперболическими функциями (tanh + tanh). По факту это потенциальная яма, которая может принимать любую форму. Выбор функций в целом был произволен, мне просто нужен был потенциальный колодец с локализованными силами. Я сейчас решил причесать свой код и оказалось, что простая смена операций при вычислениях даёт большую разницу. Я не могу сказать, что меня это сильно волнует: у меня качественная модель с ODE первого порядка. Но мне хотелось бы на будущее знать, как работать с подобными функциями, чья область значений меняется в очень больших пределах.