boundaryedges.m

this fem program is written in matlab to simulate a room temperature of a museum under the source of

function [dirichlet,neumann]= boundaryedges(coordinates, e);

%used to get dirichlet and nueuman edges
externaledge = e([1,2],:)';
dd=[];%zeros(size(externaledge,1),4);
for i = 1:size(externaledge,1)
    dd(i,[1:4])=[externaledge(i,:),sum(coordinates(externaledge(i,:),1))./2,sum(coordinates(externaledge(i,:),2))./2];
end
dd;
dirichlet = [dd(find(dd(:,3)==0),[1:2])];

neumann = [dd(find(dd(:,3)==1),[1:2]);dd(find(dd(:,4)==1),[1:2]);dd(find(dd(:,4)==0),[1:2])];