有限体积法（二维线性元）

1. 齐次Dirichlet 边值条件
2. 矩形区域

function FVM_2d_mat_poisson

%%  code style : matrix
% finite volume method :piecewise linear element
% uniform mesh & circlecentre dual mesh
% u= cos(pi/2.*x).*cos(pi