Matlab求解一阶多元偏微分方程组的欧拉显式方法

求问文献为
[1] Development of an efficient buoyant jet integral model of a bubble plume coupled with a population dynamics model for bubble breakup and coalescence to predict the transmission loss of a bubble curtain
[2] Study of the Sound Escape with the Use of an Air Bubble Curtain in Offshore Pile Driving
背景是水中喷嘴生成的气泡羽流形成过程的求解问题，需求解有六个分量的向量u，已知关于m的六个等式，和关于q的六个等式，文献使用方法是欧拉显式法，但是欧拉显式不是求解常微分方程的吗，这里偏微分方程如何求解
方程为

各向量分量