A two-dimensional time domain finite difference beam propagation method, ba
sed on the slowly varying envelope approximation is presented. Expressions
to correct this approximation and to include third-order non-linear effects
are given. The method is applied on solitary waves in a non-linear Bragg r
eflector, assuming realistic materials parameters.