How can I increase the output vector of a loop from a 7 x 7 double to a 31 x 7 double?

Hi everyone, I would like to use a for loop (employing Euler method) to solve ODEs and algebraic equations. However, loop only generate a vector of 7 x 7 double, and I am expecting a 31 x 7 double vector. Also, when I plot t-vector vs y-vector, I get an error message:

Error using plot
Vectors must be the same length.

Error in EulerMethod (line 62)
plot(t,y)

For the details, I have attached the m-file (EulerMethod.m). I have also attached the functions/files that it calls, namely Hion.m, C_Bulk.m, odes.m and Parameters.m.

I thank you.