function value = expo_int ( a, b )
%
%% expo_int returns the definite integral of exp(-x^2) over [a,b].
%
  value = 0.5 * sqrt ( pi ) * erf ( b );
  return
end
