function?[]?=?FD_matrix（）
%基礎設置
L=10;
n=11;%單邊的點數量
h=L/（n-1）;%步長
%產生場源地圖
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）=