force.m

来自「纽约州立大学的数学建模和算法的教程」· M 代码 · 共 11 行

function [f] = force(tt,P) 

% Periodic forcing with period P.

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

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