用MATLAB写了一个解扩散方程的程序（用的是老师给的程序），可是运行显示出错，哪位高手帮帮忙！！

%只修改了一点语法错误，修正了向量维数，没对数据做任何修改： function demo_MOL clc;clear all;format long; n=15;dpi=pi/n; x=dpi:dpi:pi; size(x) u=sin(x); t=0.:0.4:40.; [t u]=ode45(@myfun,t,u); uu(:,2:n)=u(:,1:n-1); uu(:,1)=0;uu(:,n+1)=0; x=[0. x];[xx yy]=meshgrid(x,t);