matlab 二维泊松 有限差分法 通过系数矩阵求解

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）=