I have $N$ differential equation for each $s \in {1,2,...,N}$ on the domain $x \in [1,2]$ of the form (e.g.)
$xf_s'(x) + \log(xs) f_s(x)+1 = 0$, with the boundary condition $f_s(1) = e^{-s}$
Now I want to solve this equation with NDSolve. Put the values for the solutions at $x = 2$, i.e. $f_s(2)$ in a table enumerated by the index $s$. Then I want to interpolate between these values in the table to get a function $g(s) = f_s(2)$.
However, since I am new to mathematica, I have no Idea how to do this. I can solve the equation for some $s$ with NDSolve:
NDSolve[{x*D[f[x], x] + Log[x*s]*f[x] + 1 == 0, f[1] == Exp[-s]}, f, {x, 1, 2}]
and this seems to work. However my main issue right now is putting the values of the output function in the table. How can I do this?