stima3.m

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

function M = stima3 ( vertices )

%% STIMA3 determines the local stiffness matrix for a triangular element.


  d = size ( vertices, 2 );

  D_eta = [ ones(1,d+1); vertices' ] \ [ zeros(1,d); eye(d) ];

  M = det ( [ ones(1,d+1); vertices' ] ) * D_eta * D_eta' / prod ( 1:d );