I want to solve a large system of 500 by 500 differential equations such that I have equations such as:
Note all the
A2 are just constants, that are not necessarily equal.
y1'[t]= A1*y1[t]+ B1*y2[t]+...+Z1*y500[t] y2'[t]= A2*y1[t]+ B2*y2[t]+...+Z2*y500[t] . . . y500'[t]=A500*y1[t]+ B500*y2[t]+...+Z500*y500[t]
In the Mathematica documentation, I have found some tutorials on how to do this for a 2x2 matrix, but not for large matrices.