bound2.m

实现对电磁层析模型的建立

function int=bound2(g);

%bound2 Computes the boundary integral of the product of two linear basis function in 2D FEM
% Function int=bound2(g) calculates the boundary integral
% of the product of two basis function from g(1,:) to g(2,:).
%
% INPUT
%
% g = integration endpoints
%
% OUTPUT
%
% int = value of the integral

% 10.5. 1996 P. Ronkanen and M. Vauhkonen
% University of Kuopio, Department of Applied Physics, PO Box 1627,
% FIN-70211 Kuopio, Finland, email: [e-mail]


w=[1/2,1/2];
ip=[1/2-1/6*sqrt(3),1/2+1/6*sqrt(3)];
dJt=sqrt((g(2,1)-g(1,1))^2+(g(2,2)-g(1,2))^2); 
int=0;
 for ii=1:2
  S=[1-ip(ii);ip(ii)];
  int=int+w(ii)*S*S';
 end
int=int*dJt;