我积极编写的matlab M文件用来求解菲克第二定律的扩散