Poisson 泊松方程的差分方法matlab实现(2)

Poisson 泊松方程的差分方法matlab实现

源代码：

function F=fivepointdiff(l,n)

h=l/n;

N=2*(n-1)*n+(3*n-1)*(n-1);

XY=zeros(2,N);%分割xy轴后每一个节点的坐标

for i=1:n

for j=1:n-1

XY(:,(n-1)*(i-1)+j)=[l+j*h;i*h];

end

end

for i=1:n-1

for j=1:3*n-1

XY(:,n*(n-1)+(3*n-1)*(i-1)+j)=[j*h;l+i*h];

end

end

for i=1:n

for j=1:n-1

XY(:,n*(n-1)+(3*n-1)*(n-1)+(n-1)*(i-1)+j)=[l+j*h;2*l+(i-1)*h]; end

end

A=zeros(N,N);

for i=1:N

for j=1:N