Issue understanding scipy.integrate.RK45 requirements

I am attempting to solve a system of first order differential equations with scipy.integrate.RK45(). I have scripted out the model function that I am hoping to plot (displacement vs time), but RK45() requires that this function takes 2 arguments, namely 't' and 'y', where 't' is a scalar and 'y' is an array in my case. This is better described below:

https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.RK45.html

My script is as follows:

import numpy as np
from scipy.integrate import RK45, RK23

# Model the calling function [X_dot]:
#------------------------------------------------------------------
def model(t, y):
    
    # Define the mass, stiffness and damping [m, c, k]:
    m = 2
    c = 10
    k = 1500

    # Define the piecewise forcing function [F]:
    if (t >= 0 and t < 0.1):
        F = 200 * t
    if (t >= 0.1 and t < 0.25):
        F = 20
    else:
        F = 0
    
    E_matrix = [[0, 1], [(-k / m), (-c / m)]]
    Q_matrix = [0, F / m]

    return E_matrix*X + Q_matrix

# Define the initial conditions and integration boundaries:
#------------------------------------------------------------------
time_step = 0.01
t_upper = 0.5
t_lower = 0

initial_conditions = [0, 0] # [displacement(t=0) = 0, velocity(t=0) = 0]

points_to_plot = RK45(fun=model(t, y), t0=t_lower, y0=initial_conditions, t_bound=t_upper, vectorized=True)

A picture of the system I'd like to solve is seen below:

I have found very few examples of this being used as most solutions use odeint().

What are these two arguments (t, y) and how would I effectively incorporate them into my function? Any help is appreciated greatly.