function f=loadff
%----形成節點荷載
global   nfp_number nfp  node_number npe_number %（定義全局變量）
f = zeros(node_number*3,1) ;%（創建f）
for inf=1:1:nfp_number
  n = nfp( inf, 1 ) ;
  d = nfp( inf, 2 ) ;
  f((n-1)*6 + d ,1) = nfp( inf, 3 ) ;
end
fe=equload;
for i=1:1:3*node_number%%%% 6*node_number
    f(i,1)=f(i,1)+fe(i,1);
end
return