Hi, can you please help me on how to solve the Lorenz equation using Runge-Kutta-Gill method in Matlab.

Lorenz Equations :

dx/dt =σ (y-x)

dy/dt = Rx-y-xz

dz/dt=xy-bz

where

σ = 10, R = 28, b = 8/3;

t, time span; t = 0 to t = 25; time increment Δt = 0.001;

initial values (x(0), y(0), z(0)) = (10, 10, 10)