Solved
When you defined D(t,X) you type a square bracket for the second entry and so you created inadvertently a 1x1 matrix .
Simply use 'normal' parentheses and its works:

BTW, there is no need to use a solve block when you use one of the stand alone ODE solvers. A solve block is needed only when you use 'odesolve'.
You may get a smoother result if you either increase N significantly or use a smaller interval [Ti; Tf], or both:

Enter your E-mail address. We'll send you an e-mail with instructions to reset your password.



