ex0902.m

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

f = inline('x^2+y^2','x','y'); 
g = inline('sqrt(x)','x','y');
x0 = 0; xf = 4; y0 = 0; yf = 4;%自变量取值范围
Mx = 50;My = 30;%等分段数
dbx=inline('x^2+y^2','x','y');
bx0 = inline('y^2','y'); %边界条件
bxf = inline('4^2*cos(y)','y'); 
by0 = inline('x^2','x'); 
byf = inline('4^2*cos(x)','x'); 
D = [x0 xf y0 yf]; MaxIter = 100; MinErr = 1e-4;
[U,x,y] = Helmholtz_Newton(f,g,dbx,bx0,bxf,by0,byf,D,Mx,My,MinErr,MaxIter);
clf, mesh(U)
xlabel('x')
ylabel('y')
zlabel('U')