Асимптотика решения разностного уравнения с малыми шагами в прямоугольной области

Строится асимптотическое разложение решения разностного уравнения
\begin{align*} &u(x,t)=\alpha(x,t)u(x-h,t-\tau)+\beta(x,t)u(x,t-\tau)+\\ &+\gamma(x,t)u(x+h,t-\tau)+f(x,t)+\\ &+hf_1(x,t,u(x,t),u(x-h,t-\tau),u(x,t-\tau),u(x+h,t-\tau))+\\ &+\tau f_2(x,t,u(x,t),u(x-h,t-\tau),u(x,t-\tau),u(x+h,t-\tau)) \end{align*}
с дополнительными условиями
$$ u(0,t)=u_1(t),\qquad u(l,t)=u_2(t),\qquad u(x,0)=u_0(x) $$
в виде суммы ряда по степеням $h$ и $x$ и ряда из функций пограничного слоя.