Add time dependent terms in a differential equation

New to Julia
5

You didn’t share data so I just set this up with random data.

using DataInterpolations, OrdinaryDiffEq
feedTime = Array(0:0.2:14.0)
feedData = [rand(1) for i in 1:length(feedTime)]
rate = CubicSpline(feedData, feedTime)

function Volume(du, u, p, t)
    du .= u + rate(t)
end

u0 = [1.0]
tspan = (0.0, 14.0)
prob = ODEProblem(Volume, u0, tspan, rate)
sol = solve(prob, Tsit5(), tstops= feedTime)