I'm sorry if this is a weird question, but I have a system of differential equations of the form \begin{align} \begin{bmatrix} S_1'=f_1(S_1, S_2, S_3,...S_n) \\ S_2'=f_2(S_1, S_2, S_3,...S_n) \\ S_3'=f_3(S_1, S_2, S_3,...S_n) \\ \vdots \\ S_n'=f_n(S_1, S_2, S_3,...S_n) \end{bmatrix} \end{align}
where $f_i(S_1, S_2, S_3,...S_n)$ is a non-linear function that is a type of logistic model equation. As such it currently stands that $0\leq S_1, S_2, S_3, ..., S_n \leq 1 \forall t>0$.
However I was wondering if I could restrict these differential equations even harder to the point in which all of the values of $S_k$ could sum to $1$, or written in math $$ \sum_{k=1}^n S_k(t) = 1 $$
I looked around, but I couldn't find anything like this. Do any of you have an idea on what I could possibly do to solve this?
Note I only am doing this numerically, an explicit solution to all equations will not be required, so if something is only numerically solvable, that is still good
Note that your constraint defines some codimention 1 surface on which your solutions must lie. In physics we call this a holonomic constraint, i'm sorry, but I am not sure whether this is a right mathematical name. Such a constraint may be used to reduce the dimension of the problem. This is done very easily. You have the constraint $$\sum_{k=1}^{n}S_{k}=1$$ Differentiate the constraint to give $$\sum_{k=1}^{n}S_{k}'(t)=\sum_{k=1}^{n}f_{k}(S_{1}, ..., S_{n})=0$$ This thing in principle can be solved for one of the $S$'s. Say, $S_{n}$, that is $$S_{n}=F(S_{1}, ..., S_{n-1})$$ Thus you have to solve $$S_{1}'=f_{1}(S_{1}, ..., S_{n-1}, F(S_{1}, ..., S_{n-1}))$$ $$....$$ $$S_{n-1}'=f_{n-1}(S_{1}, ..., S_{n-1}, F(S_{1}, ..., S_{n-1}))$$ One equation less. \ P.S. this principles are the hurt of classical mechanics. In classical mechanics, the constraints may be imposed as in your problem, or are due to intrinsic symmetries. If the number of them is $n-1$ for a system with $n$ degrees of freedom, the system is said to be completely integrable, i.e. it can be reduced to a single separable ode, which is always soluble by quadratures. A famous example of scientific importance is a two body problem, the $6$ dimensional configuration space possesses just right number of symmetries to reduce the problem to $1D$ and thus integrate it. Whereas a 3 body problem doesn't....