A  bifurcation diagram for the  Logistic map:
xn+1=rxn(1−xn){\displaystyle x_{n+1}=rx_{n}(1-x_{n})\,}
The horizontal axis is the r parameter, the vertical axis is the x variable. A starting value of x=0.25 was used, and the map was iterated 1000 times in order to stabilize the values of x. 1,000,000 x -values were then calculated for each value of r and for each x value, the corresponding (x,r) pixel in the image was incremented by one. All values in a column (corresponding to a particular value of r) were then multiplied by the number of non-zero pixels in that column, in order to even out the intensities. See File:Logistic Map Bifurcation Diagram, Matplotlib, zoomed.png for a vectorized and therefore much faster version of this Python script.
Own work
Morn