function?[]?=?FD_matrix（）
%基礎(chǔ)設(shè)置
L=10;
n=11;%單邊的點(diǎn)數(shù)量
h=L/（n-1）;%步長
%產(chǎn)生場源地圖
x=linspace（0Ln）-3;
y=linspace（0Ln）-6;
[xxyy]=meshgrid（xy）;
rr=sqrt（xx.*xx+yy.*yy）;
F=zeros（nn）;
F（rr（:）<=1）=1;
x=linspace（010n）-7;
y=linspace（010n）-3;
[xxyy]=meshgrid（xy）;
rr=sqrt（xx.*xx+yy.*yy）;
F（rr（:）<=1）=-1;
%生成邊界圖
A=zeros（nn）;
A（:1）=1;
%邊界齊次化，邊界條件向源的影射
F2=F*（h^2）;
F2（2:）=F2（2:）-A（1:）;
F2（n-1:）=F2（n-1:）-A（n:）;
F2（:2）=