This site is supported by donations to The OEIS Foundation.

User:Enrique Pérez Herrero/OEIS Package

From OeisWiki

Jump to: navigation, search

Contents

Mathematica OEIS Package

When you go on computing: if you are doing repetitive and boring things, is because you are not working in a proper way. So I decided to write some code to help myself with the worst part of my projects. And as this code began to grow and became more complex, thought it would be good idea to gather all together in the same Mathematica Package and to share it on the internet.

This code includes functions that can make your life easier if you a using the information from OEIS in your archives: for citation in documents or blogs, for programming sequences, or for many other purposes.


Package Description:

  • Version:
    • V 1.0, 26 August 2010, Testing...
    • V 2.0, 28 November 2010, Old links are no longer valid.


  • Changes:
    • 29/Aug/2010 - Added OEISbFile function and exporting to HTML
    • 28/Nov/2010 - Code rebuilt due to changes in OEIS, old links are no longer valid.
    • 04/Jun/2011 - Fixed Unicode LF, instead of Window's CRLF for bFiles.
  • Functions:
    • OEISTotalNumberOfSequences
    • OEISValidateIDQ
    • OEISImport
    • OEISURL
    • OEISFunction
    • OEISExport
    • OEISbFile

Download and Install OEIS.m

Download OEIS.m

  • Loading Package into a Mathematica notebook:
