In[]:=
s=NDSolve[{y'[t]-{1,$MachineEpsilon}*y[t],y[0]{1,1}},y,{t,0,1}]
Out[]=
yInterpolatingFunction
Domain: {{0.,1.}}
Output dimensions: {2}
