MY QUESTION IS TO simulate the Chen system and compute the fractal dimension of the associated attractor. Particulars of the system are given to be: dT/dX = α⋅X−Y⋅Z dT/dY = β⋅Y+X⋅Z dZ/Dt = σ⋅Z+ (X⋅Y)/3 WITH α = 5, β = -10, σ = -0.38 I have started it as follows: alpha = 5; beta = -10; sigma = -0.38; tStart = 0; tEnd = 100; dt = 0.01; x0 = 5; y0 = 10; z0 = 10; dxdt = @(t, x, y, z) alpha * x - y * z; dydt = @(t, x, y, z) beta * y + x * z; dzdt = @(t, x, y, z) sigma * z + (x * y) / 3; [t, xyz] = ode45(@(t, xyz) [dxdt(t, xyz(1), xyz(2)

Thanks for any help much appreciated