Issue on Runge Kutta Fehlberg algorithm

Question

I have wrote a code for Runge-Kutta 4th order, which works perfectly fine for a system of differential equations:

import numpy as np
import matplotlib.pyplot