说明：  Boundary element method function [index,g,h,a,unb,x,u,un,uin]=bemlabecom(n,c,index,in,cin) [index,g,h,a,unb,x,u,un,uin]=bemlabecom(16,[0 0 0.25 0 0.5 0 0.75 0 1 0 1 0.25 1 0.5 1 0.75 1 1 0.75 1 0.5 1 0.25 1 0 1 0 0.75 0 0.5 0 0.25 0 0],[1 0 1 0 1 0 1 0 0 200 0 200 0 200 0 200 1 0 1 0 1 0 1 0 0 100 0 100 0 100 0 100],9,[0.25 0.25 0.50 0.25 0.75 0.25 0.25 0.50 0.50 0.50 0.75 0.50 0.25 0.75 0.50 0.75 0.75 0.75]) in bemlabecom2 we change a & unb to : a1 a2,unb1 unb2 n=number of nodes=number of element=number of point c=matrix that determine the coordinate of point m=coordinate of nodes in global coordinates index: determine u or un for boundary nodes in:number of internal nodes cin: coordinate of internal nodes