Symbolic Vector ODE numerical solution

I'm trying to solve an ODE using vectorized symbolic calculus.

N = 4;
syms theta(t) [1 N]

Q = tril(ones(N));
theta = theta*Q';
eq = diff(theta, t, t) - cos(theta);
Dt = diff(theta, t);

How can I use this with ode45 or other solvers?

[Walter Roberson]

To use that with ode45 or similar, I recommend reading the first example for odeFunction() as it shows the steps to use to prepare symbolic differential equations for use with ode45 or similar.

Answer 1

[Paresh yeole]

in your equation for 'eq', if you are looking for double differentiation then it should be

eq = diff(theta, t, 2) - cos(theta);

[Walter Roberson]

diff(theta, t, 2) and diff(theta, t, t) are treated the same, both become diff(theta, t, t) on output.