function rhs=simple_rhs(t,y)
% the variable t must appear as an argument even when it is not
% explicitly used in the rhs of the ODE
global a
rhs=exp(-a*y.^2);