I have the following equations (which is a system of three ODEs). They should be solved nomerically (for example with 'ode45'). The problem is that Delta in the following equations is not constant and varies during the time. For example, if we consider 64 time elements in the time axis, there is a vector of 64 components for Delta. Each component is constant but they have different vaues. Is it possible to take such a constant (in a vector form) into consideration in 'ode45' or any othe ode solver?