The ode routines assume a vector input y0 ( 1xn ) and give a vector output y . Immediately after the call to ode15s you can convert back to a matrix, using e.g. Y = reshape(y,3 2) and the rest of your code should work.
I cant really gove you an answer,but what I can give you is a way to a solution, that is you have to find the anglde that you relate to or peaks your interest. A good paper is one that people get drawn into because it reaches them ln some way.As for me WW11 to me, I think of the holocaust and the effect it had on the survivors, their families and those who stood by and did nothing until it was too late.