2.4.3.3 Mathematica
ode=x*D[y[x],x]+y[x]==0; 
ic=y[0]==1; 
AsymptoticDSolveValue[{ode,ic},y[x],{x,0,5}]
 

Not solved