[Amplitude]=SquarePulse(t)

function [Amplitude]=SquarePulse(t)
%Input
%   t          (us)
%Output
%   Amplitude  (uA)
%   setting the Amplitude to empty indicates the end of the stimulus
%   pulse. This help speed up the calculation
%----------------------------------------
Pulsewidth  =   100;    %us
Delay       =   0;      %us
Amp         =   0.01;   %uA
%----------------------------------------

if t<=Delay
    Amplitude=0;
elseif t <= (Delay+Pulsewidth)
    Amplitude=Amp;
else
    Amplitude=[];
end