I have the equation

m*D2y +dc*Dy + sk*y = A*sin (omega t)

and I need to have numerical values to use. I'm trying to use ode45 ,but I keep hitting the error

Error using odearguments (line 93) @(T,Y)[Y(2);-150/1000*Y(1)-150/1000*Y+8*PI*SIN(5*T)/1000] returns a vector of length 3, but the length of initial conditions vector is 2. The vector returned by @(T,Y)[Y(2);-150/1000*Y(1)-150/1000*Y+8*PI*SIN(5*T)/1000] and the initial conditions vector must have the same number of elements.

I feel like i am missing something simple but i haven't figured it out yet. Can anyone show me the process of getting numerical answers, not just a symbolic answer.

My code is:

`f = @(t, y) [y(2); -150/1000*y(1) - 150/1000*y + 8*pi*sin(5*t)/1000];[xa, ya] = ode45(f, [0 50], [ 0 0 ]);`

plot(xa, ya)

## Best Answer