Hi, all For a vector of 1800 variables, I have 1800 equations like y(i)'=y(y1,y2,...,yn,t), but RHSfunction only have 5 y(j), so if I put related y(i) next to each other, the RHSfunction will become more banded, and Jacobian Matrix will be more banded also, the time consumption will go down. Am I right? Thanks! Xuan YU xxy113 at psu.edu