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