可以用来解决一维溶质运移问题matlab代码