function [y]=inital(neqn,t)
%
% Function inital sets the initial condition
% for the ODE with erf as the solution
%
% Initial condition
  y(1)=0.0;

