function r = rhs_erf(t,y)

r = 2/sqrt(pi) * exp(-t^2);
