偏微分方程matlab求解,求大神：如何使用Matlab求解偏微分方程的解析解

function [u,x,y] = poisson(f,g,bx0,bxf,by0,byf,D,Mx,My,tol,MaxIter)

%solve u_xx + u_yy + g(x,y)u = f(x,y)

% over the region D = [x0,xf,y0,yf] = {(x,y) |x0 <= x <= xf, y0 <= y <= yf}

% with the boundary Conditions:

% u(x0,y) = bx0(y), u(xf,y) = bxf(y)

% u(x,y0) = by0(x), u(x,yf) = byf(x)

% Mx = # of subintervals along x axis

% My = # of subintervals along y axis

% tol : error tolerance

% MaxIter: the maximum # of iterations

x0 = D(1); xf = D(2); y0 = D(3); yf = D(4);

dx = (xf - x0)/Mx; x = x0 + [0:Mx]*dx;

dy = (yf - y0)/My; y = y0 + [0:My]’*dy;

Mx1 = Mx + 1; My1 = My + 1;

%Boundary conditions

for m = 1:My1, u(m,[1 Mx1])=[bx0(y(m)) bxf(y(m))]; end %left/right side

for n =