matlab code . How can i find Lyapunov Exponent four dimension matlab code ?

How can i find Lyapunov Exponent four dimension matlab code ?

Lyapunov Exponent of Continuous chaotic system of four dimension. Also how to plot its spectrum. (matlab code) or Mathematica code example

dw/dt= 23*(x -w)

dx/dt= 9*w - w*y

dy/dt= w*x+x*z

dz/dt=x*y- y

for initial conditions (1, 1, 1, 1),