This file is from Wikimedia Commons and may be used by other projects.
The description on its file description page there is shown below.
Summary
DescriptionBenfordNarrow.gif
English: A narrow probability distribution plotted on a log scale, to demonstrate a case where Benford's law does not hold. This is to be contrasted with BenfordBroad.gif.
This graph image could be re-created using vector graphics as an SVG file. This has several advantages; see Commons:Media for cleanup for more information. If an SVG form of this image is available, please upload it and afterwards replace this template with {{vector version available|new image name}}.
It is recommended to name the SVG file “BenfordNarrow.svg”—then the template Vector version available (or Vva) does not need the new image name parameter.
Licensing
Public domainPublic domainfalsefalse
I, the copyright holder of this work, release this work into the public domain. This applies worldwide. In some countries this may not be legally possible; if so: I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law.
(* Source code written in Mathematica 6.0, by Steve Byrnes, 2008. I release this code into the public domain. *)
LabelsList = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 70, 100, 200, 300, 500, 1000, 2000, 5000, 10000};
LineList = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600,
700, 800, 900, 1000, 2000, 3000, 4000, 5000, 6000, 7000, 8000, 9000, 10000};
gg[xx_] := Exp[-Abs[(xx - 60)/25]^1.8] xx
firstdig1[xx_] :=
(1 < xx && xx < 2) || (10 < xx && xx < 20) || (100 < xx && xx < 200) || (1000 < xx && xx < 2000);
firstdig8[xx_] :=
(8 < xx && xx < 9) || (80 < xx && xx < 90) || (800 < xx && xx < 900) || (8000 < xx && xx < 9000);
plotB = LogLinearPlot[{gg[x], If[firstdig1[x], gg[x], 0], If[firstdig8[x], gg[x], 0]}, {x, .995, 10100},
Ticks -> {LabelsList, {}}, PlotRange -> All, TicksStyle -> Directive[Black, 10],
GridLines -> {LineList, {}}, ImageSize -> 401, AspectRatio -> .3,
Filling -> {2 -> Bottom, 3 -> Bottom}, PlotPoints -> 1000, PlotStyle -> {Black, Red, Blue},
FillingStyle -> {2 -> Red, 3 -> Blue}]
Export["plotB.gif", plotB, ImageSize -> 500]