На базе непрерывного аналога метода Ньютона и метода сплайн-коллокации создан пакет программ для решения краевых задач для нелинейных систем обыкновенных дифференциальных уравнений (ОДУ), заданных на вложенных интервалах изменения независимой переменной. В качестве конкретного примера приложения рассматривается задача о статических распределениях магнитных потоков в геометрически симметричном двухслойном джозефсоновском стеке, отдельные субконтакты которого имеют разную длину. Работа выполнена в Лаборатории информационных технологий ОИЯИ.