我在编程时想用分布差分法求解方程，对于三维的空间坐标，时间指标k就会是k——k+1/3----k+2/3----k+1这样分步实现，请教一下，编程时怎么才能办得到阿？
非常感谢各位！

哪位高手帮帮忙吧，小妹初涉matlab，很是迷茫！

先谢谢各位!
问题我没有说清楚，具体太复杂说不清楚，概括一下：方程是关于坐标x，y，z和时间t的非线性偏微分方程，在迭代求解时，可以从初始时刻t的值，求得在t+1/3时刻在x坐标方向的分布，再可由此求得t+2/3时刻在y坐标方向上的分布，然后可再由此求得t+1时刻在z坐标方向的分布，这样便有t时刻的值求出了t+1时刻的值，依此类推，可以迭代求得任意时刻方程的解。这个过程时间方向的迭代步长就是分数1/3，请教一下编程中怎么处理？