ElemStiffMatP1.m

Go to the documentation of this file.

function AElem=ElemStiffMatP1(q1,q2,q3,area)
% function AElem=ElemStiffMatP1(q1,q2,q3,area)
%   Computation of the element stiffness matrix for a triangle with
%   vertices 'q1,q2,q3'  and area 'area'
%
% Parameters:
%  q1: `2\times 1` array, first vertex of the triangle,
%  q2: `2\times 1` array, second vertex of the triangle,
%  q3: `2\times 1` array, third vertex of the triangle,
%  area: triangle area
%
% Return values:
%  AElem: `3\times 3` matrix
%
M=[q2-q3, q3-q1, q1-q2];
AElem=(1/(4*area))*M'*M;