I would like to solve a system of differential equations with the command

 p:=dsolve([sys,start], type = numeric, range = 0..10)

sys are containing the  differential equations and start the initial conditions. The equations are of the form

 f _n ' (t) = sum for i = n to 1000 of f _i (t)

where n is between 1 and 1000. My problem is that I can’t figure out how to define sys in this case. Probably I should use a for-loop, but can someone help me with the details?

Thanks in advance

Torben Svendsen