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

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

if(i==j)

A(i,j)=4;

else if(((XY(1,i)-XY(1,j))^2+(XY(2,i)-XY(2,j))^2)<2*h*h)%若是相邻点 择系数为-1

A(i,j)=-1;

end

end

end

end

f=zeros(N,1);%就是等号右边F

for i=1:N

f(i,1)=h*h;

end

U=bicg(A,f,0.1,100);%求解Au=F

F=[XY;U'];%输出

命令框中输入：

fivepointdiff(1,25);

x=ans(1,:);

y=ans(2,:);

z=ans(3,:);

plot3(x,y,z)

得到的结果：