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
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 (2012-01-29, VBScript, 1KB, 下载4次)