IGEM:Melbourne/2008/BioClock/David Version/SGO pulse

function sg0 = getSG0(time)

global tau_on tau_off;

%starting time

% %%% exponentially decaying signal % % time_0 = 0; % lambda = 0.05; % % if time < time_0; % sg0 = 0 ; % end % % if time >= time_0 % sg0 = exp(-lambda*(time-time_0)); % end

% %%% pulsed signal % t_ratio=time/(tau_on+tau_off); j= floor(t_ratio)+1; if (time >=(j-1)*(tau_on+tau_off)) && (time <= j*tau_on+(j-1)*tau_off)

  sg0 = 1;   

elseif (time >= j*tau_on + (j-1)*tau_off) && (time <= j*(tau_on+ tau_off))

  sg0 = 0;

%else % DPTH(k) = 'error' % stop end

end