0001 function areas=area(md)
0002
0003
0004
0005
0006
0007 areas=zeros(md.numberofelements,1);
0008
0009
0010 index=md.elements;
0011 x=md.x; y=md.y; z=md.z;
0012 x1=x(index(:,1)); x2=x(index(:,2)); x3=x(index(:,3)); y1=y(index(:,1)); y2=y(index(:,2)); y3=y(index(:,3));
0013
0014
0015 if strcmpi(md.type,'3d')
0016
0017 thickness=mean(z(index(:,4:6)))-mean(z(index(:,1:3)));
0018 areas=abs(0.5*((x2-x1).*(y3-y1)-(y2-y1).*(x3-x1)));
0019 else
0020
0021 areas=abs(0.5*((x2-x1).*(y3-y1)-(y2-y1).*(x3-x1)));
0022 end