How to solve a system of ODE involving symbolic matrices?

I am trying to solve a second order diff equation and decomposed it into two first order ode

Equations involving 2 variables x(t) and v(t):

           Dx = v ;

           Dv = (-I/A) - (C/A)*x + F/A',

Initial Conditions:

x(0)=0

v(0)=0

A and C are known constants F is a sinusoidal function of time and can be evaluated for all times in the timespan. The difficulty arises because of 'I' being a convolution integral containing unknowns 'v' at every time and is evaluated in form I = const*v(1)+const2*v(2)+....constt*v(t) I is of size symbolic 1x1 containing unknown variable v(1),...v(t) How do I use Matlab to solve the ODE with this symbolic I?

I may have to march forward in time using initial conditions finding x,v at every instant of time and using that as the input for the next instant. How do I do this in MATLAB.

I will be very grateful to help in this direction.