Rayleigh-Taylor 不稳定建模的困惑

关于最近Matlab解方程建立模型的困惑

最近在利用Matlab模拟Rayleigh-taylor 不稳定性原理的产生和发展过程，通过一步步解质量守恒方程，poison方程和动量守恒方程，单独解各个方程的时候并无大问题，但根据需要，一个方程的解需要带入下一个方程，得到第二个解再带入第三个方程，然后又代回第一个方程，如此一次循环，得到不断更新的解。在这个过程中，怎么确保每个方程的解是正确的，很伤脑筋，因为这是环环相扣的过程，前一个解的正确性需要后一个去验证，而后一个的解又需要新的解去验证，在这种情况下，有好的办法去验证自己的程序中间有没有出错吗？