Computer algebra systems, mathematical representation, and the DLMF

Size: px

Start display at page:

Download "Computer algebra systems, mathematical representation, and the DLMF"

Harry Parrish
5 years ago
Views:

1 Computer algebra systems, mathematical representation, and the DLMF Richard Fateman, Bruce Char, Jeremy Johnson University of California, Berkeley Drexel University, Philadelphia National Institute of Standards and Technology DLMF Seminar Series, November 6, 2000

2 Desiderata for the Digital Library of Mathematical Functions l Traditional usage l New modes of interaction Examples l New ambitions Computer Algebra and DLMF 2

3 Non-digital tradition: Finding Out Stuff l Individually owned reference works l Access to libraries references works l Access to colleagues by letter, phone, l Paper and pencil exploration l Numerical experimentation Computer Algebra and DLMF 3

4 Wolfram Research s Special Functions site: 3 versions l Huge posters l Interactive web site/ Mathematica notebooks l Printed form (or the equivalent PDF) Computer Algebra and DLMF 4

5 The posters Computer Algebra and DLMF 5

6 The web site (here, the Arcsin page) Computer Algebra and DLMF 6

7 WRI s Categories/ Some Subcategories primary definition specific values general characteristics series representations generalized power series at various points q-series exponential fourier series dirichlet series asymptotic series other series integral reprsentations on the real axis contour integrals multiple integral representation analytic continuations product representations limit representations continued fractions generating functions group representations differential equations difference equations transformations addition formulas etc operations integral transforms identities representations through more general functions relations with other functions zeros inequalities theorems other information history and applications references Computer Algebra and DLMF 7

8 Click on Series Representations Computer Algebra and DLMF 8

9 This is not very useful l These are blurry pictures of math formulas. l The most plausible next step seems to be to copy them down on paper and check by hand. l There is a possibility of making typos or fresh algebra mistakes. l The notation might be different from what you are using. l Sparse (or no) info on singularities, regions of validity. l To run some numbers through, you need to write a computer program (Fortran? Matlab? C++?) Computer Algebra and DLMF 9

10 Notebook form (I) Input form ArcSin[z] == z^3/6 + z + (3*z^5)/40 + \[Ellipsis] == Sum[(Pochhammer[1/2, k]*z^(2*k + 1))/((2*k + 1)*k!), {k, 0, Infinity}] == z*hypergeometric2f1[1/2, 1/2, 3/2, z^2] /; Abs[z] < 1 One could imagine that a system independent language such as proposed by the OpenMath consortium would replace this language. Note however that agreement on the semantics of \[Ellipsis] would be difficult. Computer Algebra and DLMF 10

11 Notebook form (II) Displayed form (one version) In reality, Mathematica does not look quite as good as this in the interactive mode. Computer Algebra and DLMF 11

12 Notebook form (III) TeX form {Condition}(\arcsin (z) = {\frac{{{\mfunction{z}}^3}}{6}} + z + {\frac{3\,{z^5}}{40}} + \ldots = \Mfunction{\sum}_{k = 0}^{\infty } {\frac{\mfunction{pochhammer}({\frac{1}{2}},k)\, {{\Mfunction{z}}^{2\,k + 1}}}{\left( 2\,k + 1 \right) \,k!}} = \Mfunction{z}\,\Mfunction{Hypergeometric2F1}( {\frac{1}{2}},{\frac{1}{2}},{\frac{3}{2}},{z^2}), \Mfunction{Abs}(z) < 1)) Useful in case you wanted to paste/edit this into another paper, using Mathematica TeX macros. Computer Algebra and DLMF 12

13 Computing Inside the Notebook How good is the 3-term approximation at z= ½? ArcSin[z] == z + z^3/6 + (3*z^5)/ /. z -> 1/2 Í Pi/6 == 2009/ Surprised? N[ Pi/6 == 2009/ ] Í == N[ Pi/6 == 2009/ , 30] Í == N[ Pi/6 == 2009/ , 30] Í == Computer Algebra and DLMF 13

14 Simplification Inside the Notebook In[30] := z* Hypergeometric2F1[1/2, 1/2, 3/2, z^2] Note: this is how Mathematica interactive output looks. This should be the same as ArcSin[z] for z <1. And yes, z/sqrt[z^2] is not the same as 1. Computer Algebra and DLMF 14

15 All commercial computer algebra systems (CAS) have essentially the same notebook paradigm l Macsyma l Maple l Mathematica l Axiom l MuPad l Scientific Word / Maple l Derive Computer Algebra and DLMF 15

16 Advice on coding a reference chapter Computer Algebra and DLMF 16

17 What about legacy knowledge? Can we convert from scanned text? Example from integral table l l In practice, we can do some parsing using OCR if we know about the domains. But in general, we cannot read with understanding. Computer Algebra and DLMF 17

