|   
  
 Gottfried Wilhelm Leibniz
 (1646 - 1716)
 
 
 Why always do complicated?
 
	
   (simply !) 
 and a derived sequnce:
 
 
 Parts of his life:
 
	
  Leibniz was born in Leipzig in 1646. 
    Although German mathematician and philosopher, all his work was written in 
    French or Latin. Son of a philosophy professor, he received very little mathematical 
    teaching, but he discovered a hidden interest for mathematics thanks to Huygens 
    with who he forms a friendship... (Thanks Huygens!). Leibniz was especially interested in numerical 
    series - That is interesting ! - and claims to be the father of the differential 
    calculus, which will lead him to a violent dispute with Newton and darken 
    the remaining of his life. We do not always know this, but this good 
    old Leibniz was a great notation inventor. we owe him the  (integral) sign, = , dx , . for the multiplication and division. But Leibniz, although he was a genius, did 
    not have a neat idea about complex numbers for example, and wrote lines and 
    lines of equation without any real meaning, like: 
 with which he astonished Huygens. His philosophic conceptions influenced a lot 
    his mathematical vision of things... He considered for example these complex 
    numbers as between existence and non existence! Ignored at his death by the German and British 
    scientific communities, the famous French hundred-year-old Fontenelle actually 
    praised him in Paris... He was worth it! Around   
  Here I am exaggerating a bit, because Leibniz 
    was not the real discoverer of this formula. James Gregory (1638-1675) 
    had in fact calculated the whole sequence expansion: arctan(x)= for x between-1 and 1. It would 
    be astonishing if Grégory had not seen the particular case x=1 which gives the following formula: 
 But here, as we will see, it's convergence 
    is more than execrable and Grégory must have realized how useless this 
    formula is practically... But admit that it is beautiful, and extremely simple! 
    So let us give Grégory what he deserves, even if the first publication 
    of this formula is only explicitly shown in Leibniz's works...On the other hand, this formula from Leibniz/Grégory reserves a few 
    surprises, discovered only recently by Roy North... Go and see about them 
    in applications !
 Demonstration 
 
	
  What can we say about this demonstration, 
    it seem evident... We can remind ourselves that we have arctan:: R->]- /2,  /2 is defined as the inverse function of the bijection tan : ]-  /2,  /2[->R. But we have tan'(x)=1+tan2(x) (using cos and sin ) and hence, 
    from the formula of the derivative of an inverse function . But we also know the whole sequence
  defined on]-1,1[. We know that we can integrate it term by term on 
    this interval from the properties of whole sequences. After seeing that the 
    result  converges at 1 and -1 which insures the uniform convergence 
  of the sequence over[-1,1] towards arctan, we take x=1 and that's it!... We can also use x=1/31/2, and that is more interesting: 
 Applications 
 
	
  Behold, this is a great moment!For the Leibniz/Grégory formula:
 
 
 
	
		| n=10 | 3,2323
			(0) |  
		| n=100 | 3,151493401 (1) |  
		| n=1 000 | 3,14259165 (2) |  
		| n=10 000 | 3,1416926435905432 (3) |  
		| n=100 000 | 3,14160265348979398846014336 (3) |  
		| n=500 000 | 3,141590653589793240462643383269502884197 (5) |  
		| n=1 000 000 | 3,14159365358879323921264313 (5) |  If there are a few digits in brown, it is because they are incorrect! But then, 
  how come the following digits be right? Is the convergence not logarithmic from 
  its form (approx. log(n)) ? That is the problem that Roy North asked the Borwein brothers a few years ago... I have not got the solution, 
  but it involves Euler 's numbers and two 
  summing formulae. Still waiting for further answers...
 
 For the second formula (just over the Application section), we have an interesting 
  linear convergence of about n/2 :
 
 
 
	
		| n=10 | 3,1415933045 |  
		| n=100 | 49 digits correct |  Acceleration of the convergence : 
  
  If there is a formula for which Aitken's Delta2 
    is really really useful, then it is Leibniz's... And other acceleration methods 
    work as well! (a logarithmic convergence is not very convincing)...
 1) With 
    the Delta2 : classical acceleration, and faster and faster normally... requires a computation 
    with a lot of digits because of the unstable numerical side of the process.
 
 2) With 
    an average:
 As it is a modifies sequence, I had the idea of using of applying an average 
    value to Leibniz's series. What should of brought us in theory 1 digit in 
    the best case seems to accelerate the convergence in a quite surprising way!:
 let
  be the la moyenne pondérée associée. 
 The results appear in the table below for the 2 
    types of acceleration:
 
 
	
		|  | Leibniz | Delta2 | Average value |  
		| n=3 | 2,8952 (0) | 3,13333
			(1) | 3,161904
			(1) |  
		| n=5 | 2,976 (0) | 3,13968
			(1) | 3,143434
			(2) |  
		| n=10 | 3,2323
			(0) | 3,141839
			(3) | 3,14150053
			(4) |  
		| n=100 | 3,15149
			(1) | 3,141592905
			(6) | 3,14159264593
			(7) |  
		| n=1000 | 3,142591
			(2) | 3,141592653839792
			(9) | 3,141592653589041
			(12) |  
		| n=10000 | 3,14169264
			(3) | 3,1415926535900
			? (10) | 3,1415926535897931634
			(15) |  The last result for Delta2 seems suspicious to me, according to the average 
  value found, although I calculated it with 100 digits...
 
 One good thing about series made up from others, is that we can repeat the process.
 Immediate application!
 
 
	
		|  | Delta2 repeated 2 times | Average repeated 2 times |  
		| n=3 | 3,13888 (1) | does not exist... |  
		| n=5 | 3,141450 (2) | 3,1215 (1) |  
		| n=10 | 3,141595655 (5) | 3,141598653 (5) |  
		| n=100 | 3,1415926536094 (9) | 3,141592653589922 (12) |  
		| n=1000 | 3,1415926535897934269 (15) | 3,14159265358979323847405 (19) |  
		| n=10000 | too suspicious! | 3,14159265358979323846264338440 (26) |  
 Just note that combinations between Delta2 and means 
  are less efficient.
 3) we can also rearrange the terms as Leibniz had done, but it is not much more 
  efficient.
 
 All right, this time I think that's all concerning this sequence that proved 
    in the end to be quite rich!!!
 
 back to home page
 |