求解泊松方程的function Finite_element_tri(Imax)% 用有限元法求解三角形形区域上的Possion方程Jmax=2*Imax;% 其中Imax Jmax分别表示x轴和y轴方向的网格数，其中Jmax等于Imax的两倍% 定义一些全局变量global ndm nel na % ndm 总节点数% nel 基元数% na 表示活动节点数V=0; J=0;X0=1/Imax;Y0=X0;%V=0为边界条件domain_tri % 调用函数画求解区域X,Y,NN,NE=setelm_tri(Imax,Jmax); % 给节点和三角形元素编号，并设定节点坐标% 以下求解有限元方程的求系数矩阵T=zeros(ndm,ndm);for n=1:nel n1=NE(1,n); n2=NE(2,n); n3=NE(3,n);%整体编号 s=abs(X(n2)-X(n1)*(Y(n3)-Y(n1)-(X(n3)-X(n1)*(Y(n2)-Y(n1)/2;%三角形
