А как найти оценку граничного условия?
По общему правилу: недостающие начальные значения (в Вашем случае это

) надо аппроксимировать с тем же порядком точности, что и глобальный порядок точности метода. Т.е. достаточно для их нахождения использовать метод на единицу меньшего порядка, что у основного метода. Соответственно, в Вашем случае достаточно использовать для нахождения

метод Эйлера (случайно, кстати, решение получится абсолютно точным). Условие же

откровенно имеет порядок погрешности

(даже не вдаваясь в детали -- просто потому, что

, а производная в нуле ненулевая) и, следовательно, не подходит.