摘要

此graph图片可使用矢量图形重新创建为SVG文件。这具有很多好处；更多信息请参见共享资源:待清理媒体。如果存在此图片的SVG格式，请将其上传，然后将此模板替换为{{vector version available|新图片名称}}。 建议您将SVG命名为“BenfordBroad.svg”，这样在使用Vector version available（或Vva）模板时就不需要填写新图片名称参数。

许可协议

我，本作品著作权人，释出本作品至公有领域。这适用于全世界。 在一些国家这可能不合法；如果是这样的话，那么： 我无条件地授予任何人以任何目的使用本作品的权利，除非这些条件是法律规定所必需的。

(* 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};

ff[x_] := 
  x Exp[-.1 x] .1 + x Exp[-.03 x] .03 +  x Exp[-.01 x] .01 + 
   x Exp[-.003 x] .003 + x Exp[-.001 x] .001 + x Exp[-.0003 x] .0003;

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);

plotA = LogLinearPlot[{ff[x], If[firstdig1[x], ff[x], 0], If[firstdig8[x], ff[x], 0]}, {x, .99, 10100}, 
  Ticks -> {LabelsList, {}}, PlotRange -> {0, 1}, 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["plotA.gif", plotA, ImageSize -> 500]

原始上传日志