<<OEIS`

PROGRAMMING EXAMPLES

IMPORTING DATA FROM OEIS:

Importing Sequence Description:

In[1]:= OEISImport["A014321","Description"]
Out[1]= The next new gap between successive odd primes (divided by 2).

Importing Sequence Data:

 In[2]:= OEISImport["A080321","Data"]
Out[2]= {{0,0},{1,4},{2,428},{3,58747},{4,1758473969},{5,340760636747},{6,14410746502617904},{7,673020037180959760871},
{8,1736622037158911219053381734718},{9,400539731185759177060380684246324}}

Importing Sequence creation date:

In[3]:= OEISImport["A117321","Date"]
Out[3]= {Mar,07,2006}

Importing Sequence Author:

In[4]:= OEISImport["A014321","Author"]
Out[4]= {Hynek Mlcousek}

Importing Image:

In[5]:= OEISImport["A030337","Image"]
Out[5]=

Importing bFile:

In[6]:=ListPlot[OEISImport["A094358","bFile"],Joined->True]
Out[6]= (Plot Image)

Importing Range of Available data:

In[7]:= OEISImport["A094358","MinData"]
Out[7]= 1
In[8]:= OEISImport["A094358","MaxData"]
Out[8]= 1314


EXPORTING DATA FROM OEIS:

  • Exporting Formats:
OEIS.m exporting formats
File Extension Exported data:
*.gp Code to add help for PARI/GP functions
*.m Code to add help for Mathematica
*.txt, *.xls, *.csv, *.tsv, *.dat Data files with sequence values.
*.jpg, *.jpeg, *.jp2, *.j2k, *.bmp, *.pgn, *.gif, *.tiff, *.tif Image plot
*.bib BibTeX citation to OEIS Sequences
*.htm, *.html HTML citation to OEIS Sequences

Exporting BibTeX citation:

The input for the function can be a String or a List of Strings:

In[1]:=OEISExport[{"A026588", "A111084"}, "example.bib"]
Out[1]=example.bib

Content of file: example.bib

@MISC{oeisA026588,
AUTHOR={},
TITLE={The {O}n-{L}ine {E}ncyclopedia of {I}nteger {S}equences},
HOWPUBLISHED={\href{http://oeis.org/A026588}{A026588}},
MONTH={},
YEAR={},
NOTE={Duplicate of A026572.}
}

@MISC{oeisA111084,
AUTHOR={N. J. A. Sloane},
TITLE={The {O}n-{L}ine {E}ncyclopedia of {I}nteger {S}equences},
HOWPUBLISHED={\href{http://oeis.org/A111084}{A111084}},
MONTH={Nov},
YEAR={2005},
NOTE={Where records occur in A111267.}
}

Exporting Help Code for PARI/GP:

In[2]:= OEISExport[{"A026588","A111084","A131473"},"code_for_PARIGP.gp"]
Out[2]= code_for_PARIGP.gp
  • Content of file: code_for_PARIGP.gp
addhelp(A026588,"A026588: Duplicate of A026572.");
addhelp(A111084,"A111084: Where records occur in A111267.");
addhelp(A131473,"A131473: a(n)=n^6-n.");

Exporting Help Code for Mathematica:

In[3]:= OEISExport[{"A026588","A111084","A131473"},"code_for_Mathematica.m"]
Out[3]= code_for_Mathematica.m
  • Content of file: code_for_Mathematica.m
A026588::usage="A026588[n]: Duplicate of A026572.";
A111084::usage="A111084[n]: Where records occur in A111267.";
A131473::usage="A131473[n]: a(n)=n^6-n.";

Exporting all available data

  • Example with EXCEL:
In[4]:= OEISExport[{"A026588","A111084","A131473"},"data.xls"]
Out[4]= {A026588_data.xls,A111084_data.xls,A131473_data.xls}

Exporting citation to HTML code

In[5]:=OEISExport[{"A040414", "A050414", "A002323", "A026588"}, "test.html"]
Out[5]= test.html
  • Content of file: test.html
<a name="oeisA040414">[1]-</a> N. J. A. Sloane, The On-Line Encyclopedia of Integer Sequences. 
<a href=http://oeis.org/A040414">A040414</a>: Continued fraction for sqrt(435).
<a name="oeisA050414">[2]-</a> Jud McCranie, The On-Line Encyclopedia of Integer Sequences. <a href=http://oeis.org/A050414">A050414</a>: Numbers n such that 2^n - 3 is prime.
<a name="oeisA002323">[3]-</a> N. J. A. Sloane, The On-Line Encyclopedia of Integer Sequences. <a href=http://oeis.org/A002323">A002323</a>: Fermat remainders.
<a name="oeisA026588">[4]-</a> The On-Line Encyclopedia of Integer Sequences. <a href=http://oeis.org/A026588">A026588</a>: Duplicate of A026572.


Exporting citation to Wiki code

In[6]:=OEISExport[{"A026518", "A111284"}, "example.wiki"]
Out[6]:=example.wiki
  • Content of file: example.wiki
* A026518: a(n) = length of n-th run of identical symbols in A026517.
* A111284: Number of permutations avoiding the patterns {2143,2341,2413,2431,3142,3241,3412,3421,4123,4213,4231,4321,4132,4312};
number of strong sorting class based on 2143.

FUNCTIONS CREATED AUTOMATICALLY: Function: OEISFunction[ID,Opts]:

The options Opts, are all of them True by default, and they can be modified as:

To skip the addition of help to the function:

OEISFunction[ID,AddHelp->False]

To skip the data from the bFile:

OEISFunction[ID,bFile->False]

To do not print on the screen the information outputed:

OEISFunction[ID,Output->False]
  • Example of use:
In[5]:= OEISFunction["A094358"]
A094358(n): Squarefree products of factors of Fermat numbers (A023394).
Link-1: http://oeis.org/A094358
Link-2: http://oeis.org/b094358.txt
Help Added: True
Loaded data from bFile: True
Loaded data from Sequence: False
Data available: Table of n, A094358(n) for n=1..1314

From now on, we can use the function with all data already loaded on it:

In[6]:= A094358[1314] 
Out[6]= 2917004348489729

GENERATING b-Files:

Function: OEISbFile[ID,Vmax,filename]


  • Example:

We can define the function first, but if we do not define it, the code creates the function and preloads all available data from OEIS:

In[5]:=A040414[0]:=20;
       A040414[n_]:=Part[{1,5,1,40},1+Mod[n-1,4}};

In[6]:=OEISbFile["A040414", 200]
Out[6]=b040414.txt

Links:

Personal tools