test5.m

matlab中PDE工具箱的帮助文件pde.pdf的例子全部用代码实现,花了我一个月的时间,有限元应用程序

[p,e,t]=initmesh('circleg','Hmax',1);
error=[]; err=1;
while err > 0.001,
[p,e,t]=refinemesh('circleg',p,e,t);
u=assempde('circleb1',p,e,t,1,0,1);
exact=-(p(1,:).^2+p(2,:).^2-1)/4;
err=norm(u-exact',inf);
error=[error err];
end
pdemesh(p,e,t)
pdesurf(p,t,u)
pdesurf(p,t,u-exact')