[求助]二维非线性薛定谔方程的matlab求解 [复制链接]

光学空间孤子在非线性材料中传播时需遵循如下的“非线性薛定谔方程”： E>b2+;Jv  
!pTJ./  
     YB<*"HxM)}  
                       fmqb`%  
C+[%7vF1  
r={c,i  
为横向的拉普拉斯算符，A为光场的电场分量的慢变振幅，k=k0n0，k0是信号光的真空波矢，n为折射率，n0 为本底折射率，z为光束的传播方向。 `Z:R Ce^  
根据上式,我们可以得到描述光学空间孤子传播动力学的二维非线性薛定谔方程: GSclK|#tE  
5=Xy,hmnC  
7SDFz}  
m!<i0thJ  
其中 5oEV-6  
j/=iMq  
Δ(n2)=4πST0α dr)YzOvba  
T|J9cgtS  
α为一常数 *B`Zq)
  
求用matlab对此二维非线性薛定谔方程求解的matlab程序，加程序解释，求教。
NBl+_/2'w  
Q!- 0xlx  

同样想求教这个问题