How to solve ODE y'' - 4y' + y - x = 0 using rk4?

edit

I got a way to solve the differential equation. If wrong please correct;)

To solve this differential equation we must first make a variable substitution to reduce the differential equation for a first order and thus create a system of ODEs. variable change:

w1 = y w2 = y w3 = y '

Thus, deriving the above variables, we have a system of equations of the form:

w1 '= w2 w2 '= f (x, w1, w2) = 4 * w2 + x

desolve_system_rk4([w2, 4*w2 + x ],[w1 , w2], ics = [0,1,0],ivar = x, step = 0.1, end_points = 2)