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.