How to solve Lorenz 96 model using Runge–Kutta method?

Question

I have written the below code to solve Lorenz 96 model using Runge–Kutta method, but the results are not reliable as can be seen in the below figure:

The corr