finite difference time domain method f

Hello everyone, 

Does anyone have a suggestion for a good starting point for implementing FDTM method in Igor? Specifically I am interested in solving 2D GPE coupled equations. I have the 0d model working fine with  and I want now to add the spatial dependence. 

any help much appreciated. 

Cheers