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