stima3.m

来自「this fem program is written in matlab to」· M 代码 · 共 11 行

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 );