function [yt]=derv(neqn,t,y)
%
% Function derv computes the derivative
% of the ODE with erf as the solution
%
% ODE
  yt(1)=exp(-t^2);
