function?[ufvf]?=?sfn（uv）
global?e?b?k;
ny?=?size（u1）;
nx?=?size（u2）;
uer?=?[u（:1）uu（:nx）];
uec?=?[u（1:）;u;u（ny:）];
ul?=?uec（3:ny+2:）+uec（1:ny:）+uer（:1:nx）+uer（:3:nx+2）-4*u;
u3?=?u.*u.*u;
uf?=?（u-u3/3-v）/e?+?k^2*ul;
vf?=?e*（u+b-0.5*v）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????330??2009-04-20?14:06??sfncont.m

?????文件????????258??2008-03-09?06:25??sfn.m

-----------?---------??----------?-----??----

??????????????????588????????????????????2