Ein[z_] := NIntegrate[(1 - Exp[-t])/t, {t, 0, z}];
EinV[z_] = ExpIntegralE[1, z] + Log[z] + EulerGamma;
SER = FullSimplify[Series[EinV[z], {z, -1, 2}], Assumptions -> z \[Element] Reals] // Normal;
Plot[{Ein[z], EinV[z], SER}, {z, -1, 1},
PlotStyle -> {Red, {Dashed, Black}, Blue},
PlotLegends -> "Expressions"]
For z complex values:
EinV1[z_] = ExpIntegralE[1, z] + Log[z] + EulerGamma;
SER = FullSimplify[Series[EinV1[z], {z, -1 + 2*I, 10}]] // Normal;
Plot[{EinV1[z] // ReIm // Evaluate}, {z, -1, 1}, PlotLegends -> "Expressions"]
Plot[{SER // ReIm // Evaluate}, {z, -1, 1}, PlotLegends -> "Series"]
Regards M.I.