El sistema en el que estoy interesado tiene la siguiente forma:
$$dx_i = x_i\left(b_i-\sum_{j=1}^n a_{ij}x_j^2\right) dt$$
Para propósitos de programación, sería más fácil escribir este conjunto de EDOs en forma de matriz. Desafortunadamente, no he podido hacerlo. Sin el cuadrado, es fácil:
$$dx_i = x_i\left(b_i-\sum_{j=1}^n a_{ij}x_j\right) dt \implies dx = diag\left(x_1,\cdots,x_n\right)\left[b-Ax\right]dt$$
¿Existe alguna manera elegante de extender esto al caso con el cuadrado?
