I have an ordinary integro-differential equation of the form $$y''(t)+C_1y(t)+C_2\int_{0}^{t}f(t-\tau)y''(\tau)d\tau +C_3\int_{0}^{t}f(t-\tau)y(\tau)d\tau=0$$ where $C_1$,$C_2$ and $C_3$ are constants.
I know that finding an analytical solution can be hard. How can I solve this numerically?
Any suggestion would be greatly appreciated.