bound.m

是个好东西

function bound
%  -----  边界条件
global node element material K bc bc_number  %（定义全局变量）
if bc_number==0
    return
end
for ibc=1:1:bc_number
  n = bc(ibc, 1 ) ;
  d = bc(ibc, 2 ) ;
  m = (n-1)*3 + d ;
  f(m) = bc(ibc, 3)* K(m,m) * 1e20 ;
  K(m, m) = K(m, m) * 1e20 ;
end
return