Is there anything like this?
sol1 = Table[
NDSolveValue[{D[xtraj[t],
t] == -(Sinh[2 xtraj[t] (-t)]/(Cosh[2 xtraj[t] (-t)])),
xtraj[0] == n}, xtraj[t], {t, -4, 8}], {n, -.1, .1, .02}]; p1 =
Plot[sol1, {t, -4, 8}, PlotRange -> {{-1, 1}, {-.15, .15}},
PlotStyle -> Automatic, AxesStyle -> Thickness[.001],
LabelStyle -> {Black, Medium},
PlotLegends -> Table[n, {n, -.1, .1, .02}]]