Есть метод ортогонализации Годунова решения при прогонке. Метод очень древний и применим только для обыкновенных дифуров - суть метода состоит в том,
Цитата:
что в процессе численного интегрирования ( например, на каждом шаге) делаются остановки, при которых восстанавливается линейная независимость решений.
. У Вас при малом шаге по границе те же проблемы.