asube.m

有限元平面问题的求解相关程序。 包括质量矩阵

function y =Asube(K,k,tum)
n=size(tum,1);
for i=1:n
    for j=1:n
        K(tum(i),tum(j)) = K(tum(i),tum(j)) + k(i,j);
    end
end
y=K;