18 What about using LaTeX as source and then converting to OpenMath/ CAS? Generally speaking: not automatically TeX does not distinguish semantically between 1*2*3 and 123. Or between x cos x and xfoox. It has no notion of precedence of operators Gradshteyn and Rhyzik, Table of Integrals and Series (Academic Press) was re-typeset completely in TeX TWICE, because the first version did not reflect semantics. MathML, XML, and OpenMath are inadequate. Computer Algebra and DLMF 18

19 Using OpenMath as original source is pretty much out of the question. Intent is to code: x cos x <OMOBJ> <OMA> <OMS cd = "arith1" name="times"/> <OMV name="x"/> <OMA> <OMS cd="transc1" name="cos"/> <OMV name="x"/> </OMA> </OMA> </OMOBJ> Computer Algebra and DLMF 19

20 Using MathML as original source is pretty much out of the question, too. <math> <msqrt> <mfrac> <mrow><mn>2</mn><mi>π</mi></mrow> <mrow><mi>κ</mi></mrow> </mfrac> <mfenced open="(" close=")"> <mn>1</mn> <mi> </mi> <mi>β</mi> <msup> <mrow><mn>2</mn></mrow> </msup> <mi>/</mi><mn>2</mn></mfenced></msqrt></math> Computer Algebra and DLMF 20

21 What can a CAS do better? Semantics for what makes sense to the CAS is immediate. Presentation for what doesn t make sense to the CAS Advantage: There is an immediate computational ontology Immediate syntactic disambiguation Easy translation into MathML for display Easy translation into OpenMath, if anyone else cares. Computer Algebra and DLMF 21

22 What about using Java Applets? l Pro: an applet provides more intimate interaction with a user. Examples from Math Forum l Con: Only Java-enabled clients can use such applets. Java standardization is problematical High quality numerical software in Java? Symbolic computation in Java? Poor access to underlying computer. Computer Algebra and DLMF 22

23 What about Server Side software? l Pro: arbitrarily powerful could be huge database and super-fast computer always up-to-date controlled by validation team? Can collect / re-distribute new data Computer Algebra and DLMF 23

24 What about Server Side software? l Con: Risk/cost of computation at server Communication requirement Cost Connection reliability Computer Algebra and DLMF 24

25 What about no software? l Pro: You can run DLMF without a computer You can work on a desert island l Con: Anyone with a computer or electricity will be disappointed Computer Algebra and DLMF 25

26 What about only browser software? l Pro: You can run DLMF on an appliance ($300) l Con: Loss in the marketplace of ideas In some respects it will suffer from comparison with software, some of which is free today Students who are used to (say) Mathematica will use other resources, even if less authoritative Computer Algebra and DLMF 26

27 What about numeric-only software? l Pro: You print numeric tables as needed l Con: Symbolic data is endlessly tabulated instead. etc Computer Algebra and DLMF 27

28 What about symbolic software l Now we can consider including algorithms for trig(n/m p) Indefinite integrals (replacing 10-20,000) Summation, Limits Definite integrals (a challenge still) Implicit application of identities, reduction of argument or order by recursion, etc. Generation of any number of terms in series Expansion in Chebyshev or other polynomials Exact arithmetic or bigfloat arithmetic Computer Algebra and DLMF 28

29 A challenge: Include a CAS in DLMF l Free Macsyma (c. 1982) l (buy) Commercial system Macsyma, Maple, Mathematica, Axiom alternatively l Require the user to have a CAS separately (like requiring a Fortran compiler to use GAMS) Computer Algebra and DLMF 29

30 What do we want? What can we attempt? l Exhaustive hierarchical hyperlinks to everything known l Human readable form l Computer usable form + ALGORITHMS l Searchable form /Unique identifiers for formulas l Provenance of information l Annotations from other users l Corrections current or past l Cross-reference hyperlinks to all of mathematics l Applications l A bridge across paper/pencil ÄÅcomputer gap l A tireless, accurate and efficient robot to help us Computer Algebra and DLMF 30

31 Computers do more than arithmetic Many persons who are not conversant with mathematical studies imagine that because the business of [Babbage s Analytical Engine] is to give its results in numerical notation, the nature of its processes must consequently be arithmetical and numerical, rather than algebraical and analytical. This is an error. The engine can arrange and combine its numerical quantities exactly as if they were letters or any other general symbols; and in fact it might bring out its results in algebraic notation, were provisions made accordingly. -- Ada Augusta, Countess of Lovelace, (1844) Computer Algebra and DLMF 31

Extracting Math from PostScript Documents

Extracting Math from PostScript Documents Michael Yang Univ. Calif., Irvine Richard Fateman Univ. Calif, Berkeley ISSAC-2004 1 Why Extract Math from Documents? The current and recent past publications