function [y]=Heaviside(t,t0);  
y=[(t-t0)>=0];
%定義一個階躍函數
%使用說明：本函數要求給出一個t的范圍，以及t0這樣一個階躍點
%如果你要u(t)，那么寫成Heaviside(t,0)
%如果你要u(t+3),那么寫成Heaviside(t,-3),因為這個階躍點為-3
%如果你要u(-t+3),那么寫成Heaviside(-t,-3),
%也就是說t的符號不變，而值的符號變化。