intet.m

别人的边界元程序

function zx=intet(numbs,xp,yp,zp)

%% calculate the potential of inner point 
global gi ome
global fjacob cosbx cosby cosbz
global yy
global mtype nord
global x y z
global icod kcd fi tk tq tkt tqt ml
global at bt sum
global l cx cy cz sol
global tp

xk=zeros(4,1);
yk=xk;
zk=xk;
nd=xk;
%%
tp=0;
for j=1:numbs
    sum=0;
    at(:)=0;
    bt(:)=0;
    for i=1:4
        nd(i)=i;
        xk(i)=x(nord(j,i));
        yk(i)=y(nord(j,i));
        zk(i)=z(nord(j,i));
    end
    coft(nd,xp,yp,zp,xk,yk,zk);
    for i=1:4
        tp=tp+bt(i)*tq(j,i)-at(i)*tk(j,i);
    end
end
%%
zx=1;
return