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: Marko.Vauhkonen@uku.fiw=[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'; endint=int*dJt;