近期开发了一套数值解一个方程的算法。方程是非线性光学中光传播的非线性薛定谔方程，方程具有以下形式：

其中U是关于时间t, 传播距离z, 横向(x, y), T=(1-i∂_t/ω_0)算符，D\hat为时间色散项k^((2) )/2! (i∂_t )^2。 n2, xk, k0等等为常数。

遇到类似问题的朋友可以交流。