force.m

纽约州立大学的数学建模和算法的教程

function [f] = force(tt,P) 

% Periodic forcing with period P.

ttt=tt;
while ttt > P,
  ttt=ttt-P;
end

f=(ttt<(P/8));