Origin of Irrational Numbers and Their Approximations

: In this article a sincere effort has been made to address the origin of the incommensu-rability/irrationality of numbers. It is folklore that the starting point was several unsuccessful geometric attempts to compute the exact values of √ 2 and π . Ancient records substantiate that more than 5000 years back Vedic Ascetics were successful in approximating these numbers in terms of rational numbers and used these approximations for ritual sacriﬁces, they also indicated clearly that these numbers are incommensurable. Since then research continues for the known as well as unknown/expected irrational numbers, and their computation to trillions of decimal places. For the advancement of this broad mathematical ﬁeld we shall chronologically show that each continent of the world has contributed. We genuinely hope students and teachers of mathematics will also be beneﬁted with this article.


Introduction
Almost from the last 2500 years philosophers have been unsuccessful in providing satisfactory answer to the question "What is a Number"? The numbers 1, 2, 3, 4, · · · , have been called as natural numbers or positive integers because it is generally perceived that they have in some philosophical sense a natural existence independent of man. We will never know if there existed a genius who invented or introduced these natural numbers, but it is generally accepted that these numbers came down to us, ready-made, from an antiquity most of whose aspects are preserved in folklore rather than in historical documents. For primitive man and children natural number sense is an inherent ability. There are several recorded incidences of birds, animals, insects, and aquatic creatures who show through their behavior a certain natural number sense. While natural numbers are primarily used for counting finite collections of objects, there is hardly any aspect of our life in which natural numbers do not play a significant-though generally hidden-part. In fact, natural numbers are building blocks of all sciences and technologies. Number Theory which mainly deals with properties and relationships of natural numbers for their own sake has been classified as pure mathematics. Since antiquity, number theory has captivated the best minds of every era. An important feature of number theory is that challenging problems can be formulated in very simple terms; however, hidden within their simplicity is complexity. Some of these problems have been instrumental in the development of large parts of mathematics. Amateurs and professionals are on an almost equal footing in this field. The set of all natural numbers is denoted as N .
A positive rational number is defined as the exact ratio/fraction/quotient of two positive integers p/q, where q = 0. It is very likely that the notion of rational numbers also dates to prehistoric times. Around 4000 BC, rational numbers were used to measure various quantities, such as length, weights, and time in the Indus river valley (which was home to more than five million people). Thus, then rational numbers were sufficient for all practical measuring purposes. The Babylonians used elementary arithmetic operations for rational numbers as early as 2000 BC. We also find ancient Egyptians texts describing how to convert general fractions into their special notation. Classical Greek and Indian mathematicians made studies of the theory of rational numbers, as part of the general study of number theory, see Euclid's Elements (300 BC) and Sthananga Sutra (around 3rd century). The set of all positive rational numbers is denoted as Q + .
Throughout the ancient history negative solutions of linear and quadratic equations have been called as absurd solutions. First systematic use of negative numbers in mathematics for finding the solutions of determinate and indeterminate systems of linear equations of higher order with both positive and negative numbers appeared in Chinese work much before Han Dynasty (202 BC-220 AD). In appreciation, the historian Jean-Claude Martzloff (1943 theorized that the importance of duality in Chinese natural philosophy made it easier for the Chinese to accept the idea of negative numbers. Brahmagupta (born 30 BC, India) in his treatise Brahmasphutasiddhanta treated negative numbers in the sense of 'fortunes' and 'debts', he also set rules for dealing with negative numbers. Most importantly, he treated zero as a number in its own right, and attempted to define division by zero. For a long history of zero, its role in life, and mathematics, see Sen and Agarwal [1]. Unfortunately, in Britain pessimistic attitude towards negative numbers continued till 18th century, in fact, William Frend (1757-1841, England) took the view that negative numbers did not exist, whereas his contemporary Francis Maseres (1731-1824, England) in 1759 wrote that negative numbers "darken the very whole doctrines of the equations and make dark of the things which are in their nature excessively obvious and simple". He came to the conclusion that negative numbers were nonsensical. However, in the 19th century negative numbers received their relevance logically across the world. The set of integers including positive, negative, and zero is denoted as Z, and the set of all rational numbers is represented by Q.
Numbers which cannot be expressed as ratios of two integers are called incommensurable or irrational (not logical or reasonable). The earliest known use of irrational numbers is in the Indian Sulbasutras. For ritual sacrifices there was a requirement to construct a square fire altar twice the area of a given square altar, which lead to find the value of √ 2 (in the literature it has been named as Pythagoras number). Indian Brahmins also needed the value of π (the ratio of the circumference to the diameter of a circle). They were successful in finding reasonable rational approximations of these numbers, keeping in mind the success of ritual sacrifices depending on very precise mathematical accuracy. In Sulbasutras there is also a discussion that these numbers cannot be computed exactly. Thus, the concept of irrationality was implicitly accepted by Indian Brahmins. We also find approximations of √ 2 in Babylonians tablets using sexagesimal fractions. In Greek geometry, two magnitudes a and b of the same kind were called commensurable if there is another magnitude c of the same kind such that both are multiples of c, that is, there are numbers p and q such that pc = a and qc = b. If the two magnitudes are not commensurable, then they are called incommensurable. While decimal fractions and decimal place value notation, a gift from India to whole world, has a long history, decimal fraction approximations of √ 2 and π appeared during 200-875 AD, in the Jain School of Mathematics (India). In terms of decimal expansions unlike a rational number, an irrational number never repeats or terminates. In fact, it is only the decimal expansion which immediately shows the difference between rational and irrational numbers. Irrational numbers have also been defined in several other ways, e.g., an irrational number has nonterminating continued fraction whereas a rational number has a periodic or repeating expansion, and an irrational number is the limiting point of some set of rational numbers as well as some other set of irrational numbers. In what follows, we will correct the speculations that incommensurability of √ 2 was proved by Pythagoras himself (and for all nonsquare integers by Theodorus), by reveling that the first (fully geometric) proof appeared in the Meno (Socratic dialogue by Plato). Here we will see an infinite process arise in an attempt to understand irrationals.
Since then over the period of 2400 years many different proofs of the irrationality of √ 2 have been given, we will demonstrate few of these, and furnish several algorithms to find its rational approximations. The proof of the irrationality of π had to wait almost two millennia, it was proved only in 1768 by Johann Heinrich Lambert (1728-1777, Switzerland). In 1683 another important number e was introduced by Jacob Bernoulli (1654-1705, Switzerland), whose irrationality was proved by Leonhard Euler (1707-1783, Switzerland) in 1748. Thus, the numbers √ 2, π, and e have infinite number of decimal places. Since the invention of computer technology, these numbers have been approximated to trillions of decimal places, we shall report these accomplishments. It is to be noted that such extensive calculations besides human desire to break records, have been used to test supercomputers and high-precision multiplication algorithms, the occurrence of the next digit seems to be random, and the statistical distribution expected to be uniform. We list here first 100 digits of these numbers, which are more than sufficient (in fact, not even first twenty) for each and every real world problem. The set of all irrational numbers is denoted as Q . The union of the sets of all rational and irrational numbers make up the set of real numbers denoted as IR. Thus, this large set contains all decimal representations of numbers terminating, repeating, nonterminating, and nonrepeating.
Euler in his work noted that e is of a different kind of irrational number, which lead to transcendental numbers (not the roots of nonzero polynomials with rational coefficients). While the existence of transcendental numbers have been proved to be uncountable, only for very few numbers their transcendence (one by one) has been established. As it stands, even to prove irrationality of a number no general method exists, proving transcendence (or otherwise) of a number is considered as life's great achievement. We shall provide a detailed account of this field.
From the 9th century, Arabic mathematicians started treating irrational numbers as algebraic objects, and initiated the idea of merging the concepts of number (algebra) and magnitude (geometry) into a more general idea of real numbers. Specially, in the 10th century they provided a geometric interpretation of rational numbers, on a horizontal straight line. This work was completed for all real numbers only in the 19th century, which is now known as Dedekind-Cantor axiom.

Sulbasutras
The meaning of the word sulv is to measure, and geometry in ancient India came to be known by the name sulba or sulva. The Sulbasutras are the appendices to four Vedas (means wisdom, knowledge or vision): Rigveda, Samaveda, Yajurveda, and Atharvaveda. Sulbasutras were codified by Krishna Dwaipayana or Sage Veda Vyasa (born 3374 BC) along with his disciples Jaimani, Paila, Sumanthu, and Vaisampayana. Only seven Sulbasutras are extant, named for the sages who wrote them: Apastamba, Baudhayana (born 3200 BC), Katyayana, Manava, Maitrayana, Varaha, and Vidhula. The four major Sulbasutras, which are mathematically the most significant, are those composed by Baudhayana, Manava, Apastamba, and Katyayana. These Sulbasutras contain a large number of geometric constructions for squares, rectangles, parallelograms and trapezia; the problem of solving quadratic equations of the form ax 2 + bx + c = 0; several examples of arithmetic and geometric progressions; a method for dividing a segment into seven equal parts; solutions of first degree indeterminate equations; and (without any proofs) remarkable approximations of √ 2 (the sign √ was introduced by Christoff Rudolff 1499-1545, Austria) and π (the Computation 2021, 9, 29 4 of 49 ratio of the circumference of a circle to its diameter), the Greek symbol π was used first by the Welshman William Jones (1675-1749, UK, in 1706). In three Sulbasutras Baudhayana, Apastamba, and Katyayana for the approximation of √ 2 the recipe is "increase the measure by its third and this third by its own fourth less the thirty-fourth part of that fourth. This is the value with a special quantity in excess". If we take 1 unit as the dimension of the side of a square, then this in modern terms can be written as (1) and, similarly, if we take the radius of the circle as 1 unit, then the approximation formula for π is These approximations were used for the construction of altars, particularly, √ 2 in an attempt to construct a square altar twice the area of a given square altar. For a successful ritual sacrifice, the altar had to conform to very precise measurements, so mathematical accuracy was seen to be of the utmost importance. Bibhutibhushan Datta (1888-1958 in his most trusted treatise [2] on Sulbas on page 27 writes "The reference to the sacrificial altars and their construction is found as early as the Rigveda (before 3000 BC). ... It seems that the problem of the squaring of the circle and the theorem of the square of the hypotenuse are as old in India as the time of Rigveda. They might be older still". Approximation (1) gives √ 2 1.414215686, which is correct to five decimal places. Perhaps the approximation (1) was used in π 18(3 − 2 √ 2), to obtain π 105/34 3.088235294. George Gheverghese Joseph (born 1928, India) in his book [3] mentions about his correspondence with Takao Hayashi (born 1949, Japan) who pointed out that the approximation of √ 2 could also be used for constructing a right-angled triangle and a square. To show (1), Datta on pages 193,194, and subsequently by several others, e.g., Joseph on pages 235,236 have provided the following reasoning which is in line with Sulbasutra's geometry. Consider two squares, ABCD and PQRS, each of 1 unit as the side of a square (see Figure 1). Divide PQRS into three equal rectangular strips, of which the first two are marked as 1 and 2. The third strip is subdivided into three squares, of which the first is marked as 3. The remaining two squares are each divided into four equal strips marked as 4 to 11. These eleven areas are added to the square ABCD as shown in Figure 1, to obtain a larger square less a small square at the corner F. The side of the augmented square AEFG is The area of the shaded square is [1/(3 · 4)] 2 , so that the area of the augmented square AEFG is greater than the sum of the areas of the original squares, ABCD and PQRS, by Computation 2021, 9, 29 5 of 49 Now to make the area of the square AEFG approximately equal to the sum of the areas of the original squares ABCD and PQRS, imagine cutting off two very narrow strips, of width x, from the square AEFG, one from the left side and one from the bottom. Then (3) Simplifying the above expression and ignoring x 2 , an insignificantly small quantity, gives The diagonal of each of the original squares is √ 2, which can be approximated by the side of the new square as just calculated, i.e., (1).
A commentator on the Sulbasutras, Rama (perhaps Rama Chandra) Vajapeyi, who lived in the middle of the fifteenth century AD in India, gave an improved approximation to (1) by adding two further terms to the equation, i.e., which gives 1.414213502 a value correct to seven decimal places. In Sulbasutras we also find approximation of √ 3, which can be written as Approximation (5) gives √ 3 1.732051282, which is correct to five decimal places. In (Datta [2], pp. 194-195), a geometric construction similar to that of (1) for (5) is also given. A simple algebraic method to get (5) is to take 5/3 as an approximation of √ 3 and put √ 3 = (5/3 + x), where x is unknown. Now square both sides of this expression, neglect x 2 , and solve the resulting linear equation for x, to get x = 1/3 · 5, thus the new approximation of √ 3 is 26/15. Repeating this procedure once more, we find x = −1/3 · 5 · 52 and the new approximation of √ 3 as 1351/780. For (1) several other descriptions have been proposed, e.g., Radha Charan Gupta (born 1935, India), in [4] uses linear interpolation to obtain the first two terms of (1), he then corrects the two terms so obtaining the third term, then correcting the three terms obtaining the fourth term. In The first indentity gives √ 2 7/5 = 1.4, whereas the second gives √ 2 17/12 = 1.41666666 · · · .
For an excellent detailed discussion of √ 2 up to 2006, see the book of Flannery [5]. Bonnell and Nemiroff on the Website https://apod.nasa.gov/htmltest/gifcity/sqrt2.1mil (accessed on 4 March 2021) have posted one million digits of √ 2, and in 2009 five million digits, see Bonnell and Nemirof [6]. Other records are by Yasumasa Kanada (1949-2020, Japan) in 1997 to 137,438,953,444 decimal places; Shigeru Kondo (born 1959, Japan) in 2010 to one trillion decimal places; Alexander Yee in 2012 to two trillion; Ron Watkins in April, 2016 to five trillion, and in June 2016 to ten trillion.
In Sulvasutras, the priests gave the following procedure for finding a circle whose area was equal to a given square. In the square ABCD, let M be the intersection of the diagonals ( Figure 2). Draw the circle with M as center and MA as radius, let ME be the radius of the is the same as (2). r the converse problem, that of squaring the circle, we are given the following rule: If you wish to le into a square, divide the diameter into 8 parts, and again one of these 8 parts into 29 parts; of t rts remove 28, and moreover, the sixth part (of the one left) less the eighth part (of the sixth p eaning is: side of the required square is the diameter of given circle. It gives the value of π = 3.088326491. ll the Sulbasutras contain a method to square the circle. It is an approximate method based ructing a square of side 13/15 times the diameter of the given circle as in the Figure 3. This corresp king the value of π as π = 4 × (13/15) 2 = 676/225 = 3.00444. For the converse problem, that of squaring the circle, we are given the following rule: If you wish to turn a circle into a square, divide the diameter into 8 parts, and again one of these 8 parts into 29 parts; of these 29 parts remove 28, and moreover, the sixth part (of the one left) less the eighth part (of the sixth part). The meaning is: side of the required square is 7 8 times the diameter of given circle. It gives the value of π = 3.088326491. All the Sulbasutras contain a method to square the circle. It is an approximate method based on constructing a square of side 13/15 times the diameter of the given circle as in Figure 3. This corresponds to taking the value of π as π = 4 × (13/15) 2 = 676/225 = 3.00444.  rth noting that many different values of π appear in the Sulbasutras, even severa same text. This is not surprising that whenever an approximate construction is g s implied. The authors thought in terms of approximate constructions, not in term s with π but only having an approximate value for it. For example, in Baudhayana t values of π are given as 676/225, 900/289, and 1156/361. In other Sulbasutras 3.004, 3.029, 3.047, 3.088, 3.1141, 3.16049 and 3.2022 can all be found. Particula basutra, see Gupta [23], the value of π ≃ 25/8 ≃ 3.125, also see interesting work of Ka 9]. For an extensive history of π (calculating up to ten trillion decimal places) till the l, et. al. [2]. In 2019, a Google cloud developer Emma Haruka Iwao from Japan set a alculating π to 31.4 trillion decimal places. She used the same software as her succe trillion, 2016) but had the advantage owing to her access to Google servers. The calcu -five cloud-based computers and a hundred and twenty one days to complete. On 29t thy Mullican of USA has broken all previous records by calculating π to 50 trillion di George Thibaut (1848-1914) translated a large portion of the Sulvasutras, which sh priests possessed significant mathematical knowledge. Thibaut was a Sanskrit schol jective was to make the mathematical knowledge of the Vedic Indians available to t firmly believed that Hindus had knowledge of irrationality, in particular, of √ 2. there is a discussion of the irrationality of π. According to Datta [14, p 195  It is worth noting that many different values of π appear in the Sulbasutras, even several different ones in the same text. This is not surprising that whenever an approximate construction is given some value of π is implied. The authors thought in terms of approximate constructions, not in terms of exact constructions with π but only having an approximate value for it. For example, in Baudhayana Sulbasutra the different values of π are given as 676/225, 900/289, and 1156/361. In other Sulbasutras the values 2.99, 3.00, 3.004, 3.029, 3.047, 3.088, 3.1141, 3.16049 and 3.2022 can all be found. Particularly, in the Mayana Sulbasutra, see Gupta [7], the value of π 25/8 3.125, also see interesting work of Kak [8] and Kulkarni [9]. For an extensive history of π (calculating up to ten trillion decimal places) till the year 2013, see Agarwal, et al. [10]. In 2019, a Google cloud developer Emma Haruka Iwao from Japan set a new world record for calculating π to 31.4 trillion decimal places. She used the same software as her successor (Peter Trueb-22.4 trillion, 2016) but had the advantage owing to her access to Google servers. The calculation took over twenty-five cloud-based computers and a hundred and twenty one days to complete. On 29 January 2020, Timothy Mullican of USA has broken all previous records by calculating π to 50 trillion digits.
In 1875, George Thibaut (1848-1914) translated a large portion of the Sulvasutras, which showed that the Indian priests possessed significant mathematical knowledge. Thibaut was a Sanskrit scholar and his principal objective was to make the mathematical knowledge of the Vedic Indians available to the learned world. He firmly believed that Hindus had knowledge of irrationality, in particular, of √ 2. In fact, in Apastamba there is a discussion of the irrationality of π. According to Datta

Aryabhata's Method for Extracting Square and Cube Roots
The legacy of this Indian genius (born 2765 BC) continues to baffle mathematicians and astronomers, for details of his astonishing contributions, see Agarwal and Sen [11] and Keller [12]. Although, Aryabhata does not provide details to find square and cube roots, it has been concluded that his method is based on decimal place-value system, and the equalities (ab) 2 = (10a b. An important feature of his method is that it finds each Computation 2021, 9, 29 8 of 49 digit of the root successively, from left to right. His method is still taught in schools. We shall summarize his method in simplified terms through the following examples. To find the square root of 625, we group it in two's from right to left as 6 25. Now search largest possible integer a such that a 2 ≤ 6, which is obviously a = 2. This will be the first digit of the required square root. The next step is to find 6 − a 2 = 6 − 2 2 = 2, and with this adjoin 25, i.e., 225. Now find largest possible integer b such that [2(10a) + b]b = (40 + b)b ≤ 225, which is obviously b = 5. This will be the next digit of the required square root. Since 225 − (45)5 = 0, it follows that √ 625 = 25. To find the square root of 474721, we group it in two's from right to left as 47 47 21. Search largest possible integer a such that a 2 ≤ 47, which is 6. Now, we find 47 − a 2 = 47 − 6 2 = 11, and with this adjoin 47, i.e., 1147 and find largest possible integer b such that [2(10a) + b]b = (120 + b)b ≤ 1147, which is 8. Next, we find 1147 − 128 × 8 = 1147 − 1024 = 123. Finally, with this we adjoin 21, i.e., 12321 and find largest possible integer c such that [2(10ab) + c]c ≤ 12321, i.e., (1360 + c)c ≤ 12321, which is 9, and the equality holds. Thus, √ 474721 = 689. Francois Viéte (1540-1603, France) noted that if one needs to calculate the square root of 2 to a high degree of accuracy, one should add as many zeros as necessary, and calculate the square root of, for example, 20, 000, 000, 000, 000, 000, 000, 000, 000, 000, 000, 000. That root he shows to be 141, 421, 356, 237, 309, 505, and thus the square root of 2 is approximately 1 41, 421, 356, 237, 309, 505 100, 000, 000, 000, 000, 000 .
We note that Aryabhata's Method explained above for 625 and 474721, combined with Viéte's observation easily computes the same approximation of √ 2, except instead of the last digit 5, we get 4; however, if we compute one more digit (which is 8) and then round it, then it is indeed 5.
To find the cube root of 1728, we group it in three's from right to left as 1 728. We search largest possible integer a such that a 3 ≤ 1, which is 1. This will be the first digit of the required cube root. Since 1 − 1 3 = 0, for the next digit we consider 728 and find largest possible integer b such that [3(10a) 2 + 3(10a)b + b 2 ]b ≤ 728, which is 2, and the equality holds. Thus, 3 √ 1728 = 12. To find the cube root of 12977875, we group it in three's from right to left as 12 977 875. We search largest possible integer a such that a 3 ≤ 12, which is 2. This will be the first digit of the required cube root. Now, we find 12 − a 3 = 12 − 2 3 = 4, and with this adjoin 977, i.e., 4977 and find largest possible integer b such that [3(10a) 2 + 3(10a)b + b 2 ]b ≤ 4977, which is 3. This will be the second digit of the required cube root. Next, we calculate 4977 − [3(10a) 2 + 3(10a)b + b 2 ]b = 4977 − 4167 = 810, and with this we adjoin 875, i.e., 810875. Finally, we find largest possible integer c such that [3(10ab) 2 + 3(10ab)c + c 2 ]c ≤ 810875, which is 5, and the equality holds. Thus, 3 √ 12977875 = 235. To find the cube root of 961504803, we group it in three's from right to left as 961 504 803. We search largest possible integer a such that a 3 ≤ 961, which is 9. This will be the first digit of the required cube root. Now, we find 961 − a 3 = 961 − 9 3 = 232, and with this adjoin 504, i.e., 232504 and find largest possible integer b such that [3(10a) 2 + 3(10a)b + b 2 ]b ≤ 232504, which is 8. This will be the second digit of the required cube root. Next, we calculate 232504 − [3(10a) 2 + 3(10a)b + b 2 ]b = 232504 − 212192 = 20312, and with this we adjoin 803, i.e., 20312803. Finally, we find largest possible integer c such that [3(10ab) 2 + 3(10ab)c + c 2 ]c ≤ 20312803, which is 7, and the equality holds. Thus, 3 √ 961504803 = 987. As for the square root, we can add as many zeros as necessary, and calculate the cube root with desired accuracy.
To find an approximate value of π, Aryabhatta gives the following prescription: Add 4 to 100, multiply by 8 and add to 62,000. This is "approximately" the circumference of a circle whose diameter is 20,000. This means π = 62,832/20,000 = 3.1416. It is important to

Great Pyramid at Gizeh and Rhind Mathematical Papyrus
From the dimensions of the Great Pyramid (erected around 2600 BC) it is possible to derive the two irrational numbers, namely, π, and the 'Golden Number' or 'divine proportion' Φ = (1 + √ 5)/2 1.61803398 · · · . A golden rectangle is a rectangle whose sides are in the ratio of Φ, it has dimensions pleasing to the eye and was used for the measurements of the facade of the Parthenon and other Greek temples, for details, see Sen and Agarwal [13]. Rhind Mathematical Papyrus (scribed 1650 BC) was scribed by Ahmes (1680-1620 BC, Egypt). It contains 87 problems, in particular, problem number 50 states that a circular field with a diameter of 9 units in area is the same as a square with sides of 8 units, i.e., π(9/2) 2 = 8 2 , and hence

Babylonians Tablet YBC 7289
There are numerous examples suggesting that Babylonians assembled large number of tables consisting of squares and square roots, and cubes and cubic roots. It has been suggested by several historians of mathematics, e.g., Victor Joseph Katz, (born 1942, USA) in his book [14] that "when square roots are needed in solving problems, the problems are arranged so that the square root is one that is listed in a table and is a rational number. However, where an irrational square root is needed, in particular, for √ 2, the result is generally written as 1; 25 = 1 5 12 ". On a fascinating tablet from Yale Babylonian Collection (YBC) number 7289 (around 1800-1600 BC), there is a scatter diagram of a square with side indicated as 30 and two numbers, see    The product of 30 by 1; 24, 51, 10 is exactly 42; 25, 35. Therefore, it is justifiable to presume that the number 42; 25, 35 represents the length of the diagonal and the number 1; 24, 51, 10 is √ 2. This confirms that Babylonians had enormous computational skills. The mathematical significance of this tablet was first recognized by the historians Otto Eduard Neugebauer (1899-1990, USA) and Abraham Sachs (1915( -1983. This tablet provides the correct value of √ 2 to six decimal digits. For further details, see Fowler and Robson [15]. The same Babylonian approximation of √ 2 was used later by Alexandrian Claudius Ptolemaeus, Ptolemy in English, (around 90-168 AD, Greek) in his Almagest, but he did not mention from where this approximation came, perhaps it was well known by his time. Carl Benjamin Boyer (1906Boyer ( -1976, USA) in his book [16] writes 1; 24, 51, 10 1.414222, which actually corresponds to 1; 24, 51, 12.
As in Sulbasutras there is no record how Babylonians obtained the approximations 1; 25 or 1; 24, 51, 10, of √ 2; however, definitely they must have realized that the exact value of √ 2 cannot be achieved. Thus, the methods which have been suggested by the historians are merely speculative. For example, Katz [14] believes that Babylonians used the algebraic identity (x + y) 2 = x 2 + 2xy + y 2 , which they might have perceived geometrically. Mathematically, the problem is for a given square of area N, we need to find its side √ N. For this, as a first step we select a regular number (evenly dividable of powers of 60) a close to, but less than, √ N (a good guess). Letting b = N − a 2 , the next step is to find c so that 2ac + c 2 is as close as possible to b, see Figure 5. If a 2 is "close enough" to N, then c 2 will be small in relation to 2ac, so c can be chosen equal to b/2a, that is, A similar argument shows that if a is greater than √ N, then √ N (a good guess). Letting b = N − a 2 , the next step is to find c so that 2ac + c 2 is as close as possible to b, see Figure 5. If a 2 is "close enough" to N, then c 2 will be small in relation to 2ac, so c can be chosen equal to b/2a, that is, A similar argument shows that if a is greater than √ N , then For √ 2, we begin with a = 1; 20 (= 4/3) < √ 2, (see (1)), to obtain a 2 = 1; 46, 40 (= 16/9), b = 0; 13, 20 (= 2/9) and b/2a = 0; 05 (= 1/12). Thus, from (8)  which is same as (1). Thus, we get all steps for √ 2 given in (1). Next, since 577/408 > √ 2, we again use (9), to obtain For √ 2, we begin with a = 1; 20 (= 4/3) < √ 2, (see (1)), to obtain a 2 = 1; 46, 40 (= 16/9), b = 0; 13, 20 (= 2/9) and b/2a = 0; 05 (= 1/12). Thus, from (8)  which is same as (1). Thus, we get all steps for √ 2 given in (1). Next, since 577/408 > √ 2, we again use (9), to obtain which is correct to 11 decimal places.
Thus, when we choose a < √ N, after applying (8), for further improvement we have to proceed to (9). Having this in mind, and looking (8) and (9), we can write the following algorithm (a process or set of rules to be followed in calculations) to compute √ N, also see (Boyer [16]), and (Ernst Sondheimer and Alan Rogerson [17]): a n+1 = 1 2 a n + N a n , n ≥ 0 where a 0 > 0 is any number (greater than or smaller than √ N), known as the initial approximation. Today algorithm (11) is derived by using Isaac Newton's (1642-1727, England) method: With appropriate x 0 the iterative scheme converges quadratically to a root of the general equation f (x) = 0. In our case the equation is f (x) = x 2 − N = 0. For N = 2 this is perhaps one of the oldest known algorithms. Historians Neugebauer and Sachs believed that the Babylonians obtained this algorithm for N = 2 based on the following principle: Suppose a is a guess which is too small (large), then 2/a will be a guess which is too large (small). Hence, their average (a + 2/a)/2 is a better approximation. This assumption that "divide and average" seems to be a general procedure of Babylonians for approximating square roots.
In the literature the algorithm (11) is also known as Heron's method after Heron of Alexandria (about 75 AD, perhaps an Egyptian) who gave the first explicit description of the method in his treatise Metrica which was discovered as recently as 1896 in Constantinople in a manuscript (the very word manuscript comes from the Latin words meaning 'written by hand') form dating from the eleventh or twelfth century. Heron used the formula (9), i.e., to calculate the square roots: "Since 720 has not a rational root, we shall make a close approximation to the root in this manner. Since the square nearest to 720 = N is 729 = a 2 , having a root a = 27, divide 27 into 720, i.e., N/a the result is N/a = 26 2 3 ; add a = 27, the result is N/a + a = 53 2 3 . Take half of this, i.e., 1 ; the result is 26 5 6 . Therefore the square root of 720 will be very nearly 26 5 6 . For 26 5 6 multiplied by itself gives 720 1 36 ; so that the difference is 1/36. if we wish to make the difference less than 1/36, instead of 729 we shall take the number now found 720 1 36 , and by the same method we shall find an approximation differing by much less than 1/36". Heron also found approximate square root of 63 also. The algorithm (11) generates a sequence {a n }, for which the concept of convergence was not existing even during the time of Heron. For the convergence of the Computation 2021, 9,29 12 of 49 sequence {a n } the following result is well-known, for example, see Agarwal et al. [18]: For the sequence {a n } the following hold 13) and the fact that |a 0 − √ N|/|a 0 + √ N| < 1. The convergence (quadratic) of this sequence to √ N immediately follows from (13). From (13), we also note that a n − √ N > 0 for all n ≥ 1. It also follows directly from the arithmetic-geometric mean inequality, in fact, for all n ≥ 1, we have a n = 1 2 a n−1 + N a n−1 ≥ a n−1 · N a n−1 with equality if and only if a 2 the sequence {a n } is decreasing. Thus, the sequence {a n }, in fact, converges monotonically. Jöran Friberg (born 1934, Sweden) in his book [19] mentions that Babylonian tablets (such as MS 3051) contain computations of areas of hexagons and heptagons, which involve the approximation of more complicated algebraic numbers (zeros of polynomials with integer coefficients) such as √ 3. The answer given there leads to the simple approximation √ 3 7/4. This does not mean they could not have calculated better approximations.
In Table 1, we use (11) to compute first three iterates for √ 2 and √ 3.  From Table 1 it is clear that the algorithm (11) gives both Sulbasutras approximations (1) and (5) of √ 2 and √ 3. It also gives Babylonian approximation √ 3 7/4. Unfortunately, from (11) we cannot get the Babylonians approximation (7) of √ 2. In fact, reversing a step in (11) leads to the equation which has only complex roots. Another simple explanation is 30547/21600 < √ 2, whereas a n > √ 2, n ≥ 1. We also note that Boyer in his book [16] has made a false assertion that a 2 with a 0 = 3/2 for √ 2 gives (7). In conclusion, Babylonians obtained (7) by some other unknown technique rather than (11), as has been claimed. A probable explanation for (7) is that Babylonians from their tables of n 2 and 2n 2 , n ≥ 1 noticed that 933119209 = (30547) 2 2(21600) 2 = 933120000. Algorithm (11) at n-th iteration requires division by a n , to avoid this we consider the equation f (x) = (1/x 2 ) − N = 0 and apply Newton's method (12), to get We multiply this by N and let a n = Nx n , to obtain a n+1 = a n 2N (3N − a 2 n ), n ≥ 0 which converges quadratically to √ N. For N = 2 with a 0 = 3/2 the above scheme gives a 1 = 45/32, a 2 = 185355/131072. These approximations of √ 2 are different from the corresponding entries in Table 1.
Problem xviii from the combined Babylonian tablet fragments BM 96957 and VAT 6598 gives two methods for calculating the diagonal d of a rectangle with sides of length a = 40 and b = 10 units. The first leads (in specific numbers) to the approximation d a + 2ab 2 3600 , (14) and the second method to the approximation is From Pythagorean Theorem d = √ 40 2 + 10 2 = 41.231056 · · · . Formulas (8), (9), (14) and (15), respectively, give the approximations The so-called Cairo Mathematical Papyrus, unearth in 1938 and first examined in 1962, dating from the early Ptolemeic dynasties (founded in 305 BC), contains 40 problems of mathematical nature. The problem particularly interesting in modern terms is to find the solution of the system of equations The scribe's method of solution amounts to adding and subtracting 2xy = 120 from the equation Furthermore, now employing (8), to obtain the approximations In an old Babylonian tablet (about 2000 BC) found in 1936 in Susa (Iraq), for the irrational number π the following expression appears 3 π = 57 60 + 36 (60) 2 , which yields π = 3 1/8 = 3.125. Babylonians were also satisfied with π = 3.

Shatapatha Brahmana
It means Priest manual of 100 paths (about 900 BC) is one of the prose texts describing the Vedic ritual. It survives in two recensions, Madhyandina and Kanva, with the former having the eponymous 100 brahmanas in 14 books, and the latter 104 brahmanas in 17 books. In these books π is approximated by 339/108 = 3.138888 · · · .

Pythagoreans (Followers of Pythagoras) Crisis of Incommensurability
Pythagoras (around 582-481 BC, Greece) is one of the most unexplained personalities in the history. He is among those individuals given the status of becoming a myth/omnipotent in his own lifetime. Since he followed typical oriental tradition (the knowledge was passed from one generation to the next mainly by word to mouth) whatever little we know about them is from the imaginations and a great many anecdotal fables thrown in by the historians who wrote (frequently contradict one another) and painted his picture hundreds of years after him, which continues even today. He has been called as mystic philosopher, master among masters, blend of genius and madness, mysterious, divinity, god-like figure, whereas some have shown doubt that such a person ever existed. In the book India in Greece, John J Griffin & Company, England, 1852, by the Greek historian Edward Pococke reports that Pythagoras, who taught Buddhist philosophy, was a great missionary. His name indicates his office and position; Pythagoras in English is equivalent to putha-gorus in Greek and Budha-guru in Sanskrit, which implies that he was a Buddhist spiritual leader. Note that Lord Gautama Buddha was during (1887-1807 BC), historians have misled the world by claiming that he flourished around (450 BC). Pythagoras is also considered to be a remarkably significant figure in the advancement of mathematics, science, and pre-Socratic philosophy (the study of the fundamental nature of knowledge, reality, and existence, the word philosophy is due to Socrates, around 469-399 BC), even though we know comparatively little about his mathematical achievements. In any case, for his many accomplishments in mathematics for which he is being credited, in recent years it has been shown that these were already known several centuries before him. For example, for the origin of Pythagorean Theorem (which made Pythagoras immortal) and Pythagorean Triples see Agarwal [20,21]. Still, the Pythagorean legacy lasted well over more than two and a half millenniums, and continues to be present in the modern day, starting from high school students. His philosophy appeared suddenly and unexpectedly in Albert Einstein's (1879-1955) formulation of the general theory of relativity. Today, Pythagoras is revered as a prophet by the Ahl-al-Tawhid or Druze (a concept, upon which a Muslim's entire faith rests) along with Greek philosopher Plato of Athens (around 427-347 BC). Plato (meaning broad) is a nickname, his real name was Aristocles, he died at a wedding feast.
Pythagoras gave 'divine significance' to most natural numbers, and attempted to find mathematical explanations for everything in the universe in terms of rational numbers "possibly the most mischievous misreading of nature in the history of human error" (Eric Temple Bell, 1883-1960, USA, Britain). He paid homage to every numerical relationship such as equation and inequality (arithmetic then). His motto was All is Number, "numbers Rules the Universe", "number is the ruler of forms and ideas and the cause of gods and demons". He identified some human attribute to most numbers, such as even numbers he regarded as feminine, pertaining to the earthly; odd numbers as masculine, partaking of celestial nature. However, the hypotenuse of a most obvious right-angle triangle with the same legs lead to the the number √ 2, which Pythagoreans could not write as a rational number. The discovery of incommensurability of √ 2 caused tremendous crisis/confusion/devastation/surprise/shattering effect among the Pythagoreans, for it challenged the adequacy of their basic philosophy that number was the essence of everything. In fact, in the numerical sense, the universe was seen to be irrational. This logical calamity enforced them to maintain the pledge of strict secrecy. To incommensurable numbers they named as "the unutterable", (Greeks used the term logos, meaning word or speech, for the ratio of two integers, when incommensurable lengths were described as alogos, the term carried a double meaning: not a ratio and not to be spoken) as it was a dangerous secret to possess. According to a legend, Hippasus of Metapontum (about 500 BC, Greece), a Pythagorean was murdered-thrown off a ship to drown at sea by fanatic Pythagoreans, because he uttered the unutterable to an outsider (some historians have speculated that Hippasus had first proof of the existence of irrational numbers), whereas others say he lost his fortune and tried to recoup his losses by teaching the doctrine of irrational numbers. Anyway it is hard to keep a secret in science. This revelation/achievement of Pythagoreans, that not all numbers are rational marked, is considered one of the most fundamental discoveries in the entire history of science (it evolved the number concept by filling the gaps which were there between rationals). Historians have also argued that this major discovery also helped in the development of deductive reasoning. However, it seems to be inexplicable as we have noted in Section 2 that irrationality of √ 2 was already conjectured in Sulbasutras, and several later Sanskrit scholars decisively claimed that irrationality of √ 2 was first discovered by ancient Hindus. In fact, from the following quotes it confirms that Pythagoras learnt about the irrationality of √ 2 in India. Francois Marie Arouet Voltaire (1694-1778), one of the greatest French writers and philosophers: "I am convinced that everything has come down to us from the banks of Ganga-Astronomy, Astrology, and Spiritualism. Pythagoras went from Samos to Ganga 2600 years ago to learn Geometry. He would not have undertaken this journey had the reputation of the Indian science had not been established before." Thomas Stearns Eliot (1888-1965), American-British poet, Nobel Laureate (1948): "I am convinced that everything has come down to us from the banks of the Ganga-Astronomy, Astrology, Spiritualism, etc. It is very important to note that some 2500 years ago at the least Pythagoras went from Samos to the Ganga to learn Geometry but he would certainly not have undertaken such a strange journey had the reputation of the Brahmins' science not been long established in Europe".
In 2007, Borzacchini [22] has asserted that Pythagorean music theory is the origin of incommensurability.

Theodorus of Cyrene (about 431 BC, Libya, Greece)
He is said to have been Plato's teacher. From the dialogues of Plato, we know that Theodorus demonstrated geometrically that the sides of squares represented by , and √ 17, are incommensurable with a unit length. That is, he showed the irrationality of the square roots of nonsquare integers from 3 to 17, 'at which point', says Plato, "for some reason he stopped", see Figure 6.
Theodorus of Cyrene (about 431 BC, Libya, Greece). He is said to have been Plato's teacher. m the dialogues of Plato, we know that Theodorus demonstrated geometrically that the sides of squares resented by , and √ 17, are incommensurable with nit length. That is, he showed the irrationality of the square roots of nonsquare integers from 3 to 17, 'at ch point,' says Plato, "for some reason he stopped," see Figure 6. It has been speculated that Theodorus constructed his spiral based on right triangles with a common tex, where in each triangle the side opposite the common vertex has length 1. The hypotenuse of the triangle then has length √ n + 1, follows immediately by Pythagorean Theorem. His spiral also suggest sible reason Theodorus stopped at √ 17: On summing of the vertex angles for the first n triangles, we e n = 16 (which gives √ 17) this sum is 351.15 0 , while for n = 17 the sum is 364.78 0 . Thus for n > 16 his al started to overlap itself (i.e., cuts the initial axis for the first time) and the drawing became "messy." aetetus (around 417-369 BC, Greece), who was a pupil of Theodorus and a member of Plato's school in ens, extended the result, demonstrating that the square root of any nonsquare integer is irrational, and cube root of any number that is not a perfect cube is irrational. Of course, today, by induction one can w √ n for any n. Also, if n is an odd integer, then √ n can be represented by the leg of a right triangle se hypotenuse is (n + 1)/2 and whose leg is (n − 1)/2, i.e., ( Further, is an even integer, then √ n can be represented by half of the leg of a right triangle whose hypotenuse + 1 and whose other leg is n − 1, i.e., (2 √ n) 2 = (n + 1) 2 − (n − 1) 2 . Plato himself also showed that a ional number could be the sum of two irrationals. In Figure 7, we provide the construction of √ 5 and √ 6 metrically. It has been speculated that Theodorus constructed his spiral based on right triangles with a common vertex, where in each triangle the side opposite the common vertex has length 1. The hypotenuse of the nth triangle then has length √ n + 1, follows immediately by Pythagorean Theorem. His spiral also suggest possible reason Theodorus stopped at √ 17: On summing of the vertex angles for the first n triangles, we have For n = 16 (which gives √ 17) this sum is 351.15 0 , while for n = 17 the sum is 364.78 0 . Thus, for n > 16 his spiral started to overlap itself (i.e., cuts the initial axis for the first time) and the drawing became "messy". Theaetetus (around 417-369 BC, Greece), who was a pupil of Theodorus and a member of Plato's school in Athens, extended the result, demonstrating that the square root of any nonsquare integer is irrational, and the cube root of any number that is not a perfect cube is irrational. Of course, today, by induction one can draw √ n for any n. Furthermore, if n is an odd integer, then √ n can be represented by the leg of a right triangle whose hypotenuse is (n + 1)/2 and whose leg is (n − 1)/2, i.e., ( Further, if n is an even integer, then √ n can be represented by half of the leg of a right triangle whose hypotenuse is n + 1 and whose other leg is n − 1, i.e., (2 √ n) 2 = (n + 1) 2 − (n − 1) 2 . Plato himself also showed that a rational number could be the sum of two irrationals. In Figure 7, we provide the construction of √ 5 and √ 6 geometrically.  Although there are speculations that incommensurability of √ 2 was proved by Pythagoras himself and for all nonsquare integers by Theodorus, the first fully geometric proof appeared in the Meno (Socratic dialogue by Plato).
Following the Website http://mitp-content-server.mit.edu:18180/books/content/ sectbyfn?collid=books_press_0&id=1043&fn=9780262661829_schh_0001.pdf (accessed on 3 March 2021), in the square ABCD we use a compass to cut off AF = AD along the diagonal CA. At F draw the perpendicular EF (see Figure 8). Then the ratio of CE to CF (hypotenuse to side) will be the same as the ratio of AC to AD, since the triangles CDA and EFC are similar. Suppose that DC and CA were commensurable. Then there would be a segment δ such that both DC and CA were integral multiples of δ. Since AF = AD, then CF = CA − AF is also a multiple of δ. Note also that CF = EF, because the sides of triangle EFC correspond to the equal sides of triangle CDA. Further, EF = DE because (connecting A and E) triangles EDA and EFA are congruent. Thus, DE = CF is a multiple of δ. Then CE = CD − DE is also a multiple of δ. Therefore, both the side CF and hypotenuse CE are multiples of δ, which therefore is a common measure for the diagonal and side of the square of side CF. The process can now be repeated as follows: on EC cut off EG = EF and construct GH perpendicular to CG. The ratio of hypotenuse to side will still be the same as it was before and hence the side of the square on CG and its diagonal also share δ as a common measure. Because we can keep repeating this process, we will eventually reach a square whose side is less than δ, contradicting our initial assumption. Therefore, there is no such common measure δ. The demonstration given here has been named as The Method of Infinite Descent, and it has been credited to Pierre de Fermat (1601-1665, France). In fact, in 1879, a paper was found in the library of Leyden, among the manuscript of Christiaan Huygens (1629-1695, Netherlands), in which Fermat describes this method by which he may have made many of his discoveries. The method is particularly useful in establishing negative results, but often difficult to apply.  wing the Website http://mitp-content-server.mit.edu:18180/books/content/sectbyfn?collid=books p 0&id=1043&fn=9780262661829 schh 0001.pdf, in the square ABCD we use a compass to cut off AF = long the diagonal CA. At F draw the perpendicular EF (see Figure 8). Then the ratio of CE to CF otenuse to side) will be the same as the ratio of AC to AD, since the triangles CDA and EF C are ar. Suppose that DC and CA were commensurable. Then there would be a segment δ such that both and CA were integral multiples of δ. Since AF = AD, then CF = CAAF is also a multiple of δ. also that CF = EF, because the sides of triangle EF C correspond to the equal sides of triangle CDA. er, EF = DE because (connecting A and E) triangles EDA and EF A are congruent. Thus, DE = CF ultiple of δ. Then CE = CDDE is also a multiple of δ. Therefore, both the side CF and hypotenuse re multiples of δ, which therefore is a common measure for the diagonal and side of the square of side The process can now be repeated as follows: on EC cut off EG = EF and construct GH perpendicular . The ratio of hypotenuse to side will still be the same as it was before and hence the side of the square • The following inquisitive geometric proof of Apostol [23] (also for similar proofs see earlier books by Kiselev [24], and Conway and Guy [25]) is in line with the above proof. A circular arc with center at the uppermost vertex and radius equal to the vertical leg of the triangle intersects the hypotenuse at a point, from which a perpendicular to the hypotenuse is drawn to the horizontal leg (see Figure 9). Each line segment in the diagram has integer length, and the three segments with double tick marks have equal lengths. (Two of them are tangents to the circle from the same point). Therefore the smaller isosceles right triangle with hypotenuse on the horizontal base also has integer sides.
The following inquisitive geometric proof of Apostol [7] (also for similar proofs see earlier books by Kiselev 8], and Conway and Guy [13]) is in line with the above proof. A circular arc with center at the uppermost ertex and radius equal to the vertical leg of the triangle intersects the hypotenuse at a point, from which perpendicular to the hypotenuse is drawn to the horizontal leg (see Figure 9). Each line segment in the iagram has integer length, and the three segments with double tick marks have equal lengths. (Two of em are tangents to the circle from the same point.) Therefore the smaller isosceles right triangle with ypotenuse on the horizontal base also has integer sides. is contributions include: a mathematical theory of "magnitudes"-such as lengths, areas, volumes; addition f numerous results on the study of golden section; invention of a process known as the method of exhaustion; nd the theory of proportion, partly to place the doctrine of incommensurables upon a thoroughly sound asis. The irrationality of the square root of two Eudoxus phrased as "a diagonal and a side of a square have o common measure". He realized that an irrational is known by the rational numbers less than it, and the tional numbers greater than it. This task was done so well that Greek mathematicians made tremendous rogress in geometry and it survived as Book V of Euclid's Elements. It still continues, fresh as ever, after e great arithmetical reconstructions of Julius Wilhelm Richard Dedekind (1831-1916, Germany) and Karl heodor Wilhelm Weierstrass (1815-1897, Germany) during the nineteenth century.
2. Aristotle (around 384-322 BC, Greece). The first semi-geometric proof of the irrationality of √ 2 is ue to Aristotle which appeared in his Analytica Priora. He concludes that if the side and the diagonal are ssumed commensurable, then odd numbers are equal to even numbers. For this, he used the method of ontradiction: Suppose that the side EH and the diagonal HF, see Figure 10, are commensurable, i.e., each an be expressed by the number of times it is measured by their common measure. Now it can be assumed at at least one of these numbers is odd, if not there would be a longer common measure. Then the squares EF G and ABCD on the side and diagonal, respectively, represent square numbers. From the Figure 10, is clear that the area of the latter square is clearly double the former, thus it represents an even square umber. Consequently, its side AB = HF is also an even number, and thus the square ABCD is a multiple f four. Finally, since HEF G is half of ABCD, it must be a multiple of two, i.e, it is also an even square. herefore its side EH must also be even. However, this contradicts the original assumption that one of 17 Figure 9. Incommensurability of √ 2 by Apostol.

Eudoxus of Cnidus (around 400-347 BC, Greece)
He was the most celebrated mathematician. His contributions include: a mathematical theory of "magnitudes"-such as lengths, areas, volumes; addition of numerous results on the study of golden section; invention of a process known as the method of exhaustion; and the theory of proportion, partly to place the doctrine of incommensurables upon a thoroughly sound basis. The irrationality of the square root of two Eudoxus phrased as "a diagonal and a side of a square have no common measure". He realized that an irrational is known by the rational numbers less than it, and the rational numbers greater than it. This task was done so well that Greek mathematicians made tremendous progress in geometry and it survived as Book V of Euclid's Elements. It still continues, fresh as ever, after the great arithmetical reconstructions of Julius Wilhelm Richard Dedekind (1831-1916, Germany) and Karl Theodor Wilhelm Weierstrass (1815-1897, Germany) during the nineteenth century.

Aristotle (around 384-322 BC, Greece)
The first semi-geometric proof of the irrationality of √ 2 is due to Aristotle which appeared in his Analytica Priora. He concludes that if the side and the diagonal are assumed commensurable, then odd numbers are equal to even numbers. For this, he used the method of contradiction: Suppose that the side EH and the diagonal HF, see Figure 10, are commensurable, i.e., each can be expressed by the number of times it is measured by their common measure. Now it can be assumed that at least one of these numbers is odd, if not there would be a longer common measure. Then the squares HEFG and ABCD on the side and diagonal, respectively, represent square numbers. From Figure 10, it is clear that the area of the latter square is clearly double the former, thus it represents an even square number. Consequently, its side AB = HF is also an even number, and thus the square ABCD is a multiple of four. Finally, since HEFG is half of ABCD, it must be a multiple of two, i.e, it is also an even square. Therefore its side EH must also be even. However, this contradicts the original assumption that one of HF, EH is odd. In conclusion, the two lines EH and HF are incommensurable. Thus, Aristotle in number theory succeeded in proving the existence of irrationals.
HF, EH is odd. In conclusion, the two lines EH and HF are incommensurable. Thus, Aristotle in number theory succeeded in proving the existence of irrationals.  Figure 10, it is clear that the area of ABCD is the same as two times the area of HEF G. This construction is due to Socrates (around 469-399, Greece) in the Meno. Socrates is considered as one of the founders of Western philosophy, he was sentenced to death by the drinking of a mixture containing poison hemlock, because he was found guilty of corrupting the minds of the youth of Athens and of impiety "not believing in the gods of the state." 13. Euclid of Alexandria (around 325-265 BC, Greece, Egypt). His masterpiece work The Elements is divided into 13 books (each about the length of a modern chapter) and contains 465 propositions on plane and solid geometry, and number theory. In compiling the Elements, Euclid organized deductively on the basis of explicit axioms the experience and achievements of his predecessors of three centuries just past. Euclid's semi-geometrical demonstration by the method of contradiction of the irrationality of √ 2 is given in Book 10, Proposition 27. Though it is less perspicuous than the strictly arithmetical proof current today, it is more suggestive historically, and more precise than Aristotle's proof, see Section 12. The argument goes as follows: If the diagonal AC and side AB of the square ABCD (see Figure 10) have a common measure, say δ, then there exist p, q ∈ IN satisfying AC = pδ, AB = qδ. The ratio of these segments is In what follows, we can assume that common factors of p and q have been cancelled, i.e., gcd(p, q) = 1. Thus, at least one them is odd. Squaring the identity (16), we have Now in view of Pythagorean theorem in the triangle ABC, we find (AC) 2 = 2(AB) 2 , so that (17) is the same as Now since 2q 2 is an even integer, p 2 must also be even. But, then p is also even, i.e., p = 2k. Substituting this in the equation p 2 = 2q 2 gives q 2 = 2k 2 . But then q 2 and hence q is also an even number. In conclusion, both p and q are even, which contradicts our initial assumption that they have no common factor, or one of them is odd.
From Figure 10, it is clear that the area of ABCD is the same as two times the area of HEFG. This construction is due to Socrates (around 469-399, Greece) in the Meno. Socrates is considered as one of the founders of Western philosophy, he was sentenced to death by the drinking of a mixture containing poison hemlock, because he was found guilty of corrupting the minds of the youth of Athens and of impiety "not believing in the gods of the state".

Euclid of Alexandria (around 325-265 BC, Greece, Egypt)
His masterpiece work The Elements is divided into 13 books (each about the length of a modern chapter) and contains 465 propositions on plane and solid geometry, and number theory. In compiling the Elements, Euclid organized deductively on the basis of explicit axioms the experience and achievements of his predecessors of three centuries just past. Euclid's semi-geometrical demonstration by the method of contradiction of the irrationality of √ 2 is given in Book 10, Proposition 27. Though it is less perspicuous than the strictly arithmetical proof current today, it is more suggestive historically, and more precise than Aristotle's proof, see Section 12. The argument goes as follows: If the diagonal AC and side AB of the square ABCD (see Figure 10) have a common measure, say δ, then there exist p, q ∈ IN satisfying AC = pδ, AB = qδ. The ratio of these segments is In what follows, we can assume that common factors of p and q have been cancelled, i.e., gcd(p, q) = 1. Thus, at least one them is odd. Squaring the identity (16), we have Now in view of Pythagorean theorem in the triangle ABC, we find (AC) 2 = 2(AB) 2 , so that (17) is the same as Now since 2q 2 is an even integer, p 2 must also be even. However, then p is also even, i.e., p = 2k. Substituting this in the equation p 2 = 2q 2 gives q 2 = 2k 2 . However, then q 2 and hence q is also an even number. In conclusion, both p and q are even, which contradicts our initial assumption that they have no common factor, or one of them is odd.

•
In the above proof we can ignore all geometric arguments, and directly proceed to algebraic Equation (18) [26] of 1920, and later Gardner [27], and Laczkovich [28] in their books assume that in √ 2 = p/q the integer q is the smallest possible such number. Their main argument is essentially to use the equality (2q − p) 2 = 2(p − q) 2 which is true if and only if (18) holds. Thus, it follows that This contradicts the minimality of q.
It is interesting to note that • Rademacher and Toeplitz in their book of 1957 ( [29], Chapter 4) assert that (18) implies p is even, so q must be odd. However, the square of an even number is divisible by 4, which leads to conclude that q must be even. Thus, we have Aristotle type contradiction.
Now to prove (18) we shall apply the Fundamental Theorem of Arithmetic (FTOA). Euclid's Elements Book VII, propositions 30, 31 and 32, and Book IX, proposition 14 substantiate the statement and proof of the FTOA. Although Euclid felt that irrational numbers simply did not belong in a work based on arithmetic, some authors claim that Euclid in Book X, Proposition 117 uses FTOA to almost show the impossibility of (18), but most of the English transactions of Elements have only 115 propositions. Fritz [30] indicates that the early Greek mathematicians did not explicitly use the FTOA to prove the irrationality of √ 2. In fact, on the Website http://people.math.harvard.edu/mazur/preprints/Eva.Nov.20.pdf, accessed on 4 March2021, posted in 2005, Mazur claims that the explicit use of the FTOA is post Karl Friedrich Gauss (1777-1855, Germany). We state the modern version of this theorem in the following: Fundamental Theorem of Arithmetic. Every integer n ≥ 2 is either prime or can be expressed as a product of primes: that is, n = p 1 p 2 · · · p m , where p 1 , p 2 , · · · , p m are primes. Furthermore, this factorization is unique except possibly for the order in which the factors occur.
Now among the primes p i and q i , the prime 2 may occur (it will occur if either p or q is even). If it does occur, it must appear an even number of times on the left side of Equation (19) (since each prime there appears twice), and an odd number of times on the right side (because 2 already appears there once). However, then we have a contradiction: since the factorization into primes is unique, the prime 2 cannot appear an even number of times on one side of the equation and an odd number on the other. Thus, Equation (18) is impossible.
• From the uniqueness of the factorization, one can argue directly that p 2 has even number of prime factors, whereas 2q 2 has odd number of prime factors, which is absurd. • Some of the above illustrations can be extended to prove the result: If N ∈ IN, then √ N is a rational number if and only if √ N is an integer. First, we model its proof due to Gardner [27]. Clearly, if √ N is an integer, then √ N is rational. Conversely, we assume that √ N is rational, i.e., it can be written as √ N = p/q, where p, q ∈ IN and q is the smallest possible such integer. Let k = [ √ N], where [·] is the usual greatest integer function. Then, it follows that k < p/q < k + 1, and therefore 0 < p − kq < q. Now note that the equality (Nq − kp) 2 = N(p − kq) 2 is true if and only if p 2 = Nq 2 holds. Thus, but this contradicts the fact that q is the smallest. Now we will apply FTOA. Again if √ N is an integer, then √ N is rational. Conversely, we assume that √ N is rational, i.e., it can be written as √ N = p/q, where p, q ∈ IN and gcd(p, q) = 1. Since p/q is not an integer, q ≥ 2. Again, we have p 2 = Nq 2 . By FTOA, q has a prime factor m. Thus, m|Nq 2 and so m|p 2 , but then m|p. Hence, m|p and m|q, which contradicts our assumption that gcd(p, q) = 1.

•
Dedekind in his proof assumed that if N is not a square of an integer, then there exists a positive integer λ such that λ 2 < N < (λ + 1) 2 . Again, if N is rational, then there exist p, q ∈ IN such that p 2 − Nq 2 = 0, where q is the least possible integer possessing the property that its square multiplied by N is the square of p. Since λq < p < (λ + 1)q, it follows that the integers s = p − λq and t = Nq − λp are positive, and we have t 2 − Ns 2 = (λ 2 − N)(p 2 − Nq 2 ) = 0, which contradicts the assumption on q. • On the Website https://www.quora.com/If-p-is-a-natural-number-but-not-a-perfectnth-power-how-does-one-prove-that-the-nth-root-of-p-is-not-rational (accessed on 3 March 2021), Thomas Schürger (2019) has provided a very simple proof of the following general result: The kth, k ∈ IN, k ≥ 2 root of a nonnegative integer N ≥ 2 is rational if and only if N is a perfect kth power. One direction of this statement is clearly true: the kth root of a kth power is rational. Let us prove the other direction via proof by contradiction. Let us assume that N is not a perfect kth power, and k √ N is rational, i.e., k √ N = p q for some p, q in IN such that p/q is in lowest terms. Since and p/q is in lowest terms p k /q k is also in lowest terms, and N/1 is clearly in lowest terms. It follows that p k = N and q k = 1, which is a contradiction since we assumed that N is not a perfect kth power. Hence, k √ N must be an irrational number. • Some of the above arguments need slight modification to prove: If r and s are distinct primes, then √ rs and log r s are irrational. For example, to show log r s is irrational, we assume contrary, i.e., log r s = p/q, where p, q ∈ IN. We can assume that gcd(p, q) = 1. Then r p/q = s and so (r p/q ) q = s q . Therefore, r p = s q . Since r|r p , it follows that r|s q and so r|s, which is a contradiction. • We shall follow Dov Jarden (1911-1986, Israel) work of 1953 to show that there exist irrational numbers a and b such that a b is rational. Consider the irrational numbers a = b = √ 2. If the number a b = √ 2 √ 2 is rational, we are done. If √ 2 √ 2 is irrational, we consider the numbers a = √ 2 Note that in this proof we could not find irrational numbers a and b such that a b is rational.

Archimedes of Syracuse (287-212 BC, Greece)
He is considered as one of three complete mathematicians world has so far produced (the other two are Newton and Gauss). Archimedes developed a general method of exhaustion, specially to approximate the value of π. His method is based on the following arguments: the circumference of a circle lies between the perimeters of the inscribed and circumscribed regular polygons (equilateral and equiangular) of n sides, and as n increases, the deviation of the circumference from the two perimeters becomes smaller. If a n and b n denote the perimeters of the inscribed and circumscribed regular polygons of n sides, and C the circumference of the circle, then it is clear that {a n } is an increasing sequence bounded above by C, and {b n } is a decreasing sequence bounded below by C. Both of these sequences converge to the same limit C. For simplicity, we choose a circle with the radius 1, then from Figure 11 it immediately follows that Note that in this proof we could not find irrational numbers a and b such that a b is rational.
14. Archimedes of Syracuse (287-212 BC, Greece). He is considered as one of three complete mathematicians world has so far produced (the other two are Newton and Gauss). Archimedes developed a general method of exhaustion, specially to approximate the value of π. His method is based on the following arguments: the circumference of a circle lies between the perimeters of the inscribed and circumscribed regular polygons (equilateral and equiangular) of n sides, and as n increases, the deviation of the circumference from the two perimeters becomes smaller. If a n and b n denote the perimeters of the inscribed and circumscribed regular polygons of n sides, and C the circumference of the circle, then it is clear that {a n } is an increasing sequence bounded above by C, and {b n } is a decreasing sequence bounded below by C. Both of these sequences converge to the same limit C. For simplicity, we choose a circle with the radius 1, then from Figure  11 it immediately follows that Figure 11. Archimedes Approximation of π 1 α/2 sin α/2 α/2 1 tan α/2 a n = n sin π n and b n = n tan π n .
It is clear that lim n→∞ a n = π = lim n→∞ b n . Further, b 2n is the harmonic mean of a n and b n , and a 2n is the geometric mean of a n and b 2n , i.e., b 2n = 2a n b n a n + b n and a 2n = a n b 2n .
From (20) for the hexagon, i.e., n = 6 it follows that a 6 = 3, b 6 = 2 √ 3. Then, Archimedes successively took polygons of sides 12, 24, 48, and 96, used the recursive relations (21), and the inequality to obtain the bounds 3.140845 · · · = 3 10 71 The approximation 22/7 is often called the Archimedean value of π, and it is good for most purposes. Archimedes' polygonal method remained unsurpassed until 18 centuries, see Agarwal et. al. [2]. The inequality (22) is of paramount interest because the bounds 265/153 and 1351/780 are best rational approximations up to the respective denominators. The following rational bounds for π where either the lower 21 Figure 11. Archimedes Approximation of π. a n = n sin π n and b n = n tan π n .
It is clear that lim n→∞ a n = π = lim n→∞ b n . Further, b 2n is the harmonic mean of a n and b n , and a 2n is the geometric mean of a n and b 2n , i.e., b 2n = 2a n b n a n + b n and a 2n = a n b 2n .
From (20) for the hexagon, i.e., n = 6 it follows that a 6 = 3, b 6 = 2 √ 3. Then, Archimedes successively took polygons of sides 12, 24, 48, and 96, used the recursive relations (21), and the inequality to obtain the bounds 3.140845 · · · = 3 10 71 The approximation 22/7 is often called the Archimedean value of π, and it is good for most purposes. Archimedes' polygonal method remained unsurpassed until 18 centuries, see Agarwal et al. [10]. The inequality (22) is of paramount interest because the bounds 265/153 and 1351/780 are best rational approximations up to the respective denominators. The following rational bounds for π where either the lower bound or the upper bound is the best k-digit rational approximation are obtained in Sen et al. [31]  One of the most frequently debated questions in the history of mathematics is the "puzzling" approximation of √ 3, appeared in his book Measurement of a Circle, namely, the inequality (22) which Archimedes presented without a justification. On the Website https://mathpages.com/home/kmath038/kmath038.htm (accessed on 3 March 2021), for the inequality (22) several reviews which appeared in the popular history of mathematics books have been summarized, for example: Walter William Rouse Ball (1850-1925, Britain) in 1908 "it would seem...that [Archimedes] had some (at present unknown) method of extracting the square root of numbers approximately", Thomas Little Heath (1861-1940, Britain) in 1921 "the successive solutions in integers of the equations x 2 − 3y 2 = 1 and x 2 − 3y 2 = −2 may have been found...in a similar way to...the Pythagoreans", Bell in 1937, "...he also gave methods for approximating to square roots which show that he anticipated the invention by the Hindus of what amount to periodic continued fractions", Boyer in 1968, "his method for computing square roots was similar to that used by the Babylonians", Morris Kline (1908-1992, USA) in 1972, without any explanation claimed that if N = a 2 ± b where a 2 is the rational square nearest to N, larger or smaller, and b is the remainder, then the following inequalities can be used to obtain (22) As we have seen the right side bounds of the inequality (24) lead to the algorithm (11) which indeed gives the upper bound of (22) (see Table 1, N = 3, a 0 = 5/3), the left side bounds of (24) give us two new iterative schemes and a n+1 = a n − a 2 For (25), by induction, we shall show that a n ≤ √ N < a n + 1 implies that a n+1 ≤ √ N < a n+1 + 1, n ≥ 0. For this, it suffices to show that a 2 n + a n + N 2a n + 1 ≤ √ N < a 2 n + a n + N 2a n + 1 + 1 or (a n − √ N)(a n + 1 − √ N) ≤ 0 < (a n − √ N) 2 + (a n + 1 − √ N) + 2a n , which in view of a n ≤ √ N < a n + 1 is obvious. From (25), we also have a n ≤ a n+1 . Thus, the sequence {x n } generated by (25) is monotonically increasing, and bounded above, and hence converges to √ N. For the sequence {x n } generated by the iterative scheme (26) numerical evidence suggests that the convergence is oscillatory. Further, from (25) as well as (26) we could not get the lower bound of (22), see Table 2. Since the error is negated and reduced by a factor of nearly 52 in each iteration. Iterative scheme (27)  Thus, a 2 and a 3 , respectively, give the lower and upper Archimedes bounds of √ 3.
• On the same Website and on the Website https://www.mathpages.com/home/ kmath190/kmath190.htm (accessed on 3 March 2021), following Babylonians' the basic ladder rule for generating a sequence of integers to yield the square root of a number N the following recurrence relation has been discussed where a is the largest integer such that a 2 is less than N. Letting q = √ N + a, or (q − a) 2 = N, it follows that and hence s 0 = 1, s 1 = q, s n = q n , n ≥ 2 satisfies (31). Now since q = s n−1 /s n−2 and q 2 = s n /s n−2 from (31) it immediately follows that q = s n+1 /s n , n ≥ 0. However, since exactly q is unknown, we can begin with arbitrary (initial) integer values of s 0 , s 1 and generate the sequence of the ratios {s n+1 /s n }, which must converge to the solutions of (31), namely, q = √ N + a. Thus, {(s n+1 /s n ) − a} converges to √ N. We also note that {(N − a 2 )(s n /s n+1 )} converges to (N − a 2 )/q = (N − a 2 )/( √ N + a) = √ N − a, and hence {a + (N − a 2 )(s n /s n+1 )} converges to √ N. Now we shall show that for both the sequences {(s n+1 /s n ) − a} and {a + (N − a 2 )(s n /s n+1 )} convergence is oscillatory. For the first sequence it suffices to show that if (s n+1 /s n ) − a > √ N, which is the same as (s n /s n+1 ) < 1/( √ N + a), then (s n+2 /s n+1 ) − a < √ N. For this, from (31) we have Similarly, for the second sequence it suffices to show that if a + (N − a 2 )(s n /s n+1 ) > √ N, which is the same as (s n /s n+1 ) > 1/( However, this is the same as proving (s n+2 /s n+1 ) > √ N + a. Now from (31) it follows that For N = 2 and N = 3, we need to take a = 1, so that the recurrence relation (31), respectively, reduces to s n = 2s n−1 + s n−2 and s n = 2s n−1 + 2s n−2 , n ≥ 2. We shall consider these recurrence relations with s 0 = 0 and s 1 = 1, i.e., s n = 2s n−1 + s n−2 , n ≥ 2, s 0 = 0, s 1 = 1 (32) and s n = 2s n−1 + 2s n−2 , n ≥ 2, s 0 = 0, s 1 = 1.
Although solutions of (32) and (33) can be written explicitly as for the computation they are of little help. In Table 3, we directly use (32) and (33) to list successive approximations obtained for √ 2 and √ 3. Table 3 contains most of the data of Table 1, also it includes Archimedes' lower and upper bounds for √ 3, in fact, it is probable that Archimedes used iterative scheme (31) to establish the inequality (22).
and an argument similar to that of bisection method to compute Archimedes lower and upper bounds in (22). For this, he assumed a pair of two approximations α = v/u and β = y/x of √ 3 such that α < √ 3 < β. Now calculate γ = (v + y)/(u + x) and replace α by γ if γ < √ 3, i.e., (v + y) 2 < 3(u + x) 2 , and replace β by γ if γ > √ 3, i.e., (v + y) 2 > 3(u + x) 2 . This gives an improved pair of approximations. The procedure continues until the desired accuracy is achieved. With α = 1 and β = 2 his first sixteen pairs of approximations are 1  a n+1 = a n−1 f (a n ) − a n f (a n−1 ) f (a n ) − f (a n−1 ) , n ≥ 1 (34) where a 0 , a 1 are two initial approximations, one is less than a * and the other is greater than a * . For the root a * the secant method is superlinear, i.e., the rate of convergence is the Golden Number Φ. We note that for the equation f (x) = x 2 − N = 0 the secant method (34) simply reduces to a n+1 = a n−1 a n + N a n−1 + a n , n ≥ 1.
It is interesting to note that if in (35), we take a n−1 = a n , then it is the same as (11). Applying (35) with N = 3, a 0 = 5/3 (which is less than √ 3), and a 1 = 26/15 (which is greater than √ 3), see Table 1, we immediately get a 2 = 265/153, which is the lower bound in (22). From (35), we also compute a 3 = 13775/7953 1.73205079844, which is a better lower bound than in (22). Based on the above inequalities, we have the following three algorithms HMA, GMA, and AMA c n+1 = 2a n b n a n + b n , b n+1 = a n b n , a n+1 = a n + b n 2 , n ≥ 0 where a 0 , b 0 are positive (initial approximation) numbers. The GMA and AMA first appeared in the works of Lagrange, and their properties were further analyzed by Gauss, for their applications to approximate π see the recent monograph of Chan [33]. It is clear that c n+1 ≤ b n+1 ≤ a n+1 , n ≥ 0. From this, it immediately follows that a n+1 − a n = a n + b n 2 − a n = b n − a n 2 ≤ 0, b n+1 = a n b n ≥ b n b n = b n , and c n+1 − c n = 2a n b n a n + b n − c n ≥ 2a n b n 2a n − c n = b n − c n ≥ 0, thus the sequence {a n } is decreasing, the sequence {b n } is increasing, the sequence {c n } is also increasing and b n ≤ c n+1 ≤ b n+1 . Thus, min{a 0 , b 0 } ≤ c 1 ≤ c n ≤ b n ≤ a n ≤ a 1 ≤ max{a 0 , b 0 }. In conclusion all the three sequences {c n }, {b n }, {a n } converge to the same limit. The convergence of {c n } also follows from the relation HMA = GMA 2 /AMA. Now to find √ N we let b n = N/a n for all n ≥ 0. Then HMA, GMA, and AMA, respectively, reduce to Here a 0 is some positive rational number. Clearly, AMA is the same as (11). We note that the equation (a + 3/a)/2 = 1351/780 gives a = 26/15, and (a + 3/a)/2 = 26/15 holds for a = 5/3. Thus, if we employ AMA for N = 3 with a 0 = 5/3 (which is a reasonable choice, see (5)) then a 2 is the same as the upper bound of the inequality (22). We further note that the equation 6a/(a 2 + 3) = 265/153, which is the same as 265a 2 − 918a + 795 = 0 has no rational roots, and hence lower bound of (22) cannot be obtained from HMA for N = 3.
• A proof of (22) based on very simple inequalities is as follows:

Apollonius of Perga (around 262-200 BC, Greece)
He earned the title 'The Great Geometer.' Apollonius wrote a work on the cylindrical helix and another on irrational numbers, which is mentioned by Proclus Diadochus (410-485 AD, Greece).

Bakhshali Manuscript (about 200 BC)
It was found in 1881 in the village Bakhshali in Gandhara, near Peshawar, North-West India (present-day Pakistan). It is written in an old form of Sanskrit on birch bark. Only about 70 mutilated birch barks still exist, the greater portion of the manuscript has been lost. This manuscript gives various algorithms and techniques for a variety of problems, such as computing square roots, dealing with negative numbers, and finding solutions of quadratic equations. To find an approximate root of a non-square number it says "In case of a non-square (number), subtract the nearest square number; divide the remainder by twice (the root of that number). Half the square of that (that is, the fraction just obtained) is divided by the sum of the root and the fraction and subtract; (this will be the approximate value of the root) less the square (of the last term)". Thus, if N = a 2 + b, then In fact, to obtain (36) both (8) and (9) are used. Let a be the largest integer such that a 2 is less than N, and N = a 2 + b. Then, (8) gives Thus, we can use (9), to get Now let a be the smallest integer such that a 2 is greater than N, and N = a 2 − b. Then, (9) gives Thus, we can use (9) again, to get Relations (37) and (39) lead to the algorithm a n+1 = a 2 n (a 2 n + 6N) + N 2 4a n (a 2 n + N) , a 0 = a, n ≥ 0.
Clearly, in (40) we can take a any convenient real number so that a 2 is close to N. Further, from our considerations it is clear that the iterative scheme (40) converges quartically. In Table 4, we give few iterates for N = 2, 3, and 41, 105, 481 considered in Bakhshali Manuscript. An immediate extension of (4) for any nonlinear equation equation For this algorithm and its higher order extensions and their scope in real-word computation see Sen et al. [34].

Marcus Vitruvius Pollio (about 75-15 BC, Italy)
Commonly known as Vitruvius describes the use of √ 2 progression or ad quadratum technique. It uses geometry to double a square in which the diagonal of the original square is equal to the side of the resulting square.

Theon of Smyrna (about 70-135 AD, Turkey-Greece)
He described how prime numbers, geometrical numbers such as squares, progressions, music and astronomy are interrelated. He also formulated an algorithm (see Filep [35], and the Website http://numbers.computation.free.fr/Constants/Sqrt2/sqrt2.html (ac-Computation 2021, 9, 29 31 of 49 cessed on 3 March 2021)) to compute approximations of √ 2. His algorithm is based on the construction of two sequences {a n } and {b n } of natural numbers (he called a n as the side number and b n as the diagonal number), which satisfy the recurrence relations a n = a n−1 + b n−1 , b n = 2a n−1 + b n−1 , n ≥ 1. (41) We notice that and hence, if a n−1 , b n−1 is a solution of then a n , b n is a solution of b 2 − 2a 2 = ∓ 1.
Thus, it follows that b n a n = 2 ± 1 a 2 n and since lim n→∞ a n = lim n→∞ b n = ∞, we can make (1/a n ) 2 arbitrarily small. Hence, lim n→∞ (b n /a n ) = √ 2. In conclusion, if (a 0 , b 0 ) is an integer solution of (42) then (41) converges to √ 2 and the convergence is oscillatory. From these observations names for a n as the side number and for b n as the diagonal number become clear.
In the literature Equation (42) mistakenly known as Pell's equation. In fact, John Pell (1611-1685, Britain) has nothing to do with these equations. Euler mistakenly attributed to Pell a solution method that had in fact been found by William Brouncker (1620-1684, Britain), in response to a challenge by Fermat. In reality second order indeterminate equations, of the form Nx 2 + 1 = y 2 where N is an integer, were first discussed by Brahmagupta. For their solution, he employed his "Bhavana" method and showed that they have infinitely many solutions. Unfortunately, it has been recorded that Fermat was the first to assert that it has infinitely many solutions. Brahmagupta's celebrated work Brāhmasphutasiddhānta, was translated into English by Henry Thomas Colebrooke, (1765-1837, Britain). Now let (a n−1 , b n−1 ) be an integer solution of (42), then from the above observations (a n+1 , b n+1 ) = (3a n−1 + 2b n−1 , 4a n−1 + 3b n−1 ) is also a solution of the same Equation (42). Thus, if for the iterative scheme x n = 3x n−1 + 2y n−1 , y n = 4x n−1 + 3y n−1 , n ≥ 1 (43) (x 0 , y 0 ) is an integer solution of b 2 − 2a 2 = 1(−1) then (43) converges to √ 2, and the convergence will be monotonically decreasing (increasing).
It is easy to see that system (41) with (a 0 , b 0 ) = (2, 3) and (a 0 , b 0 ) = (1, 1), respectively, can be written as and a n+1 = 2a n + a n−1 , a 0 = 1, Now recall that in the construction of Table 3 for N = 2, we executed the recurrence relation (32) to obtain {s n }. It can easily be verified that a n and b n obtained from (45) are connected with s n by the relations a n−1 = s n , b n−1 = s n+1 − s n , n ≥ 2, and hence b n /a n , n ≥ 1 leads to the second column of Table 3. Similarly, a n and b n obtained from (44) are connected with s n by the relations a n−1 = s n+1 , b n−1 = s n+1 + s n , n ≥ 2, and hence b n /a n , n ≥ 1 leads to the third column of Table 3.
Similar to that of (41), system (43) with (x 0 , y 0 ) = (2, 3) and (x 0 , y 0 ) = (1, 1), respectively, can be written as and Again, looking at Table 3, we find that x n and y n obtained from (47) are connected with the same s n by the relations x n = s 2n+2 , y n = s 2n+3 − s 2n+2 , n ≥ 0, and hence y n /x n , n ≥ 0 leads to the second column of Table 3 with * , and monotonically decreasing. Similarly, x n and y n obtained from (46) are connected with s n by the relations x n = s 2n+1 , y n = s 2n+1 + s 2n , n ≥ 1, and hence y n /x n , n ≥ 1 leads to the third column of Table 3 with * , and monotonically increasing.
• For n ≥ 0, explicit solutions of the system (45) are a n = 1 Now we define T 0 = 0, T 1 = 1, T n = a n−1 b n−1 , n ≥ 2 (recall a n−1 , b n−1 , n ≥ 2, respectively, are the denominator and numerator of column 2 in Table 3) then from the above expressions it follows that which is the solution of the recurrence relation In 1778, Euler showed that are the only (infinite) numbers that are both perfect squares and triangular (t k = k(k + 1)/2). Clearly, compare to the above explicit representation of T 2 n , for the computation of T n algorithm (48) is very simple. Now to find T 2 n corresponds to which t k we need to find solutions of T 2 n = k(k + 1)/2, which is the same as finding positive integer solutions of Pell's equation b 2 − 2a 2 = 1, where b = (2k + 1) and a = 2T n . Since solutions (y k , x k ), k ≥ 2 of the system (46) computed in the second column of Table 3 with * (respectively, numerator and denominator) are first few positive integer solutions of b 2 − 2a 2 = 1, the corresponding k can be easily obtained with the relation 2k + 1 = y k . Some perfect square triangular numbers T 2 n (obtained from (48)) and the corresponding t k are as follows: · · · 1 2 6 2 35 2 204 2 1189 2 6930 2 40391 2 235416 2 · · · t 1 t 8 t 49 t 288 t 1681 t 9800 t 57121 t 332928 · · · For more details on this work see Website https://en.wikipedia.org/wiki/Square_ triangular_number (accessed on 3 March 2021).
A generalization of (41) for any integer N ≥ 2 is straightforward. In fact, for the recurrence relations a n = a n−1 + b n−1 , b n = Na n−1 + b n−1 , n ≥ 1 (49) Now since {a n } is a strictly increasing sequence, and a 2 = a 1 + b 1 = a 1 + Na 0 + b 0 > N + 1, the right side of (50) tends to zero. This means the sequence {b n /a n } converges to √ N, and the convergence is oscillatory. From (49) it also follows that b n a n − In particular, for N = 3 if we choose fundamental solution of b 2 − 3a 2 = 1 which is (a 0 , b 0 ) = (1, 2) then (49) leads to the algorithm a n = a n−1 The sequence {b n /a n } generated from (51) gives the fourth column of Table 3. We note that system (51) can be written as and its solution is Again, for N = 3 if we choose fundamental solution of b 2 − 3a 2 = −2 which is (a 0 , b 0 ) = (3, 5) then (49) leads to the algorithm a n = a n−1 The sequence {b n /a n } generated from (54) gives the fifth column of Table 3. Next, we consider the nonlinear recurrence relations a n = 2a n−1 b n−1 b n = Na 2 n−1 + b 2 n−1 , n ≥ 1 (55) and note that b 2 n − Na 2 n = (b 2 n−1 − Na 2 n−1 ) 2 = · · · = (b 2 0 − Na 2 0 ) 2 n . Thus, if (a 0 , b 0 ) is the fundamental solution (in fact, any integer solution) of b 2 − 3a 2 = 1, then the sequence {b n /a n } generated by (55) decreases monotonically to √ N. From (55), we also have b n a n − √ N = 1 2(b n−1 /a n−1 ) In Table 5, we provide first three iterates to approximate N = 2, 3, 5, and 7 with the corresponding fundamental solutions of b 2 − Na 2 = 1 as (2, 3), (1, 2), (4,9), and (3,8).  For N = 2 and 3 all entries in Table 5 are the same as in Table 3. Table 5 also indicates superiority of the nonlinear algorithm (55) compared to all linear algorithms we have discussed above. However algorithm (40) appears to have superiority. Now we will consider the recurrence relations a n = (p + q)a n−1 + 2qb n−1 b n = 2pa n−1 + (p + q)b n−1 , n ≥ 1 (56) where p = q and a 0 , b 0 are positive integers. For (56) it follows that Since a n ≥ (p + q)a n−1 implies a n ≥ (p + q) n a 0 , we find Thus, the sequence {b n /a n } generated by (56) converges to p/q, furher if b 0 /a 0 > p/q (b 0 /a 0 < p/q) the convergence is monotonically decreasing (increasing). For p = 11, q = 5 we list first few terms of {b n /a n }.

Liu Hui (around 220-280, China)
He wrote an extremely important commentary on the Jiuzhang suanshu or, as it is more commonly called, Nine Chapters on the Mathematical Art, which is believed to have been originally written around 1000 BC. This work contains approximation of π as 3.141014, and Chapter 4 Shao guang (Short width) suggest algorithms to find square and cube roots of numbers. For square roots the method is a combination of completing squares iteratively, and geometry, i.e., something like Figure 12 always in mind, see Burgos and Beltrán-Pellicer [36], Katz [14], and Yong [37]. We explain the method by considering the problem 12, where square root of 55225 is calculated. We begin with finding the integers a, b, c so that the answer can be written as 100a + 10b + c. We calculate the largest integer a so that (100a) 2 < 55225. Clearly, a = 2 is the right choice. The difference between the large (given) square (55,225) and the square with side 100a = 200, i.e., (40,000) in Figure 12 is the large gnomon with area 55,225 − 40,000 =15,225. Now if we ignore the outer thin gnomon, then b must satisfy 15225 > 2(100a)(10b) = 4000b, which gives the largest integer b = 3. To verify that the choice b = 3 is correct, i.e., when the square on 10b included, the area of the large gnomon is still less than 15,225, it is necessary to check that 2(100a)(10b) + (10b) 2 = 12900 <15,225. Since this is true, we can continue to find c. For this, we need 55,225 − 40,000-30(2×200 + 30) > 2 × 230c or 2325 > 460c. An easy check shows that the largest integer which satisfies this is c = 5. Finally, since (100a + 10b + c) 2 = (200 + 30 + 5) 2 = (235) 2 , the exact square root of 55,225 is 235.
Clearly, a = 2 is the right choice. The difference between the large (given) square (55225) and the square with side 100a = 200, i.e. (40000) in Figure 12 is the large gnomon with area 55225 − 40000 = 15225. Now if we ignore the outer thin gnomon, then b must satisfy 15225 > 2(100a)(10b) = 4000b, which gives the largest integer b = 3. To verify that the choice b = 3 is correct, i.e., when the square on 10b included, the area of the large gnomon is still less than 15225, it is necessary to check that 2(100a)(10b) + (10b) 2 = 12900 < 15225. Since this is true, we can continue to find c. For this, we need 55225 − 40000 − 30(2 × 200 + 30) > 2 × 230c or 2325 > 460c. An easy check shows that the largest integer which satisfies this is c = 5. Finally, since (100a + 10b + c) 2 = (200 + 30 + 5) 2 = (235) 2 , the exact square root of 55225 is 235.   Similar to square roots, having cubes in mind there are examples in Jiuzhang suanshu to find cube roots of numbers. For example, it is shown that the cube root of 1,860,867 is the exact number 123. In case, answer is not an exact number, the procedure continues using decimal fractions. Later Chinese extended their procedure to find roots of polynomial equations up to degree ten.

Bhaskara II or Bhaskaracharya (Working 486, India)
His contributions to mathematics include: first visual proof of the Pythagorean theorem; solutions of quadratic, cubic and quartic indeterminate equations; solutions of indeterminate quadratic equations; integer solutions of linear and quadratic indeterminate equations; a cyclic Chakravala method for solving indeterminate equations, and solutions of quadratic equations with more than one unknown, including negative and irrational solutions.

Abu Kamil, Shuja ibn Aslam ibn Muammad ibn Shuja (850-930, Egypt)
He contributed to algebra and geometry. His Book of Algebra contains a total of 69 problems. Kamil was probably the first mathematician who used irrational numbers as coefficients of an algebraic equation, and also accepted irrational numbers as solutions of the equation. In the literature often he is known as "The Reckoner from Egypt".

Abu Abd Allah Muhammad ibn Isa Al-Mahani (about 820-880, Iran-Iraq)
He wrote commentaries on parts of Euclid's Elements. In particular, for book X, Al-Mahani examined and classified quadratic irrationals and cubic irrationals. He provided definitions for rational and irrational magnitudes, which he treated as irrational numbers. He dealt with them freely but explains them in geometric terms.

Abu Ja'far al-Khazin (900-971, Iran)
He provided a meaningful definition of rational and irrational magnitudes.

Al-Hashimi (10th Century, Iraq)
He provided general proofs (rather than geometric demonstrations) for irrational numbers, as he considered multiplication, division, and other arithmetical functions. He also gave a method to prove the existence of irrational numbers.

Abu Abdallah al-Hassan ibn al-Baghdadi (10th Century, Iraq)
In his influential book Treatise on Commensurable and Incommensurable Magnitudes he related the concepts of number and magnitude by establishing a correspondence between numbers and line segments, which continues today. Given a unit magnitude a, each whole number N corresponds to an appropriate multiple Na of the unit magnitude. Parts of this magnitude, such as (p/q)a, then correspond to parts of a numbers (p/q). Al-Baghdadi considered any magnitude expressible this way as a rational magnitude. He showed that these magnitudes relate to one another as numbers to numbers. Magnitudes that are not parts he considered as irrational numbers. He also attempted to imbed the rational numbers into a number line. Al-Baghdadi also proved a result on the density of irrational magnitudes, namely, that between any two rational magnitudes there exist infinitely many irrational magnitudes. In the late nineteenth century it was proved that between any two real numbers there are infinitely many rational and irrational numbers, further irrational numbers are infinitely more numerous than rational numbers.
To see Al-Baghdadi's geometric interpretation of rational numbers, on a horizontal straight line mark two distinct points O and A, where A is right of O. Now choose the segment OA as a unit of length and let O and A represent the numbers 0 and 1,, respectively. Then the positive and negative integers can be represented by a set of points on the line spaced at unit intervals apart, the positive integers being represented to the right of O and the negative integers to the left of O. The fraction with denominator q may then be represented by the points that divide each of the unit intervals into q equal parts. Thus, each rational number can be represented by a point on the line. In Figure 13, the point P corresponds to the irrational number √ 2, which is between two rational numbers.
gave a method to prove the existence of irrational numbers.

Abu Abdallah al-Hassan ibn al-Baghdadi (10th century, Iraq). In his influential bo
Treatise on Commensurable and Incommensurable Magnitudes he related the concepts of number and ma nitude by establishing a correspondence between numbers and line segments, which continues today. Giv a unit magnitude a, each whole number N corresponds to an appropriate multiple N a of the unit mag tude. Parts of this magnitude, such as (p/q)a, then correspond to parts of a numbers (p/q). Al-Baghda considered any magnitude expressible this way as a rational magnitude. He showed that these magnitud relate to one another as numbers to numbers. Magnitudes that are not parts he considered as irration numbers. He also attempted to imbed the rational numbers into a number line. Al-Baghdadi also prov a result on the density of irrational magnitudes, namely, that between any two rational magnitudes th exist infinitely many irrational magnitudes. In the late nineteenth century it was proved that between a two real numbers there are infinitely many rational and irrational numbers, further irrational numbers a infinitely more numerous than rational numbers. To see Al-Baghdadi's geometric interpretation of rational numbers, on a horizontal straight line mark t distinct points O and A, where A is right of O. Now choose the segment OA as a unit of length and let O a A represent the numbers 0 and 1, respectively. Then the positive and negative integers can be represented a set of points on the line spaced at unit intervals apart, the positive integers being represented to the rig of O and the negative integers to the left of O. The fraction with denominator q may then be represent by the points that divide each of the unit intervals into q equal parts. Thus, each rational number can represented by a point on the line. In Figure 13, the point P corresponds to the irrational number √ 2, wh is between two rational numbers.

Omar Khayyám (1048-1131, Iran)
He is considered one of the major mathematicians and astronomers of the medieval period. His major contributions include the length of the year 365.24219858156 days, commentary on Euclid's Elements, Euclid's parallel postulate, and his classification to nineteen types of cubic equations. He believed that for cubic equations arithmetic solutions were impossible. To the Western world Omar is known as the author of The Rubaiyat (Persian poetry). Omar considered the problems of irrational numbers and their relations to rational numbers. He called irrational magnitudes as numbers themselves. He writes that methods for calculating square and cube roots came from India, and he has extended them to the determination of roots of any order.

Nilakanthan Somayaji (around 1444-1544, India)
His work includes the inductive mathematical proofs, a derivation and proof of the arctangent trigonometric function, and in Sanskrit poetry the series In the literature (5) is known as Gregory-Leibniz after James Gregory (1638-1675, Scotland) and Gottfried Wilhelm von Leibniz (1646-1716, Germany) series. He also gave sophisticated explanations of the irrationality of π, the correct formulation for the equation of the center of the planets, and a heliocentric model of the solar system. The following expansion of π is also due to him This series converges faster than (57).

Nicolas Chuquet (around 1445-1488, France)
He is famous for his work Triparty en la science des nombres, which was not published in his lifetime. In this work, he introduced our familiar numerical terms billion, trillion, quadrillion, etc. He also showed that √ 5 2 161 682 and √ 6 2 89 198 .

Michael Stifel (1486-1567, Germany)
He invented logarithms independently of Napier, using a totally different approach. His most famous work is Arithmetica Integra which was published in 1554. For irrational numbers, he wrote: "We are moved and compelled to assert that they are numbers, compelled that is, by the results which follow from their use. On the other hand ... just as an infinite number is not a number, so an irrational number is not a true number, but lies hidden in some sort of cloud of infinity". He and at the same time Johannes Scheubel (1494-1570, Germany) gave a method to find higher order roots using appropriate row of the Blaise Pascal (1623-1662, France) triangle, a method which Chinese have already used several centuries earlier to extract roots of polynomial equations up to degree ten.

Guillaume Gosselin (1536-1600, France)
He summarized the rules that supported the conventional arithmetic, algebraic calculations on geometric progressions, extraction of roots, calculations of irrational expressions, and notations of the objects of algebra, then the rules that resolved equations of first and second degrees with one unknown with numerical coefficients.

Zhu Zaiyu (1536-1611, China)
In 1604 wrote a New Explanation of the Theory of Calculation in which he derived values of the roots of 2. He was so attracted to √ 2 that he used nine abacuses to compute it to 25-digit accuracy!

Francois Viéte (1540-1603, France)
He has been called the father of modern algebra and the foremost mathematician of the sixteenth century. In his 1593 book Supplementum geometriae, he showed 3.1415926535 < π < 3.1415926537, i.e., gave the value of π correct to 9 places. He also represented π as an infinite product 2 π = cos π 4 cos π 8 cos π 16 cos π 32 · · · = √ 2 2 This formula is one of the milestones in the history of π. The convergence of Vieta's formula was proved by Ferdinand Rudio (1856-1929, Germany) in 1891.

Simon Stevin (1548-1620, Elgium)
In 1585, he published a 36-page booklet, La Thiende (The Tenth), which contains a brief account of decimal fractions. He showed that it is possible to add and multiply infinite decimals, so that infinite decimals not only look like numbers, but also behave like them. He compared fractions with an unknown island having beautiful fruits, pleasant plains, and precious minerals. In the defence of irrational numbers, he objected to the very idea of calling them "irrational", because all numbers are equally concrete from a geometric point of view.

John Wallis (1616-1703, Britain)
Wallis was the most influential English mathematician before Newton. In 1662, he led the formation of the Royal Society of London, in 1655 he published a treatise Arithmetica infinitorum that defined conic sections analytically, in Calculus he extended the works of all his predecessors, and presented π as an infinite product · · · (60) without proof. The convergence of Wallis formula immediately follows by using squeeze theorem of calculus. On the night of 22 December 1669, he occupied himself with finding the integral part of the square root of 3 × 10 40 while in bed, and several hours afterward he wrote down the result from memory. Two months later, he was challenged to extract the square root of a number of 53 digits, which he performed mentally; one month later, he dictated the answer that he had not committed to writing.

Jacob Bernoulli: Introduction of e
He was first of the eight prominent mathematicians in the Bernoulli family. Jacob wrote on infinite series, studied many special curves, invented polar coordinates, developed logarithmic differentiation, and introduced the Bernoulli numbers. he formulated the basic principle in the theory of probability known as Bernoulli's theorem or the law of large numbers. In 1683, he introduced the number e.

Continued Fractions
There is a sufficient evidence that Aryabhata used continued fraction to solve a linear indeterminate equation. A simple manipulation of Euclid's Algorithm which is mainly used to find gcd of the ratio of two geometric magnitudes leads to a finite (infinite) continued fraction provided the ratio is rational (irrational and pointed out how the partial fractions are successively larger and smaller than , and the process converges to . While Brouncker was not kind enough to provide details of his expansion, Wallis in his book Opera Mathematica of 1695 detailed basic facts and properties of continued fractions (this term is also coined by him). Later Huygens, Euler, Lambert, and Joseph Louis Lagrange (1736-1813, France-Italy) enriched the theory and applications of continued fractions to the extend that it became a subject in its own right. Especially, Euler showed that every rational number can be expressed as a terminating simple continued fraction, and consequently every infinite continued fraction is irrational. Continued fractions play dominate role in finding best rational approximations of irrational numbers. An equation of the type which Alfred Pringsheim (1850-1941, Germany) wrote as and when all a i = 1, i ≥ 1 Gauss wrote as is called the k-th convergent r k of a continued fraction, where a i , b i are integers (functions) determined from a given irrational number (function). It is well known, for example, see Agarwal [38], and Sen and Agarwal [13] that r k can be written as r k = u k /v k , where the numerator u k and the denominator v k satisfy the recurrence relations Now we shall use the algorithm (61) to find rational approximations of a given positive number √ N. For this, again we assume that a to be an initial guess of √ N so that N = a 2 + r.
Thus, (61) reduces to which is exactly the same as (45) with u k = b k , v k = a k , k ≥ 1 and therefore {r k = u k /v k } k=1 leads to the second column of Table 3, and forms the sequence of best approximations of √ 2. Similarly, for N = 3 and a = 1, we have r = 2 and Thus, (61) reduces to whose solution apperas as From (53) and (62) it follows that u k = 2b k−1 , v k = 2a k−1 , k ≥ 1 and therefore {r k = u k /v k } k=1 gives the fourth column of Table 3.

Leonhard Euler: Irrationality of e
He was probably the most prolific mathematician who ever lived. Euler's energy and capacity for work were virtually boundless. His collected works form about 80 quarto sized volumes and it is believed that much of his work has been lost. What is particularly astonishing is that Euler became virtually sightless in his right eye during the mid-1730s, and was blind for the last 17 years of his life, and this was one of the most productive periods! In 1748, Euler considered the functions cos x and sin x/x, x = 0 which have roots at (2n ± 1)π/2, n ≥ 0 and ±nπ, n ≥ 1, respectively. Then, he wrote these functions in terms of infinite products and Substituting x = π/4 in (63), we get the following infinite product of 1/ √ 2, which gives √ 2 = 2 · 2 1 · 3 6 · 6 5 · 7 10 · 10 9 · 11 · · · = ∏ k≥0 (4k + 2) 2 (4k + 1)(4k + 3) .
It is interesting to note that (64) with x = π/2 immediately gives Wallis' formula (47). Now in (64) on equating the coefficients of x 2 , we get which is the same as 1 The above demonstration of Euler is based on manipulations that were not justified at the time, and it was not until 1741 that he was able to produce a truly rigorous proof. In the literature summing the series in (66) has become famous as Basel problem and currently for this several proofs are known. The convergence of the infinite products (59), (60), (65), and series (66) are extremely slow, but they are immensely of theoretical interest.
However, p(q − 1)! − q! ∑ q k=0 1/k! is a positive integer. Thus, e is irrational. In 1998, Martin Aigner (born 1942, Austria) and Günter Matthias Ziegler (born 1963, Germany) provided details about the irrationality of e x for any nonzero rational x. We also remark that employing several different algorithms, a massive details about the approximations of e have been given in Sen and Agarwal [13].

Johann Heinrich Lambert: Irrationality of π
Lambert was the first to introduce hyperbolic functions into trigonometry. He generalized Euler's method to show that continued fractions of e x and tan x are irrational if x is a nonzero rational. Lambert's following continued fractions of e x and tan x of 1761 are of great historical importance To prove the irrationality of π, in 1768, Lambert substituted x = π/4 in (68), so that the left side of (68) is simply one. Then he assumed that there exist integers p and q such that π/4 = p/q, i.e., π/4 is rational and then showed that the right side of (68) is irrational. The complete Lambert's proof is available on the Website https://math.stackexchange. com/questions/895611/lamberts-original-proof-that-pi-is-irrational (accessed on 3 March 2021). After Lambert's proof several prominent mathematicians gave alternative proofs (claimed to be simpler) to prove the irrationality of π, for example, in 1794, Charles Hermite (1822-1901, France) showed that π 2 is irrational, from which the irrationality of π follows immediately; in 1945, Dame Mary Lucy Cartwright (1900-1998, England) set as an example in an exam at the Cambridge University a new proof of the irrationality of π (the origin of the proof is not yet known); this is followed by the proof of Nicolas Bourbaki (born 1939, France) in 1949; then Ivan Morton Niven (1915Niven ( -1999, Canadian-American) in 1947 (also see his book [39]), and Miklós Laczkovich (born 1948, Hungary) in 1997. For details see the Website https://en.wikipedia.org/wiki/Proof_that_piis_irrational (accessed on 3 March 2021).
For each natural number q and each nonnegative integer n, let A n (q) = q n π 0 x n (π − x) n n! sin(x)dx.
Since A n (q) is the integral of a function which is defined on [0, π], takes the value 0 at the lower and upper limits and positive in (0, π), A n (q) > 0. Further, since x(π − x) ≤ (π/2) 2 , we have A n (q) ≤ πq n 1 n! π 2 2n = π (qπ 2 /4) n n! and hence A n (q) < 1 for sufficiently large n. On the other hand, recursive integration by parts leads to the fact that, if p and q are natural numbers such that π = p/q and f is the polynomial function from [0, π] to IR defined by f (x) = x n (p − qx) n /n!, then Since f is a polynomial of degree 2n, the last term is 0. Now since each function f (k) , 0 ≤ k ≤ 2n as well as sin(x) and cos(x) take integer values at 0 and π, this shows that A n (q) is an integer. Since it is greater than 0, it must be a positive integer. However, we have seen that A n (q) < 1 if n is sufficiently large. This contradiction shows that π = p/q is impossible.

Joseph Liouville (1809-1882, France)
The fascinating and difficult theory of transcendental numbers originated in the work of Liouville. The irrationality of e and π, which is equivalent to the fact these numbers are not roots of any linear equation of the form ax + b = 0 whose coefficients are integers, had been proved by Euler and Lambert. In 1844 Liouville showed that e is also not a root of any quadratic equation with integral coefficients. This led them to conjecture that e does not satisfy any polynomial equation a n x n + a n−1 x n−1 + · · · + a 1 x + a 0 = 0 with integral coefficients. However, all his efforts to prove this conjecture failed.

Karl Theodor Wilhelm Weierstrass: Sequential Definition of Irrationality
He has been acknowledged as the "father of modern analysis". In 1872 his rigorous work led them to discover a function that, although continuous, had no derivative at any point. Following the work of Augustin-Louis Cauchy (1789-1857, France) of 1821, Weierstrass continued the sequential definition of irrational numbers based on convergent series. This definition was further extended to classes of equivalent sequences by Heinrich Edward Heine (1821-1881, Germany) in 1872.
In Ramanujan's first letter on 31 January 1913 to Godfrey Harold Hardy (1877-1947, England), among several theorems the above series was one of the results.

Charles Hermite: Transcendence of e
He made two very significant contributions to mathematics. In 1858, acknowledging Niels Henrik Abel's (1802-1829, Norway) result that the general polynomial equation of the fifth degree cannot be solved by functions involving only rational operations and root extractions, Hermite showed that such equations can be solved by elliptic functions. In 1873, he proved the conjecture of Liouville and affirmed that e is indeed a transcendental number. According to him, we are servants rather than masters in mathematics. Here we shall follow an elegant demonstration of Richard Schwartz (born 1966, USA) to prove the

Aleksander Osipovich Gelfond (1906-1968, Russia)
In 1934, he gave an affirmative answer to Hilbert's seventh problem, and this was followed by Theodor Schneider (1911( -1988. This result in the literature is known as Gelfond-Schneider theorem. This result does not help determine whether numbers such as e e , π π , or π e are transcendental, since both the bases and exponents are transcendental numbers and therefore do not satisfy the conditions of the Gelfond-Schneider theorem. Log 2 (base 10) can be shown to be transcendental using the Gelfond-Schneider theorem. The transcendental number 2 √ 2 is known as the Gelfond-Schneider constant (or Hilbert number), and the transcendental number e π = (e iπ ) −i = (−1) −i is known as Gelfond's constant.
In 1979 Apéry published an unexpected proof of the irrationality of ζ(3). In the literature ζ(3) is known as Apéry constant. The irrationality of ζ(2n + 1), n ≥ 2 is expected but not yet established.

More on Numbers
The trigonometric functions sin x, cos x, tan x, csc x, sec x, cot x, and their hyperbolic counterparts, for any nonzero algebraic number x, expressed in radians are transcendental. The number e π √ n is also transcendental. It is not yet known if the numbers π + e, π − e, πe, π e . π π are rational, algebraic, irrational, or transcendental. However, it is certain that both π + e and πe cannot be rational (or algebraic). In fact, if both are rational then (π + e) 2 − 4πe is rational. However, this gives (π − e) 2 , and so π − e is algebraic. However, then adding and subtracting π − e with π + e, we find that both π and e are algebraic, which contradicts the fact that both are transcendental.

Conclusions
One of the greatest discovery in whole of mathematics is the invention of irrational numbers, and then their understanding. In this article we have demonstrated that Vedic Ascetics more than 5000 years back were unsuccessful in finding exact values of the numbers √ 2 and π. The ancient records (supported by great philosophers, mathematicians, and historians) stipulate that Vedic Ascetics were also definite that these numbers are incommensurable/irrational. We have exhibited that the claim of the historians of mathematics that Pythagoras proved the irrationality of √ 2 is only conjectural. In fact, the first geometric proof of the irrationality of √ 2 appeared only in Meno (Socratic dialogue by Plato) almost two hundred years after Pythagoras. Since then several different proofs of the irrationally of √ 2 and in general for √ N for any natural number N which is not not a perfect square have been given. We have provided some of these important proofs. The next major understanding of irrational numbers came from the scholars of the Islamic Middle East towards the end of the first millennium CE. They started treating irrational numbers as algebraic objects, and most importantly provided a geometric interpretation of rational numbers on a horizontal straight line. Since then research continues for the known as well as unknown/expected irrational numbers, their subset of transitive numbers, and their computation to trillions of decimal places, we have detailed some of these advancements. We have also discussed Dedekind-Cantor axiom of the 19th century which provides geometric interpretation of all real numbers, and thus completes the Islamic work. We have arranged individuals contributions chronologically to show that each continent of the world has contributed in this fascinating field of mathematics.