I need a Python code on how to plot bifurcation diagrams of Lorenz system with varying fractional order.

The one I’ve been able to come up with is this image and I am not sure it is correct

def getAndPlotBif():

  Z = []

  ordss = genOrds(0.8, 1.0, 0.001)

  ordz = []

  xs0, ys0, zs0 = x0

  for ord in ordss:

    x, y, z = ABM(ord, 3000, 0.004)

    for i in range(len(z)):

      Z.append(x[i])

      ordz.append(ord)

    xs0, ys0, zs0 = x[-1], y[-1], z[-1]

  Z = np.array(Z)

  ordz = np.array(ordz)

  return ordz, Z

This is the plotter code here:

def bif_plot(ords, Z, title="Bifurcation diagram of Lorenz system with different derivative order", x='Derivative order'):

  pl = plt

  pl.figure(figsize=(10, 6))

  pl.scatter(ords, Z, color='blue', s=0.5, alpha=0.2)

  pl.title(title)

  pl.xlabel(x)

  pl.ylabel("Z values")

  pl.show()

This is the result image