It seems, like you are not using ODE45, but ODE15S.
Currently you divide the simulation in 401 parts by using [tstart:tspan:tfinal]. (btw. you can omit the square brackets. a:b:c is a vector already). Now you want to divide it to 2000 seconds. If your tfinal is measured in seconds, why not simply use tspan=2000 ?
