matlab 二维泊松 有限差分法 逐次迭代求解

function [] = FD（）
%逐次更新发
%场源图
% n=52;
% h=0.1;
% A=zeros（nn）;
% A（1:）=10000;
% F=zeros（nn）;
% B=A;
% k=10000;%迭代次数
% A=FD2（AFh1000）
% surf（A）;

n=52;
h=0.2;
A=zeros（nn）;%边界图
x=linspace（010n）-3;
y=linspace（010n）-6;
[xxyy]=meshgrid（xy）;
rr=sqrt（xx.*xx+yy.*yy）;
F=zeros（nn）;%场源图
F（rr（:）<=1）=1;
x=linspace（010n）-7;
y=linspace（010n）-2;
[xxyy]=meshgrid（xy）;
rr=sqrt（xx.*xx+yy.