% 接続行列と慣性行列の計算例

point = [0,0; 1,0; 2,0; 0,1; 1,1; 2,1];
triangle = [1,2,4; 5,4,2; 2,3,5; 6,5,3];
h = 2;

[Jlambda, Jmu] = connection_matrices(point,triangle,h);

Jlambda
Jmu

rho = 12;

M = inertia_matrix(point,triangle,h,rho);

M
