prime1-20-16

GET THE FULL STORY

photo: “Digits in largest prime by year 2014” by Chiph588Mathematica source code:primeData = {{1952, 157}, {1952, 687}, {1957, 969}, {1961, 1281}, {1961, 1332}, {1963, 2917}, {1963, 3376}, {1971, 6002}, {1978, 6533}, {1979, 6987}, {1979, 13395}, {1982, 25962}, {1983, 39751}, {1985, 65050}, {1989, 65087}, {1992, 227832}, {1994, 258716}, {1996, 378632}, {1996, 420921}, {1997, 895932}, {1998, 909526}, {1999, 2098960}, {2001, 4053946}, {2003, 6320430}, {2004, 7235733}, {2005, 7816230}, {2005, 9152052}, {2006, 9808358}, {2008, 12978189}, {2013, 17425170}};fitEq = Exp[a x + b] /. FindFit[{#1, Log[#2]} & @@@ primeData, a x + b, {a, b}, x];primes = ListLogPlot[primeData, PlotRange -> {{1950, 2014}, {90, Automatic}}, Joined -> True, Frame -> True, FrameStyle -> GrayLevel[.3],GridLines -> {Range[1960, 2010, 10], 10^Range[3, 7]}, GridLinesStyle -> Directive[Gray, Dashed, Opacity[.25]], ImageSize -> 1000,FrameTicksStyle -> Directive[FontFamily -> “Times”, FontSize -> 22], Mesh -> All, MeshStyle -> Directive[ColorData[1][1], PointSize[Medium]]];bestfit = LogPlot[fitEq, {x, 1950, 2014}, PlotStyle -> Directive[Red, Dashed, AbsoluteDashing[5], Opacity[.5]]];res = Show[primes, bestfit]. Licensed under CC BY-SA 3.0 via Commons.