ex0906.m

来自「偏微分方程matlab m文件」· M 代码 · 共 11 行

function ex0906()
A = 1e-3;
it0 = inline('0','x','y'); %(E9.3.2a)
bxyt = inline('x^2*sin(y)-y^2*sin(x)','x','y','t'); %(E9.3.2b)
D = [0 5 0 5]; T = 1000; Mx = 50; My = 50; N = 50;
[u,x,y,t] = TDE(A,D,T,it0,bxyt,Mx,My,N);
mesh(u);
xlabel('x')
ylabel('y')
zlabel('U')