function f = pushing_force_param ( t, n, un, vn, k, fpush, tpush )
    f = zeros(n,1);
    f(k) = -fpush * (t <= tpush); % [o,tpush] ‚ÌŠÔ�CŠO—Í -fpush
end
