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

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

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

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

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

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