CINXE.COM
MATHEMATICA TUTORIAL, Part 2.5; Biography of Joseph Fourier
<html xmlns="http://www.w3.org/1999/xhtml"><head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8" /> <script type="text/javascript" src="//cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"> </script> <meta http-equiv="Content-Type" content="text/html; charset=UTF-8"> <title>MATHEMATICA TUTORIAL, Part 2.5; Biography of Joseph Fourier</title> <script type="text/javascript" src="../js/jquery-1.12.4.min.js"></script> <script type="text/javascript" src="../js/bootstrap.min.js"></script> <script type="text/javascript" src="../js/main.js"></script> <link href="../css/bootstrap.min.css" rel="stylesheet"> <link type="text/css" rel="stylesheet" href="../css/style.css" /> <link type="text/css" rel="stylesheet" href="../css/am34.css" /> <!--- <script type="text/javascript"> $(document).ready(function () { document.getElementById('show').style.zIndex = '1'; }); function clickcaption() { if (document.getElementById('show').style.zIndex == '1') { $('#subshow').hide('slow'); document.getElementById('show').style.zIndex = '0'; document.getElementById('show').innerHTML = '[show]'; } else { $('#subshow').show('slow'); document.getElementById('show').style.zIndex = '1'; document.getElementById('show').innerHTML = '[hide]'; } } function doit(j) { $('.chapter').css('display','none'); var b = $('.chapter'); b[j].style.display = 'block'; } function clickit(z) { switch (z) { case 0: {$('#intro').slideToggle();} break; case 1: {$('#first').slideToggle();} break; case 2: {$('#second').slideToggle();} break; case 3: {$('#third').slideToggle();} break; case 4: {$('#forth').slideToggle();} break; case 5: {$('#fifth').slideToggle();} break; case 6: {$('#sixth').slideToggle();} break; case 7: {$('#seventh').slideToggle();} break; case 8: {$('#eighth').slideToggle();} break; case 9: {$('#ninth').slideToggle();} break; case 10: {$('#tenth').slideToggle();} break; case 11: {$('#eleventh').slideToggle();} break; case 12: {$('#twelfth').slideToggle();} break; case 13: {$('#thirteen').slideToggle();} break; case 14: {$('#fourteen').slideToggle();} break; case 15: {$('#fifteen').slideToggle();} break; case 16: {$('#sixteen').slideToggle();} break; case 17: {$('#seventeen').slideToggle();} break; case 18: {$('#eighteen').slideToggle();} break; } } </script> ----> <script type="text/javascript"> tday=new Array("Sunday","Monday","Tuesday","Wednesday","Thursday","Friday","Saturday"); tmonth=new Array("January","February","March","April","May","June","July","August","September","October","November","December"); function GetClock(){ var d=new Date(); var nday=d.getDay(),nmonth=d.getMonth(),ndate=d.getDate(),nyear=d.getFullYear(); var nhour=d.getHours(),nmin=d.getMinutes(),nsec=d.getSeconds(),ap; if(nhour==0){ap=" AM";nhour=12;} else if(nhour<12){ap=" AM";} else if(nhour==12){ap=" PM";} else if(nhour>12){ap=" PM";nhour-=12;} if(nmin<=9) nmin="0"+nmin; if(nsec<=9) nsec="0"+nsec; document.getElementById('clockbox').innerHTML=""+tday[nday]+", "+tmonth[nmonth]+" "+ndate+", "+nyear+" "+nhour+":"+nmin+":"+nsec+ap+""; } window.onload=function(){ GetClock(); setInterval(GetClock,1000); } </script> <style> .dropdown { position: relative; display: inline-block; } .dropdown-content { display: none; position: absolute; background-color: #f9f9f9; min-width: 260px; box-shadow: 0px 8px 16px 0px rgba(0,0,0,0.2); padding: 12px 16px; z-index: 1; } .dropdown:hover .dropdown-content { display: block; } </style> <style type="text/css"> .custom { font-family: Times; font-weight: bold; color: brown; font-size:150%; } </style> </head> <body> <div id="header"> <h1>MATHEMATICA <div class="dropdown"><span>TUTORIAL </span> <div class="dropdown-content"> <p>under the terms of the GNU General Public License (<a href="https://en.wikipedia.org/wiki/GNU_General_Public_License">GPL</a>)</p> </div> </div> for the Second Course. Part V: Biography of Joseph Fourier </h1> <div id="authors"> <h3> Email: <a href="mailto:Vladimir_Dobrushkin@brown.edu?Subject=Mathematica" style="font-size: 14px; color: black;">Prof. Vladimir Dobrushkin</a> (<b><span id="clockbox" style="font-family:Times; font-size:1em; color:brown;"></span></b>)</h3> </div> <hr color="#cccccc" noshade="noshade" size="1" width="100%"> </div> <div id="content"> <div class="chapter"> <!---- <h1>Preface</h1> <hr color="#cccccc" noshade="noshade" size="1" width="100%"> <p> This tutorial was made solely for the purpose of education and it was designed for students taking Applied Math 0340. It is primarily for students who have some experience using <em>Mathematica</em>. If you have never used <em>Mathematica</em> before and would like to learn more of the basics for this <a href="https://en.wikipedia.org/wiki/Computer_algebra_system">computer algebra system</a>, it is strongly recommended looking at the APMA 0330 <a href="http://www.cfm.brown.edu/people/dobrush/am33/Mathematica/index.html">tutorial</a>. As a friendly reminder, don't forget to clear variables in use and/or the kernel. The <em>Mathematica</em> commands in this tutorial are all written in <b>bold black font</b>, while <em>Mathematica</em> output is in regular fonts. </p> <p>Finally, you can copy and paste all commands into your <em>Mathematica</em> notebook, change the parameters, and run them because the tutorial is under the terms of the GNU General Public License (<a href="https://en.wikipedia.org/wiki/GNU_General_Public_License">GPL</a>). You, as the user, are free to use the scripts for your needs to learn the <em>Mathematica</em> program, and have the right to distribute and refer to this tutorial, as long as this tutorial is accredited appropriately. The tutorial accompanies the textbook<a href="https://www.crcpress.com/Applied-Differential-Equations-The-Primary-Course/Dobrushkin/9781439851043"> <em>Applied Differential Equations. The Primary Course</em> </a> by Vladimir Dobrushkin, CRC Press, 2015; <a href="http://www.crcpress.com/product/isbn/9781439851043">http://www.crcpress.com/product/isbn/9781439851043</a></p> ----> </div> <div id="menu"> <h4>Contents <a href="javascript:clickcaption();" id="show">[hide]</a></h4> <div id="subshow"> <ul> <li><a href="http://www.cfm.brown.edu/people/dobrush/am34/Mathematica/index.html">Preface</a> <ol class="submenu" start="0" type="1"> <li><a href="javascript:clickit(0)">Introduction</a> <ol id="intro" start="1" type="i"> <li><a href="../intro/3dplot.html">3D plotting</a></li> <li><a href="../intro/tube.html">Tubing</a></li> <li><a href="../intro/exist.html">Existence and Uniqueness </a></li> <li><a href="../intro/picard.html">Picard iterations </a></li> <li><a href="../intro/adomian.html">Adomian iterations </a></li> <li><a href="../intro/lagrange.html">Euler--Lagrange equations </a></li> <li><a href="../intro/multiply.html">Lagrange multipliers </a></li> <li><a href="../intro/hamilton.html">Hamilton principle </a></li> <li><a class="subli" href="../intro/convert.html">Converting to a system </a></li> <li><a class="subli" href="../intro/scale.html">Scaling </a></li> <li><a class="subli" href="http://www.dobrush.us/probability/index.php">Probability </a></li> </ol> </li> <li><a href="../ch1/part1.html">Part I:</a> <a href="javascript:clickit(1)">Matrix Algebra</a> <ol id="first" start="1" type="i"> <li><a class="subli" href="../ch1/vector.html">How to define vectors </a></li> <li><a class="subli" href="../ch1/matrix.html">How to define matrices </a></li> <li><a class="subli" href="../ch1/basic.html">Basic operations with matrices </a></li> <li><a class="subli" href="../ch1/solve.html">Linear systems of equations </a></li> <li><a class="subli" href="../ch1/det.html">Determinants and inverses </a></li> <li><a class="subli" href="../ch1/special.html">Special matrices </a></li> <li><a class="subli" href="../ch1/eigen.html">Eigenvalues and Eigenvectors </a></li> <!--- <li><a class="subli" href="ch1/geneigen.html">Generalized -- --Eigenvectors </a></li> --> <li><a class="subli" href="../ch1/diagonalization.html">Diagonalization procedure </a></li> <li><a class="subli" href="../ch1/sylvester.html">Sylvester formula </a></li> <li><a class="subli" href="../ch1/resolvent.html">The resolvent method </a></li> <li><a class="subli" href="../ch1/spectral.html">Polynomial interpolation </a></li> <li><a class="subli" href="../ch1/positive.html">Positive matrices </a></li> <li><a class="subli" href="c../h1/roots.html">Roots </a></li> <li><a class="subli" href="c../h1/application1.html">Miscellany</a></li> </ol> </li> <li><a href="../ch2/part2.html">Part II:</a> <a href="javascript:clickit(2)">Linear Systems of Ordinary Differential Equations</a> <ol id="second" start="1" type="i"> <li><a class="subli" href="../ch2/motiv.html">Motivation </a></li> <li><a class="subli" href="../ch2/variable.html">Variable coefficient systems of ODEs</a></li> <li><a class="subli" href="../ch2/floquet.html">Floquet theory</a></li> <li><a class="subli" href="../ch2/constant.html">Constant coefficient systems of ODEs</a></li> <li><a href="../ch2/portrait.html">Planar Phase Portrait</a></li> <li><a href="../ch2/eulersys.html">Euler systems of equations</a></li> <li><a class="subli" href="../ch2/exp.html">Fundamental matrices</a></li> <li><a class="subli" href="../ch2/reduction.html">Reduction to a single equation </a></li> <li><a class="subli" href="../ch2/undetcoeff.html">Method of undetermined coefficients</a></li> <li><a class="subli" href="../ch2/varofparam.html">Variation of parameters</a></li> <li><a class="subli" href="../ch2/laplace.html">Laplace transform</a> </li> <li><a class="subli" href="../ch2/second.html">Second order ODEs</a></li> <li><a class="subli" href="../ch2/spring.html">Spring-mass systems</a></li> <li><a class="subli" href="../ch2/circuit.html">Electric circuits</a></li> <li><a class="subli" href="../ch2/application.html">Applications</a></li> </ol></li> <li><a href="../ch3/part3.html">Part III:</a> <a href="javascript:clickit(3)">Non-linear Systems of Ordinary Differential Equations</a> <ol start="1" id="third" type="i"> <li><a class="subli" href="../ch3/ndsolve.html">Numerical solutions</a></li> <li><a class="subli" href="../ch3/stability.html">Stability </a></li> <li><a class="subli" href="../ch3/linear.html">Linearization</a></li> <li><a class="subli" href="../ch3/planar.html">Planar autonomous systems</a></li> <li><a class="subli" href="../ch3/conserv.html">Conservative systems</a></li> <li><a class="subli" href="../ch3/grad.html">Gradient systems</a></li> <li><a class="subli" href="../ch3/species.html">Competing species</a></li> <li><a class="subli" href="../ch3/lotka.html">Predator-Prey equations</a></li> <li><a class="subli" href="../ch3/harvest.html">Harvesting species</a></li> <li><a class="subli" href="../ch3/lyapunov.html">Lyapunov second method</a></li> <li><a class="subli" href="../ch3/hiv.html">HIV models</a></li> <li><a class="subli" href="../ch3/period.html">Periodic solutions </a></li> <li><a class="subli" href="../ch3/asynch.html">Asynchronous solutions </a></li> <li><a class="subli" href="../ch3/bendixson.html">Limited Cycles </a></li> <li><a class="subli" href="../ch3/pol.html">van der Pol equations</a></li> <li><a class="subli" href="../ch3/neuron.html">Neuroscience </a></li> <li><a class="subli" href="../ch3/enzyme.html">Biochemistry </a></li> <li><a class="subli" href="../ch3/application3.html">Miscellany </a></li> </ol> </li> <a href="../ch3/chaos.html" style="color: #0A0A3C"><b>Part III C:</b></a> <a href="javascript:clickit(13)">Chaos</a> <ol start="1" id="thirteen" type="i"> <li><a class="subli" href="../ch3/lorenz.html">Lorenz equations </a></li> <li><a class="subli" href="../ch3/rossler.html">Rössler attractor </a></li> <li><a class="subli" href="../ch3/electric.html">Electric circuits </a></li> <li><a class="subli" href="../ch3/chua.html">Chua circuits </a></li> <li> <a class="navSectionTitle" href="../ch3/pendulum.html">Pendulum</a> <a class="showDropdown" href="javascript:clickit(8)"> <b class="caret"></b></a> <ol start="1" id="eighth" type="a"> <!-- <div class="navSection"> <a class="navSectionTitle" href="ch3/pendulum.html">Pendulum</a> <a class="showDropdown" href="javascript:clickit(8)"><b class="caret"></b></a> </div> <ol start="1" id="eighth" type="a"> --> <li><a class="subli" href="../ch3/ppendulum.html">Phase Portrait</a></li> <li><a class="subli" href="../ch3/mpendulum.html">Moving pivot</a></li> <li><a class="subli" href="../ch3/npendulum.html">Numerical simulation</a></li> <li><a class="subli" href="../ch3/spendulum.html">Spring pendulum</a></li> <li><a class="subli" href="../ch3/ependulum.html">Elastic pendulum</a></li> <li><a class="subli" href="../ch3/dpendulum.html">Double pendulum</a></li> <li><a class="subli" href="../ch3/escapement.html">Escapement</a></li> <li><a class="subli" href="../ch3/ipendulum.html">Inverted pendulum</a></li> <li><a class="subli" href="../ch3/adomian.html">ADM approximation</a></li> <li><a class="subli" href="../ch3/fourier.html">Fourier series</a></li> </ol> </li> <li> <a class="navSectionTitle" href="../ch3/mass.html">Spring-mass system</a> <a class="showDropdown" href="javascript:clickit(9)"> <b class="caret"></b></a> <ol start="1" id="ninth" type="a"> <li><a class="subli" href="../ch3/anharmonic.html">Anharmonic motion</a></li> <li><a class="subli" href="../ch3/fnharmonic.html">Forced anharmonic motion</a></li> <li><a class="subli" href="../ch3/duffing.html">Duffing equations</a></li> <li><a class="subli" href="../ch3/fduffing.html">Forced Duffing equation</a></li> <li><a class="subli" href="../ch3/hard.html">Driven hard spring</a></li> <li><a class="subli" href="../ch3/soft.html">Driven soft spring </a></li> <li><a class="subli" href="../ch3/coulomb.html">Coulomb damping</a></li> <li><a class="subli" href="../ch3/quadratic.html">Quadratic damping</a></li> </ol> </li> <li><a class="subli" href="../ch3/mechanic.html">Mechanical problems </a></li> <li><a class="subli" href="../ch3/miscellany.html">Miscellany </a></li> </ol> </li> <li><a href="../ch4/part4.html">Part IV:</a> <a href="javascript:clickit(4)">Numerical Methods</a> <ol start="1" id="forth" type="i"> <li><a class="subli" href="../ch4/seidel.html">Iterative methods </a></li> <!---- <li><a class="subli" href="../ch4/newton.html">Picard Iterations </a></li> ----> <li><a class="subli" href="../ch4/ndsolve.html">Numerical solutions </a></li> <li><a class="subli" href="../ch4/picard.html">System conversion </a></li> <li><a class="subli" href="../ch4/power.html">Power series method </a></li> <li><a class="subli" href="../ch4/mdm.html">Modified decomposition method </a></li> <li><a class="subli" href="../ch4/euler.html">Euler's methods</a></li> <li><a class="subli" href="../ch4/RK4.html">Runge--Kutta methods </a></li> <li><a class="subli" href="../ch4/fdm.html">Finite difference methods</a></li> <li><a class="subli" href="../ch4/adm.html">Adomian decomposition method </a></li> <li><a class="subli" href="../ch4/vim.html">Variational iteration method </a></li> <li><a class="subli" href="../ch4/fem.html">Finite element method </a></li> <li><a href="../ch4/second.html">Second order ODEs</a></li> <li><a class="subli" href="../ch4/application4.html">Applications</a></li> </ol> </li> <li><a href="../ch5/part5.html">Part V:</a> <a href="javascript:clickit(5)">Fourier Series</a> <ol start="1" id="fifth" type="i"> <li><a class="subli" href="../ch5/sturm.html">Sturm--Liouville problems </a></li> <li><a class="subli" href="../ch5/singular.html">Singular Sturm--Liouville problems </a></li> <li><a class="subli" href="../ch5/transform.html">Fourier transform </a></li> <li><a class="subli" href="../ch5/fourier.html">Fourier series </a></li> <li><a class="subli" href="../ch5/periodic.html">Periodic extension</a></li> <li><a class="subli" href="../ch5/complex.html">Complex Fourier series</a></li> <li><a class="subli" href="../ch5/even.html">Even and odd functions</a></li> <li><a class="subli" href="../ch5/fexample.html">Examples</a></li> <li><a class="subli" href="../ch5/modes.html">Modes of convergence </a></li> <li><a class="subli" href="../ch5/gibbs.html">Gibbs phenomenon</a></li> <li><a class="subli" href="../ch5/converge.html">Convergence of Fourier series</a></li> <li><a class="subli" href="../ch5/cesaro.html">Cesàro summation</a></li> <li><a class="subli" href="../ch5/step.html">Square wave functions</a></li> <li><a href="../ch5/orthogonal.html">Orthogonal expansions </a></li> <li><a class="subli" href="../ch5/chebyshev.html">Chebyshev expansions</a></li> <li><a class="subli" href="../ch5/legendre.html">Legendre expansion</a></li> <li><a class="subli" href="../ch5/bessel.html">Bessel expansion</a></li> <li><a class="subli" href="../ch5/hermite.html">Hermite expansion</a></li> <li><a class="subli" href="../ch5/laguerre.html">Laguerre expansion</a></li> <li><a class="subli" href="../ch5/music.html">Motivated examples </a></li> </ol> </li> <li><a href="../ch6/part6.html">Part VI:</a> <a href="javascript:clickit(6)">Partial Differential Equations</a> <ol start="1" id="sixth" type="i"> <li><a class="subli" href="../ch6/first.html">First order PDEs</a></li> <li><a class="subli" href="../ch6/separate.html">Separation of variables </a></li> <li><a class="subli" href="../ch6/green.html">Green's functions</a></li> <li><a class="subli" href="../ch6/blasius.html">Blasius equation</a></li> <li><a class="subli" href="../ch6/fluid.html">Fluid problems</a></li> <li><a class="subli" href="../ch6/application6.html">Applications</a></li> <li><a class="subli" href="../ch6/miscellany.html">Miscellany</a></li> </ol> </li> <!---- <li> <a class="navSectionTitle" href="../ch6/parabolic.html">Parabolic equations</a> <a class="showDropdown" href="javascript:clickit(10)"> <b class="caret"></b></a> <ol start="1" id="tenth" type="a"> ----> <a href="../ch6/parabolic.html" style="color: #0A0A3C"><b>Part VI P: </b> </a> <a href="javascript:clickit(10)">Parabolic Equations</a> <ol start="1" id="tenth" type="i"> <li><a class="subli" href="../ch6/heat.html">Heat conduction equations</a></li> <li><a class="subli" href="../ch6/bheat.html">Boundary Value Problems for heat equation</a></li> <li><a class="subli" href="../ch6/oheat.html">Other heat transfer problems</a></li> <li><a class="subli" href="../ch6/hfourier.html">Fourier transform</a></li> <li><a class="subli" href="../ch6/fokas.html">Fokas method</a></li> <li><a class="subli" href="../ch6/resolvent.html">Resolvent method</a></li> <li><a class="subli" href="../ch6/fp.html">Fokker--Planck equation</a></li> <li><a class="subli" href="../ch6/nheat.html">Numerical solutions of heat equation </a></li> <li><a class="subli" href="../ch6/bs.html">Black Scholes model </a></li> <li><a class="subli" href="../ch6/pmonte.html">Monte Carlo for Parabolic </a></li> </ol> </li> <!---- <li> <a class="navSectionTitle" href="../ch6/hyperbolic.html">Hyperbolic equations</a> <a class="showDropdown" href="javascript:clickit(11)"> <b class="caret"></b></a> <ol start="1" id="eleventh" type="a"> ---> <br /> <a href="../ch6/hyperbolic.html" style="color: #0A0A3C"><b>Part VI H: </b> </a> <a href="javascript:clickit(11)">Hyperbolic equations</a> <ol start="1" id="eleventh" type="i"> <li><a class="subli" href="../ch6/wave.html">Wave equations</a></li> <li><a class="subli" href="../ch6/IBVP.html">IBVPs</a></li> <li><a class="subli" href="../ch6/2wave.html">2D wave equations</a></li> <li><a class="subli" href="../ch6/fwave.html">Forced wave equations</a></li> <li><a class="subli" href="../ch6/transverse.html">Transverse vibrations of beams</a></li> <li><a class="subli" href="../ch6/nwave.html">Numerical solutions of wave equation </a></li> <li><a class="subli" href="../ch6/kg.html">Klein–Gordon equation</a></li> <li><a class="subli" href="../ch6/3wave.html">3D wave equations</a></li> <li><a class="subli" href="../ch6/cagniard.html">Cagniard method</a></li> </ol> </li> <!---- <a class="navSectionTitle" href="../ch6/elliptic.html">Elliptic equations</a> <a class="showDropdown" href="javascript:clickit(12)"> <b class="caret"></b></a> <ol start="1" id="twelfth" type="a"> ---> <br /> <a href="../ch6/elliptic.html" style="color: #0A0A3C"><b>Part VI E: </b> </a> <a href="javascript:clickit(12)">Elliptic equations</a> <ol start="1" id="twelfth" type="i"> <li><a class="subli" href="../ch6/laplace.html">Laplace equation</a></li> <li><a class="subli" href="../ch6/dirichlet.html">Dirichlet problem</a></li> <li><a class="subli" href="../ch6/laplacen.html">Neumann problems for Laplace equation</a></li> <li><a class="subli" href="../ch6/laplacem.html">Mixed problems for Laplace equation</a></li> <li><a class="subli" href="../ch6/stripe.html">Laplace equation in infinite stripe</a></li> <li><a class="subli" href="../ch6/semistripe.html">Laplace equation in infinite semi-stripe</a></li> <li><a class="subli" href="../ch6/nlaplace.html">Numerical solutions of Laplace equation </a></li> <li><a class="subli" href="../ch6/polar.html">Laplace equation in polar coordinates</a></li> <li><a class="subli" href="../ch6/corner.html">Laplace equation in a corner</a></li> <li><a class="subli" href="../ch6/sphere.html">Laplace equation in spherical coordinates</a></li> <li><a class="subli" href="../ch6/poisson.html">Poisson's equation</a></li> <li><a class="subli" href="../ch6/helmholtz.html">Helmholtz equation</a></li> <li><a class="subli" href="../ch6/liouville.html">Liouville's equation</a></li> <li><a class="subli" href="../ch6/emonte.html">Monte Carlo for Elliptic </a></li> </ol> </li> <li><a href="../ch7/part7.html">Part VII:</a> <a href="javascript:clickit(7)">Special Functions</a> <ol start="1" id="seventh" type="i"> <li><a class="subli" href="../ch7/orthog.html">Orthogonal polynomials </a></li> <li><a class="subli" href="../ch7/gamma.html">Gamma function</a></li> <li><a class="subli" href="../ch7/bessel.html">Bessel functions</a></li> <a class="showDropdown" href="javascript:clickit(16)"> <b class="caret"></b></a> <ol start="1" id="sixteen" type="a"> <li><a class="subli" href="../ch7/bessel1.html">Generating functions </a></li> <li><a class="subli" href="../ch7/bessel2.html">Modified Bessel functions</a></li> <li><a class="subli" href="../ch7/bessel3.html">Hunkel functions </a></li> <li><a class="subli" href="../ch7/bessel4.html">Kelvin functions </a></li> <li><a class="subli" href="../ch7/bessel5.html">Recurrences </a></li> <li><a class="subli" href="../ch7/besselo.html">Orthogonality of Bessel functions</a></li> <li><a class="subli" href="../ch7/airy.html">Airy functions</a></li> <li><a class="subli" href="../ch7/bessel7.html">Applications </a></li> </ol> </li> <li><a class="subli" href="../ch7/chebyshev.html">Chebyshev functions</a></li> <a class="showDropdown" href="javascript:clickit(17)"> <b class="caret"></b></a> <ol start="1" id="seventeen" type="a"> <li><a class="subli" href="../ch7/cheby1.html">Generating functions </a></li> <li><a class="subli" href="../ch7/cheby2.html">Orthogonality </a></li> <li><a class="subli" href="../ch7/cheby3.html">Recurrences </a></li> <li><a class="subli" href="../ch7/cheby4.html">Zeroes</a></li> <li><a class="subli" href="../ch7/cheby7.html">Applications </a></li> </ol> </li> <li><a class="subli" href="../ch7/legendre.html">Legendre functions </a></li> <a class="showDropdown" href="javascript:clickit(18)"> <b class="caret"></b></a> <ol start="1" id="eighteen" type="a"> <li><a class="subli" href="../ch7/lege1.html">Generating functions </a></li> <li><a class="subli" href="../ch7/lege2.html">Orthogonality </a></li> <li><a class="subli" href="../ch7/lege3.html">Recurrences </a></li> <li><a class="subli" href="../ch7/lege4.html">Zeroes</a></li> <li><a class="subli" href="../ch7/lege7.html">Applications </a></li> </ol> </li> <li><a class="subli" href="../ch7/lambert.html">Lambert function</a></li> <li><a class="subli" href="../ch7/mathieu.html">Mathieu function </a></li> <li><a class="subli" href="../ch7/elliptic.html">Elliptic functions</a></li> <li><a class="subli" href="../ch7/hyper.html">Hypergeometric functions</a></li> <li><a class="subli" href="../ch7/kummer.html">Kummer's equation </a></li> <li><a class="subli" href="../ch7/application7.html">Miscellany</a></li> </ol> </li> </ol> </li> </ul> </div> </div> <p>Return to <a href="../../../am33/computing33.html">computing</a> page for the first course APMA0330<br /> Return to <a href="../../computing.html">computing</a> page for the second course APMA0340<br /> Return to <a href="../../../am33/Mathematica/index.html">Mathematica</a> tutorial for the first course APMA0330<br /> Return to <a href="../index.html">Mathematica</a> tutorial for the second course APMA0340<br /> Return to the <a href="../../../am33/index.html">main page</a> for the first course APMA0330<br /> Return to the <a href="../../index.html">main page</a> for the second course APMA0340<br /> Return to <a href="part5.html">Part V</a> of the course APMA0340 <br /> Introduction to <a href="http://www.cfm.brown.edu/people/dobrush/cs52/index.html">Linear Algebra</a> with <em>Mathematica</em> <h2 class="link"><a href="../../glossary/">Glossary</a></h2> <div style="clear: both;"></div> <div class="chapter" style="padding-top: 8px;"> <h1><a name="joseph" id="joseph"> Biography of Joseph Fourier </a></h1> <hr color="#cccccc" noshade="noshade" size="1" width="100%" /> <p> <a href="https://en.wikipedia.org/wiki/Joseph_Fourier">Jean-Baptiste Joseph Fourier</a> was born on the 21st of March 1768 in the ancient town of <a href="https://en.wikipedia.org/wiki/Auxerre">Auxerre</a>, Burgundy, Kingdom of France (now in Yonne, France). This small town has many cultural, religious, political and educational traditions. It was one of the foremost centers of teaching and learning in France in the ninth and tenth centuries. It is located on the bank of the river <a href="https://en.wikipedia.org/wiki/Yonne_%28river%29">Yonne</a> with many magnificent tall building including the famous ancient Clock Tower, the even more famous <a href="https://en.wikipedia.org/wiki/Abbey_of_Saint-Germain_d'Auxerre">Abbey St. Germain</a> founded by St. Germain himself in the fifth century A.D., and the Gothic <a href="https://en.wikipedia.org/wiki/Auxerre_Cathedral">Cathedral St. Etienne</a>. It became more famous as the birthplace of Joseph Fourier, one of the most brilliant mathematical scientists France has ever produced. <p> Joseph’s father was a master tailor of Auxerre. After the death of his first wife, with whom he had three children, he remarried and Joseph was the ninth of the twelve children of this second marriage. Joseph's mother was a housewife. When he was nine- or ten-years old, both his parents died. He first entered the Pallais elementary school to receive lessons in Latin and French. He was a gifted student in both languages, but showed extraordinary talents for mathematics and mechanics. At the age of 13, he had completed an extensive study of six volumes of Etienne Bézzout’s (1730--1783) Cours de mathématique. In 1783, he received the first prize for his study of C. Bossut’s (1730--1814) Mécanique. At a very young age, Joseph presented a research paper on algebraic equations at the Académie Royale des Sciences. In 1787, Fourier decided to train for the priesthood and entered the Benedictine abbey of <a href="https://en.wikipedia.org/wiki/Saint-Beno%C3%AEt-sur-Loire">Saint-Benoit-sur-Loire</a>. Later, however, he gave up his idea of studying theology. He was then admitted to the progressive École Royale Militaire to study science and mathematics. In 1790, he became a mathematics teacher at the Benedictine College, <a href="http://www.jacques-amyot.org/college_et_ecole_militaire-auxerre.pdf">École Royale Militaire of Auxerre</a>, where he had been a student. <p> The <a href="https://en.wikipedia.org/wiki/French_Revolution">French Revolution</a> was a watershed event in modern European history that began in 1789 and ended in the late 1790s with the ascent of Napoleon Bonaparte. In 1789, Joseph returned to Auxerre and received a teaching position in his old school. Fourier took an active role in the French Revolution, as he was keenly attracted to its egalitarian ideals. He joined his local Revolutionary Committee, but soon regretted it. During his time in the Committee, he made the mistake of confronting a rival sect in defense of his own faction while on a trip to <a href="https://en.wikipedia.org/wiki/Orl%C3%A9ans">Orléans</a>, consequences of which would later threaten his life. Certainly, Fourier was unhappy about the <a href="https://www.britannica.com/event/Reign-of-Terror">Terror</a> that resulted from the French Revolution and he attempted to resign from the committee. However this proved impossible and Fourier was now firmly entangled with the Revolution and unable to withdraw. After Fourier's return to Auxerre, he continued to work on the Revolutionary Committee and continued to teach at the College. In July 1794 he was arrested by the order of the Committee for Public Safety for his possible revolutionary and terrorist activities in the Revolution in Auxerre during 1793--1794, the charges relating to the Orléans incident, and he was imprisoned. Fourier feared that he would go to the guillotine but, after <a href="https://en.wikipedia.org/wiki/Maximilien_Robespierre">Maximilien Robespierre</a> himself went to the guillotine, political changes resulted in Fourier being freed. <p> Later in 1794, now free, Joseph was nominated to study at the <a href="https://en.wikipedia.org/wiki/%C3%89cole_normale_sup%C3%A9rieure_%28Paris%29">École Normale</a> in Paris, which opened in 1795 as a model school for training teachers in France. He joined this teacher-training school and got a rare opportunity to meet some of the foremost French mathematicians of that time including <a href="https://en.wikipedia.org/wiki/Joseph-Louis_Lagrange">J.L. Lagrange</a> (1736--1813), <a href="https://en.wikipedia.org/wiki/Pierre-Simon_Laplace">P.S. Laplace</a> (1749--1827) and <a href="https://en.wikipedia.org/wiki/Gaspard_Monge">Gaspard Monge</a> (1746--1818). In September 1795, Fourier was appointed as an assistant lecturer at the École Polytechnique to support the teaching of Lagrange and Monge. In a short period of two years, he succeeded Lagrange to occupy the Chair of Analysis and Mechanics at the Polytechnique in 1797. In spite of his cordial relationship with Lagrange, Laplace and Monge, and his exceptional teaching record, Fourier’s research record at the time was not that outstanding by any standard. <p> In 1798, Fourier was selected to join <a href="https://en.wikipedia.org/wiki/Napoleon">Napoleon</a> Bonaparte’s (1769--1821) <a href="https://en.wikipedia.org/wiki/French_campaign_in_Egypt_and_Syria">Egyptian expedition</a> as scientific adviser to liberate unhappy Egyptian people from deep social and political problems and to provide them with all the benefits of European civilization (NB Fourier himself was opposed to the abusive treatment inflicted by the colonial rulers, and during his stay in Egypt, worked almost solely within the field or pure mathematics). He remained in this position until the downfall of Napoleon in 1815. During his occupation of Egypt, Napoleon established the <a href="https://en.wikipedia.org/wiki/Institut_d'%C3%89gypte">Institute d’Egypt</a> in <a href="https://en.wikipedia.org/wiki/Cairo">Cairo</a> with Monge as its President and Fourier as Permanent Secretary. Fourier continued his position during the entire occupation of <a href="https://en.wikipedia.org/wiki/Egypt">Egypt</a> and was fully responsible for collating all literary and scientific discoveries made during that time. In addition, he was very active in administering the archaeological research of the Cairo Institute. This position led Fourier to publish the Description de l’Egypte (Description of Egypt) with an elegant preface, partly edited by Napoleon himself. <p> In 1799, <a href="https://en.wikipedia.org/wiki/Napoleon">Napoleon</a> returned to Paris to hold the absolute power of France and subsequently crowned himself Emperor of the French. In 1801, Fourier also came back to Paris to resume his position as Professor of Analysis at the <a href="https://en.wikipedia.org/wiki/%C3%89cole_Polytechnique">École Polytechnique</a>. After his return from Egypt in 1802, Fourier was appointed Prefect of Isére in <a href="https://en.wikipedia.org/wiki/Grenoble">Grenoble</a> with many major and varied administrative duties and responsibilities. As Prefect from 1802 to 1815, Fourier was totally successful in planning, management and the draining operation of the largest swamps of Bourgoin, as the swamps had been fully responsible for annual epidemics of dangerous fever causing deaths of a large number of the surrounding inhabitants at a young age. After completion of the draining of the swamps in 1812, the annual epidemics of fever stopped leading to a marked improvement in health problems of the inhabitants. This was one of the greatest public and social service accomplishments of Fourier. His other major accomplishment as Prefect was the planning of a very long, spectacular highway from Grenoble to <a href="https://en.wikipedia.org/wiki/Turin">Turin</a>, of which only the French Section was built. Despite his unfriendly relationship with Napoleon, Napoleon awarded Fourier the title of <a href="https://en.wikipedia.org/wiki/Baron">Baron</a> in recognition of his major accomplishments as Prefect. <p> During 1787--1810, a series of personal letters between Fourier and C.L. Bonard (1765--1819), a mathematics teacher at Auxerre provided the unique opportunity to establish a real and permanent friendship between them. These letters were full of their many personal and family experiences and reminiscences of the old days, especially during the revolution in Auxerre. As a family friend and true admirer of Fourier, Bonard was always very proud of Fourier’s mathematical and other achievements. Unfortunately, Bonard died in 1819 just three years before the publication of Fourier’s masterpiece work on the Theory of Heat. <p> This masterpiece, which included his work on the Fourier series, was actually finished in 1807 as part of his publication 'On the Propagation of Heat in Solid Bodies'. It was not published until 1822, however, for his first reading of it (to the Paris Institute) in 1807 was met with very mixed and critical reception. Both Lagrange and Laplace objected to the notion of what we now call Fourier series: the expansions of functions as trigonometrical series. Along with another scientist, <a href="https://en.wikipedia.org/wiki/Jean-Baptiste_Biot">Jean-Baptiste Biot</a> (1774--1862), Lagrange and Laplace also objected to Fourier’s derivation of the equations of transfer of heat. (Biot had written an earlier paper on the topic in 1804, although that paper proved incorrect.) Nonetheless, when the Paris Institute held a competition on the topic of how heat propagates in solid bodies in 1811, Fourier submitted his memoir for consideration. He won the prize, in part because only one other entry was received. </p> <p> Because of the controversy, Fourier’s memoir was not published until 1822, after his election to the Académie des Sciences in 1817, and the same year he became the Académie’s secretary. His work did contain flaws, but it also provided the basis for later work on trigonometric series and the theory of functions of a real variable, most notably the Fourier transform, an operation that turns one function of a real variable into another. It is widely used in digital signal processing, as well as in the physical study of wave motion and optics. </p> <p> In 1827, like Jean D’Alembert (1717–1783) and Laplace before him, Fourier was elected to the literary Académie Francaise. In the same year, after the death of Laplace, Fourier was also elected President of the Conseil de perfectionnement of the École Polytechnique. Several foreign scientific organizations including the Royal Society of London and the Royal Swedish Academy of Sciences elected Fourier as honorary member for his outstanding contributions to mathematical sciences. Fourier continued to publish papers on mathematics until his death in 1830, when he tripped and fell down the stairs at home. </p> <p> Throughout his life, Fourier’s general health was not very good. Towards the end of his life, he began to display some unusual symptoms which were thought to have possibly been contracted during his stay in Egypt from the extremely hot climate. In early May of 1830, Fourier suffered from a heart attack and his condition gradually deteriorated until he died on May 16 at the age of 62. He was buried at the cemetery of <a href="https://en.wikipedia.org/wiki/P%C3%A8re_Lachaise_Cemetery">Père Lachaise</a> in Paris. His gravestone was decorated with an Egyptian motif to recognize his outstanding service as Permanent Secretary of the Cairo Institute and his collation of the landmark <a href="https://en.wikipedia.org/wiki/Description_de_l'%C3%89gypte">Description de l'Égypte</a>. <p> </p> <div class="chapter"> <p class="custom">Fourier and heat conduction </p> <hr color="#cccccc" noshade="noshade" size="1" width="100%" /> </div> <p> During his stay in Grenoble, Fourier found the time and intellectual energy to discover the theory of heat conduction based on his new mathematical ideas, methods, experiments and observations. He completed his classic publication on the conduction of heat entitled ‘On the propagation of heat in solid bodies’, and submitted it to the Academy of Sciences of Paris in 1807 for a research prize. His article was judged by a great scientific committee consisting of <a href="https://en.wikipedia.org/wiki/Joseph-Louis_Lagrange">Lagrange</a>, <a href="https://en.wikipedia.org/wiki/Pierre-Simon_Laplace">Laplace</a>, <a href="https://en.wikipedia.org/wiki/Gaspard_Monge">Monge</a>, <a href="https://en.wikipedia.org/wiki/Adrien-Marie_Legendre">A.M. Legendre</a> (1752--1833) and <a href="https://en.wikipedia.org/wiki/Sylvestre_Fran%C3%A7ois_Lacroix">S.L. Lacroix</a> (1765--1843). The committee was generally impressed by the importance and novelty of Fourier’s work, but it was split on the overall quality and rigor of the mathematical treatment involving the use of trigonometric (Fourier) series expansions of functions without adequate justification. The controversy created by his work on trigonometric expansions of functions was due to his demonstration of the paradoxical property of equality in a finite interval between algebraic results of a totally different form. Another major criticism of his publication was about his derivation of the heat equation in a continuous solid was based on inadequate physical principles. Unfortunately, his memoir was rejected by the committee for the prize. Although some of his contemporary mathematicians vociferously objected to his revolutionary work, the Academy was generally impressed by the subject matter of Fourier’s memoir and encouraged him to resubmit it for a Grand Prize to be awarded in 1812. At the same time, <a href="https://en.wikipedia.org/wiki/Sim%C3%A9on_Denis_Poisson">S.D. Poisson</a> (1781--1840) and <a href="https://en.wikipedia.org/wiki/Jean-Baptiste_Biot">J.B. Biot</a> (177--1862), his rivals, expressed serious concerns with major criticisms of Fourier’s work on his mathematical and physical treatment of the theory of heat in solids. Surprisingly, Biot praised Poisson’s less original work on the same subject. Poisson’s harsh criticism is quoted below in full: <blockquote> As the partial differential equation to which it corresponds is linear and has constant coefficients, one can also satisfy it by an integral composed of an infinity of exponentials of sines and cosines containing an infinite number of arbitrary constants: this integral is contained in the preceding one; but it would be difficult to decide a priori if it has the same degree of generality and if it can replace it identically, something which necessarily throws doubt and obscurity on all solutions deduced from this second form of the integral. M. Fourier, who did not go beyond a solution of this kind, remarks himself that it is similar to that which Daniel Bernouilli gave to the problem of vibrating strings; but it is well known that Euler, d’Alembert, and Lagrange, who occupied themselves at the same time with the same problem, and who differed among themselves on various points, were at one nevertheless in regarding Bernouilli’s solutions as incomplete and less general than that containing arbitrary functions. This is not true of the formulas of M. Fourier: I am sure that all the results he obtains are correct; but against his analysis can be advanced the same objections as those advanced against that of Bernouilli and repeated in other similar cases. <p> In general it seems to me that whenever an unknown quantity depends on a partial differential equation, and when its values should reduce in fact to a sum of particular integrals, the only way of disposing of all doubts and retaining for the mathematical certainty result is not to suppose in advance such a form for the unknown quantity, but to deduce it, on the contrary, from the general integral by a succession of direct and rigorous transformations. This is what I have attempted to do in this memoir . . . <p> . . . I leave it to mathematicians to judge if I have attained the end that I have set myself. </blockquote> <p> Fourier had real talent, and extraordinary courage to defend his work in an absolutely unprecedented manner. His devastating response to Poisson and Biot’s unfounded criticisms of his solution was contained in his unpublished <i>Historical Pré́cis</i>, and in his letters to Laplace around 1808--1809 with a note on certain analytic expressions in connection with the equations of the theory of heat. Thereafter, all criticisms seemed to have subsided once, and for all, barely having tarnished his permanent reputation. <p> Fourier continued his research and resubmitted his revised article in 1811 with some addition of new material on the cooling of infinite solids, radiant and terrestrial heat, comparison of his theory with experimental observations and equations of the movement of heat in fluids. In spite of some disagreements, Fourier was awarded the Grand Prize for his work in 1812, but his work was not recommended for publication in the <i>Mé́moires</i> of the Academy. In 1817, he was elected member of the Acadé́́mie des Sciences by an overwhelming majority. He continued to revise his work for publication. In 1822, <a href="https://en.wikipedia.org/wiki/Jean_Baptiste_Joseph_Delambre">J.B. Delambre</a> (1749--1822), who was Permanent Secretary of the mathematics and physics division of the Acadé́mie des Sciences died, and Fourier was elected as his replacement. Fourier remained in this position until his death in 1830. Shortly after he became the Secretary, the Acadé́mie published his masterpiece work entitled <i>La Thé́́orie Analytique de la Chaleur</i> (The Analytical Theory of Heat) based on <a href="https://en.wikipedia.org/wiki/Isaac_Newton">Isaac Newton’s</a> (1642--1727) <a href="https://en.wikipedia.org/wiki/Newton's_law_of_cooling">law of cooling</a>, that the flow of heat between two adjacent molecules is proportional to extremely small temperature difference. In his preface of this great book, Fourier described his own assessment of the work as follows: <blockquote> First causes are not known to us, but they are subjected to simple and constant laws that can be studied by observation and whose study is the goal of Natural Philosophy . . . Heat, like gravity penetrates every substance of the universe; its rays occupy all regions of space. The object of our work is to set forth the mathematical laws that this element obeys . . . But whatever the extent of the mechanical theories, they do not apply at all to the effects of heat. They constitute a special order of phenomena that cannot be explained by principles of movement and of equilibrium . . . The differential equations of the propagation of heat express the most general conditions, and reduce physical questions to problems in pure analysis that is the proper object of the theory. </blockquote> <p> His analytical theory of heat was highly mathematically based on a general equation of propagation of heat subject to appropriate initial and boundary conditions of different kinds. Fourier’s fundamental heat equation for the temperature distribution, <i>u</i> = <i>u</i>(<b>x</b>, t) is given by <div class="math"> \[ u_t = \kappa\,\nabla^2 u , \] </div> where <b>x</b> = (x,y,z), <span class="math"> \( \kappa = \frac{k}{c\,\rho} \) </span> is the diffusivity constant, <i>k</i> is the heat conductivity, <i>c</i> is the specific heat, ρ is the density, and <span class="math"> \( \nabla^2 = \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2} \) </span> is the three-dimensional Laplacian. Thus, the heat equation is of the first order in time derivative and of the second order in space derivatives. This was a striking contrast with almost all dynamic equations, including the wave equation in strings or air that are of second order in both time and space derivatives. </p> <p> One of the French scientific philosophers and close colleagues of Fourier, <a href="https://en.wikipedia.org/wiki/Auguste_Comte">Auguste Compte</a> (1798--1857) attended a course of lecture of Fourier on the role of the positivist philosophy of science. In his published book entitled <i>Cours de Philosophie</i> positive dedicated to Joseph Fourier, Compte expressed his eloquent and complimentary views on Fourier’s celebrated theory of heat as follows: <blockquote> In fact, in this work, of which the philosophical characteristic is so eminently positive the most important and precise laws of thermal phenomena are discovered without the author having once inquired about the intimate nature of heat. </blockquote> Subsequently, <a href="https://en.wikipedia.org/wiki/Henri_Poincar%C3%A9">Henri Poincaré́</a> (1854--1912), an extremely versatile and creative mathematician and mathematical physicist, made a delightful comment on Fourier’s fundamental work as follows: <blockquote> Fourier’s theory of heat is one of the first examples of the application of analysis to physics. Starting from simple hypotheses, which are nothing but generalized facts, Fourier deduced from them a series of consequences which together make up a complete and coherent theory. The results which he obtained are certainly interesting in themselves, but what is still more interesting is the method which he used to arrive at them and which will always be a model for all those who wish to cultivate any branch of mathematical physics. </blockquote> <p> Even with lack of rigor in Fourier’s classic work on the analytical theory of heat, a eminent British physicist, <a href="https://en.wikipedia.org/wiki/William_Thomson,_1st_Baron_Kelvin">William Thomson</a>, 1st Baron Kelvin (1824--1907) regarded the work as Fourier’s permanent legacy and called it ‘a great mathematical poem.’ <p> </p> <div class="chapter"> <p class="custom">The Fourier integral theorem </p> <hr color="#cccccc" noshade="noshade" size="1" width="100%" /> </div> <p> This section includes some of Fourier’s major contributions to other areas of mathematics and physics and their foremost impact on mathematical physics, mathematical economics, probability, statistics, pure and applied mathematics. During 1817--1818, he made notable contributions to the theory of Fourier transforms and their applications to partial differential equations including the heat and the wave equations. His 1822 monumental treatise on the <i>Analytical Theory of Heat</i> provided the modern mathematical theory of conduction of heat, Fourier series and Fourier integrals with many new examples of applications. In his treatise, he formulated a remarkable theorem which is now universally known as the <a href="https://en.wikipedia.org/wiki/Fourier_inversion_theorem">Fourier Integral Theorem</a> that states that an arbitrary function <i>f</i>(<i>x</i>) has the Fourier double integral representation in the form <div class="math"> \[ f(x) = \frac{1}{2\pi} \int_{-\infty}^{\infty} {\text d}k\, e^{{\bf j}kx} \int_{-\infty}^{\infty} {\text d}\xi\, f(\xi )\, e^{-{\bf j}kx} . \] </div> This representation shows that <i>f</i>(<i>x</i>) has been decomposed into an infinite number of waves of different wavenumbers <i>k</i> (or wavelengths λ = 2π/<i>k</i>), and amplitude <span class="math"> \( \frac{1}{2\pi} \int_{-\infty}^{\infty} {\text d}\xi\, f(\xi )\, e^{-{\bf j}kx} , \) </span> whereas a given function can be represented by a Fourier series in terms of an infinite but discrete set of wave components. </p> <p> This theorem has been expressed in several slightly different forms to better adapt it for many particular applications. It has been recognized, almost from the start, however, that the form which best combines mathematical simplicity and complete generality makes use of the exponential oscillating function exp(<b>j</b>kx). Indeed, the Fourier integral formula is regarded as one of the most fundamental results of modern mathematical analysis, and it has widespread physical and engineering applications. The generality, applicability and importance of the theorem is eloquently expressed by Lord Kelvin and Scottish mathematician <a href="http://www-history.mcs.st-and.ac.uk/Biographies/Tait.html">Peter Guthrie Tait</a> (1831--1901) as follows: ‘. . . Fourier’s Theorem, which is not only one of the most beautiful results of modern analysis, but may be said to furnish an indispensable instrument in the treatment of nearly every recondite question in modern physics. To mention only sonorous vibrations, the propagation of electric signals along a telegraph wire and the conduction of heat by the earth’s crust, as subjects in their generality intractable without it, is to give but a feeble idea of its importance’. This integral formula is usually used to define the classical Fourier transform of a function and the inverse Fourier transform. No doubt, the scientific achievements of Joseph Fourier have not only provided the fundamental basis for the study of heat equation, Fourier series and Fourier integrals, but also for the modern developments of the theory and applications of the partial differential equations. <p> Although Fourier is most celebrated for his work on the conduction of heat, Fourier series and Fourier transforms, the new mathematics he created has proved to be very useful in a wide variety of ways. His superb mastery of analytical treatment and physical intuition played a fundamental role in his success. In mathematical physics, Fourier’s contributions have proved to be much more useful than those of his renowned contemporaries including Lagrange, Laplace, Poisson, Biot and Augustin Cauchy (1789–1857). Both Cauchy and Poisson successfully applied Fourier’s universal mathematical analysis to the subject of linear water waves and to other areas. It has also been a remarkable fact that Fourier analysis eventually resolved the long-standing controversy over the mathematical treatment of the famous vibrating string problem which was vigorously investigated by several eminent mathematical scientists including Daniel Bernoulli (1700–1782), Leonhard Euler (1707–1783), D’Alembert and Lagrange. Fourier’s work was also remarkably influential for the comprehensive development of James Clerk Maxwell’s (1831–1879) electromagnetic theory and kinetic theory of gases. In 1864, Maxwell first formulated the basic equations for all electromagnetic phenomena, known as the Maxwell equations. Maxwell expressed his full admiration for Fourier’s great work on the theory of heat. <p> Among his other research contributions to physical sciences, Fourier first suggested that the earth’s atmosphere helps trap the heat of sunlight. This special behavior of the atmosphere is known as the greenhouse effect because it resembles the action of the green roof of a greenhouse. In his series of articles published during 1824 and 1827, he investigated various possible sources of the additional observed heat, and suggested that interstellar radiation may be responsible for a large amount of warmth. Apart from this work, Fourier published many papers on terrestrial heat and radiant heat which led to his mathematical derivation of the sine law for emission of radiation at the surface of heated bodies based on <a href="https://en.wikipedia.org/wiki/John_Leslie_%28physicist%29">John Leslie’s</a> (1766--1832) experimentally determined the sine law. He also studied the application of the principle of virtual velocities to fluid flows, oscillations of system of bodies about the state of equilibrium, and the subject of wave motions in elastic lamina. He also took up a few aspects of Laplace’s work on probability and statistics dealing with the significance of the normal distribution and the central limit theorem. In particular, his interest in the density function <span class="math"> \( f(x) = \frac{1}{\sqrt{2\pi}} \,\exp \left( - x^2 /2 \right) \) </spam> of the 2 standard normal (or Gaussian) distribution occurs in a wide variety of mathematical contexts, from probability theory to Fourier analysis to quantum mechanics. This density function has some unique properties that are not shared by any other special functions in mathematics and mathematical physics. <p> </p> <div class="chapter"> <p class="custom">Fourier series </p> <hr color="#cccccc" noshade="noshade" size="1" width="100%" /> </div> <p> Historically, it was Fourier who first not only represented an arbitrary function in trigonometric series of cosine and sine functions, but also raised the fundamental question of convergence of the series. Once he made the following bold statements: ‘Thus, there is no function . . . which cannot be expressed by a trigonometric series . . . [or] definite integral’. ‘Regarding the researches of <a href="https://en.wikipedia.org/wiki/Jean_le_Rond_d'Alembert">d’Alembert</a> and <a href="https://en.wikipedia.org/wiki/Leonhard_Euler">Euler</a> could one not add that if they knew this expansion they made but a very imperfect use of it. They were both persuaded that an arbitrary and discontinuous function could never be resolved in series of this kind, and it does not seem that anyone had developed a constant in cosines of multiple arcs, the first problem which I had to solve in the theory of heat’. The first major result about the convergence of a series of functions was <a href="https://en.wikipedia.org/wiki/Peter_Gustav_Lejeune_Dirichlet">Peter Gustav Lejeune Dirichlet’s</a> (1805--1859) theorem of convergence of Fourier series of piecewise continuous functions. In 1829, he was one of the first mathematicians who proved the first theorem giving sufficient (and very general) conditions for the Fourier series of function <i>f</i>(<i>x</i>) to converge pointwise to <i>f</i>(<i>x</i>) after 18 years of Fourier’s discovery. The proof given by Dirichlet was a kind of refinement of that sketched by Fourier in the final section of his treatise on the theory of heat. Dirichlet’s method of proof was to calculate the limit of the <i>n</i>th partial sum of a Fourier series, <i>s</i><sub>n</sub>(<i>x</i>) as <i>n</i> → ∞. He proved that for any given value of <i>x</i>, the sum of the Fourier series is <i>f</i>(<i>x</i>) if <i>f</i>(<i>x</i>) is continuous at that point <i>x</i>, and is <span class="math"> \( \frac{1}{2} \left[ f(x+0) +f(x-0) \right] \) </span> if <i>f</i>(<i>x</i>) is discontinuous at that point. It was Dirichlet who closely followed Fourier’s work and then in 1829 established for the first time a rigorous mathematical theory of Fourier series. In 1835, he also established the Poisson summation formula for the Fourier analysis, and then modernized the concept of function that served as the fundamental basis for the nineteenth century investigations of mathematical analysis. In connection with the study of convergence of Fourier series, Dirichlet commended Fourier’s great work: ‘The sine and cosine series, by which one can represent an arbitrary function in a given interval, enjoy among other remarkable properties that of being convergent. This property did not escape the great geometer (Fourier) who began, through the introduction of representation of functions just mentioned, a new career for the applications of analysis; it was stated in the article that contains his first research on heat. But no one so far, to my knowledge, gave a general proof of it’. Inspired by Dirichlet’s work on Fourier series, <a href="https://en.wikipedia.org/wiki/Bernhard_Riemann">Bernard Riemann</a> (1826--1866) first formulated the concept of Riemann integral in his memoir on trigonometric series in 1854 that served the fundamental basis for his profound work on the convergence of Fourier series. <p> According to his fundamental theorem, for a bounded integrable function <i>f</i>(<i>x</i>), the convergence of its Fourier series at a point <i>x</i> in [-π, π] depends only on the behavior of <i>f</i>(<i>x</i>) in an arbitrarily small neighborhood of that point <i>x</i>. It was Riemann who first observed this behavior and referred it as the localization principle for the convergence of Fourier series. Both Dirichlet and Riemann recognized that more than continuity is required to ensure the convergence of Fourier series, namely the requirement of bounded variation of the function <i>f</i>(<i>x</i>) in a neighborhood of <i>x</i>. In 1854, Riemann praised Fourier’s work and said: ‘Nearly fifty years had passed without any progress on the question of analytic representation of an arbitrary function, when an assertion of Fourier threw new light on the subject. Thus a new era began for the development of this part of Mathematics and this was heralded in a stunning way by major developments in mathematical Physics’. Further investigations, continuing to the present, led to extensions of this result and revealed that the general question of convergence of a Fourier series is a very difficult and deep subject in mathematical sciences. However, the problem of finding necessary and sufficient conditions on <i>f</i>(<i>x</i>) so that its Fourier series converges to <i>f</i>(<i>x</i>) remained unsolved for a period of century. The nature of the convergence of Fourier series also received further attention after the introduction of the notion of uniform convergence by several leading mathematical scientists including George <a href="https://en.wikipedia.org/wiki/Sir_George_Stokes,_1st_Baronet">Gabriel Stokes</a> (1819--1903), <a href="https://en.wikipedia.org/wiki/Philipp_Ludwig_von_Seidel">Philipp L. Seidel</a> (1821--1896) and <a href="https://en.wikipedia.org/wiki/Karl_Weierstrass">Karl Weierstrass</a> (1815--1897). Subsequently, based on measure theory <a href="https://en.wikipedia.org/wiki/Henri_Lebesgue">Henri Lebesgue’s</a> (1875--1941) theory of integral and convergence in <i>L</i>²-norm were put Fourier analysis on the modern mathematical map. His discovery of these fundamental ideas and methods undoubtedly stimulated further development of Fourier series and integrals. Despite various attempts, it was an open question and unsolved problem for a period of a century whether a Fourier series of a continuous function converges at a point. In 1915, a Russian mathematician <a href="https://en.wikipedia.org/wiki/Nikolai_Luzin">Nikolai Luzin</a> (1883--1950) conjectured that if <i>f</i> ∈ <i>L</i>²[0, 2π], then its Fourier series converges almost everywhere. In 1966, <a href="https://en.wikipedia.org/wiki/Lennart_Carleson">Lennart Carleson</a> (born in 1928) provided an affirmative answer of Luzin’s conjecture with the deepest theorem which states that the Fourier series of any square integrable function <i>f</i>(<i>x</i>) in [0, 2π] (i.e. <i>f</i> ∈ <i>L</i>²[0, 2π]) converges almost everywhere to <i>f</i>(<i>x</i>). Shortly, after Carleson’s great work, Hunt proved that almost everywhere convergence of Fourier series of functions in <i>L</i><sup>p</sup>[0, 2π] for 1 < <i>p</i> < ∞. One year later, Sjolin gave further refinement of Hunt’s theorem to include spaces of functions that are smaller than <i>L</i>¹, but larger than <i>L</i><sup>p</sup> for every <i>p</i> > 1. Subsequently, Hunt and Taibleson proved the sharpest result on almost everywhere convergence of Fourier series on the ring of integers of a local field. However, another great mathematician of the twentieth century, <a href="https://en.wikipedia.org/wiki/Andrey_Kolmogorov">A.N. Kolmogorov’s</a> (1903--1987) 1926 classical counter example shows that the result fails when <i>p</i> = 1, that is, the Fourier series of <i>L</i>¹ functions diverges everywhere. This resolved an important question of convergence of Fourier series which became very useful in mathematical analysis and in many other areas of mathematical sciences. <p> </p> <div class="chapter"> <p class="custom">Fourier pure mathematics </p> <hr color="#cccccc" noshade="noshade" size="1" width="100%" /> </div> <p> Fourier had early strong interests in pure mathematics. In fact, Bezout’s book on <i>Théorie gé́nérale des equations algébbraiques</i> (1779) stimulated Fourier’s active research interest in classical algebra. At the age of 21 in 1789, he presented a paper on algebraic equations at the Académie Royale des Sciences. He gave a new proof of the famous rule of <a href="https://en.wikipedia.org/wiki/Ren%C3%A9_Descartes">René Descartes</a> (1596--1650) which describes how to determine the number of positive and negative roots of algebraic polynomial equations. His lifelong interest in the theory of equations was evident from his incomplete book manuscript on determinate equations which was edited by his close friend and a famous applied mathematician, <a href="https://en.wikipedia.org/wiki/Claude-Louis_Navier">Claude-Louis Navier</a> (1785--1836) and then published it posthumously with the title <i>Fourier’s Analyse des équations détterminées</i> in 1831. Fourier’s other ideas on algebra eventually led to his student <a href="https://en.wikipedia.org/wiki/Jacques_Charles_Fran%C3%A7ois_Sturm">Charles Sturm’s</a> (1803--1855) basic theorem. The solution of Fourier’s general equation of heat conduction in two or three dimensions required precise boundary conditions. Fourier’s first mathematical treatment and then by his students Charles Sturm and Joseph Liouville (1809–1882) opened up the whole new field of eigenvalues and eigenfunctions of ordinary and partial differential equations of such enormous importance for modern applied mathematics and mathematical physics. Apparently, Fourier himself realized the extraordinary nature and true significance of his discovery of the equation of heat conduction and its methods of solution by Fourier series and Fourier transforms. This work displayed his unique physical understanding and mathematical ingenuity. Other modern developments in pure mathematics which can be traced back to Fourier’s analytical theory of heat included the theory of orthogonal functions, Fourier’s expansion of an arbitrary function in terms of eigenfunctions, the ideas of different types of convergence, definite integral as the limit of a sum and theory of infinite determinants. Thus, Fourier’s most important, and, indeed, revolutionary contribution to pure mathematics was the representation of functions by trigonometric (Fourier) series, its numerous mathematical properties including the theory of convergence, term-by-term integration and differentiation, and then, the determination of the sum of numerical infinite series which is often impossible by other methods of summation of infinite series. There are examples which showed that the constructed Fourier series whose values on a given interval may not, in general, be the same as those of the function. <p> <a href="https://en.wikipedia.org/wiki/Albert_Einstein">Albert Einstein</a> (1879--1955) well-described the mutual relationship of pure mathematics, applied mathematics and physical laws as: ‘pure mathematics enables us to discover the concepts and laws connecting them, which gives us the key to the understanding of the phenomena of nature’. <p> Fourier went even further than this by insisting on the way in which situations observed in nature stimulate advanced study and research in pure mathematics: ‘The profound study of nature is the most fertile source of mathematical discoveries. Not only does this study, by offering a definite goal to research, have the advantage of excluding vague questions and futile calculations, but it is also a sure means of moulding analysis itself, and discerning those elements in it which it is still essential to know and which science ought to conserve. These fundamental elements are those which recur in all natural phenomena’. He was right to make such a profound statement: some of the most abstract pure mathematics involved in the theory of Fourier series and Fourier transforms have been used successfully to solve many applied physical problems including the vibration of elastic strings, and to develop the analytical theory of flow of heat in solid bodies of different geometric configurations. During the subsequent development of Fourier series, the English mathematician <a href="https://en.wikipedia.org/wiki/Henry_Wilbraham">Henry Wilbraham</a> (1825--1883) gave a remarkable surprise by his discovery in 1848 what is now known as the <a href="https://en.wikipedia.org/wiki/Gibbs_phenomenon">Gibbs phenomenon</a> (called by <a href="https://en.wikipedia.org/wiki/Maxime_B%C3%B4cher">Maxime Bôcher</a> after <a href="https://en.wikipedia.org/wiki/Josiah_Willard_Gibbs">Joseph Willard Gibbs</a> who discovered it in 1899) which states that near a jump discontinuity of a function, its Fourier series overshoots (or undershoots) it by approximately 9% of the jump. <p> Admittedly, the idea of convolution of two or more functions arose naturally first in the study of Fourier series and Fourier transforms. It also played a fundamental role in Fourier analysis, number theory, probability theory, mathematical statistics, harmonic analysis and served more generally in mathematical analysis of functions in other settings. <a href="Georg Cantor">Georg Cantor</a> (1845--1918) revolutionized mathematics in the second half of the nineteenth century by creating the foundation of the set theory which was motivated in part by questions of convergence and uniqueness of trigonometric and Fourier series expansions. In the 1920s, <a href="https://en.wikipedia.org/wiki/Frigyes_Riesz">Frigyes Riesz</a> (1880--1956) established the norm convergence for Fourier series, and almost simultaneously, <a href="https://en.wikipedia.org/wiki/Ernst_Sigismund_Fischer">Ernst Fisher</a> (1875--1959) unequivocally proved the mean-square convergence theorem for Fourier series. All these developments led to the theory of singular integral operators at the center of harmonic analysis. Fourier analysis, in general and Fourier transforms, in particular, have been very useful for the solution of linear partial differential equations with initial and boundary conditions of different kinds. <p> </p> <div class="chapter"> <p class="custom">Fourier’s heat equation and Brownian motion </p> <hr color="#cccccc" noshade="noshade" size="1" width="100%" /> </div> <p> It was shown by <a href="https://en.wikipedia.org/wiki/Albert_Einstein">Albert Einstein</a> that the Fourier heat equation is closely associated with the <a href="https://en.wikipedia.org/wiki/Diffusion_process">diffusion process</a> of <a href="https://en.wikipedia.org/wiki/Brownian_motion">Brownian motion</a>. In 1828, the Scottish botanist <a href="https://www.britannica.com/biography/Robert-Brown-Scottish-botanist">Robert Brown</a> (1773–1858) observed the continuous random movement of small particles of pollen suspended in water under the microscope. This never-ending very irregular and wiggly motion caused by collisions of molecules is known as Brownian motion. It was also found in colloidal particles in chemistry and in many other physical and biological phenomena. Brown pointed out that this irregular motion appeared to obey no known mathematical or physical laws. During the nineteenth century, several attempts had been made to understand the origin of the Brownian motion in microscopic living animals and colloidal particles in chemistry. It has been proposed that the Brownian motion of particles was caused by the electrostatic forces among them, and the constant changes in direction in motion took place due to the influences on particles from molecules of the surrounding medium. For sufficiently small light particles, numerous collisions could have microscopic impact on their random displacements. From a thermodynamic point of view, the thermal excitation of molecules is increased due to increase in temperature of the system, and hence, the physical explanation of Brownian motion is consistent with experimental observations. Fourier suggested that heat propagates almost like a collection of infinitesimal Brownian particles, and hence, his diffusion equation is closely related to Brownian motion. </p> <p> In one of his three famous 1905 papers, Einstein developed a major physical theory of the random motion of the individual particles in Brownian motion. He described the collective motion of Brownian particles of density ρ(<i>x, t</i>) at position <i>x</i> and time <i>t</i> which satisfies the Fourier diffusion equation, <span class="math">\( \rho_t = \kappa\,\rho_{xx} \) </span> with the solution <div class="math"> \[ \rho (x,t) = \frac{1}{\sqrt{4\pi \kappa t}} \,\exp \left\{ - \frac{x^2}{4\kappa t} \right\} , \] </div> where κ is called the <i>diffusion</i> constant which measures how rapidly the distribution spreads out with time. In order to verify Einstein’s theoretical molecular explanation of Brownian motion, <a href="https://en.wikipedia.org/wiki/Jean_Baptiste_Perrin">Jean Perrin</a> (1870--1942) performed experiments and observed that functions describing the paths of Brownian motion are nowhere differentiable, even though they are continuous. </p> <p> Using statistical methods together with heat equation, Einstein showed that if Brownian motion of particles starts at the origin, then after a fixed time <i>t</i>, its position is randomly distributed according to the three-dimensional Gaussian distribution with mean zero and variance κt. In other words, the mathematical description of the Brownian motion of particles in terms of a random process <i>B</i><sub>t</sub> with transitional probability densities is given by the diffusion (or Gaussian) kernel <i>G</i><sub>t</sub> (<i>x, m</i>) in <i>n</i> dimensions in the form <div class="math"> \[ G_t ({\bf x}, \xi ) = \frac{1}{(4\pi t)^{n/2}} \, \exp \left\{ - \frac{({\bf x} - \xi)^2}{4t} \right\} . \] </div> Since diffusion process results from the Brownian motion of particles, the exact solution <i>u</i>(<b>x</b>, <i>t</i>) of the Fourier diffusion equation in ℝ<sup>n</sup> with the initial data <i>u</i>(<b>x</b>,0) = <i>u</i><sub>0</sub>(<b>x</b>) is equal to the expected value of <span class="math">\( u_0 \left( \sqrt{2}\,B_t^x \right) , \) </span> that is, <div class="math"> \[ u ({\bf x}, t) = E \left[ u_0 \left( \sqrt{2}\,B_t^x \right)\right] , \] </div> where <span class="math">\( B_t^x \) </span> is the <i>n</i>-dimensional Brownian motion starting at <b>x</b>. It is further noted that a trajectory of Brownian motion is almost surely non-differentiable. This explains the irregular randomness of the Brownian motion. Thus, the expected value result <span class="math">\( E \left[ u_0 \left( \sqrt{2}\,B_t^x \right)\right] \) </span> describes the relation between the solution of heat equation and probability theory and stochastic processes. Subsequently, this relationship led to the celebrated <a href="https://en.wikipedia.org/wiki/Feynman%E2%80%93Kac_formula">Feynman--Kac formula</a>: <div class="math"> \[ u ({\bf x}, t) = E \left[ u_0 \left( \sqrt{2}\,B_t^x \right) \exp \left\{ - \int_0^t V \left( \sqrt{2}\,B_s^x \right) {\text d}s \right\} \right] , \] </div> which relates the Brownian motion of particles to the solution of the <i>n</i>-dimensional heat equation <div class="math"> \[ u_t ({\bf x}, t) = \nabla^2_n u ({\bf x}, t) - V ({\bf x})\, u ({\bf x}, t) , \] </div> with a given initial data, <i>u</i><sub>0</sub> (<b>x</b>), where ∇²<sub>n</sub> is the <i>n</i>-dimensional Laplacian. Subsequently, the idea of stochastic process or stochastic convergence of probability distributions came into existence when interpreted as modes of convergence of their Fourier transforms. <a href="https://en.wikipedia.org/wiki/Norbert_Wiener">Norbert Wiener</a> (1894--1964) first developed the mathematical theory of probability measure and integration in function space to construct rigorous mathematical models for Brownian motion which was till then understood only heuristically by Einstein and others. He also constructed a probability measure in the space of all continuous functions on <i>C</i>[0, ∞), universally known as the <a href="https://en.wikipedia.org/wiki/Classical_Wiener_space">Wiener measure</a>. In fact, he discovered the systematic mathematical description of Brownian motion known as the <a href="https://en.wikipedia.org/wiki/Wiener_process">Wiener process</a> (a continuous-time stochastic process) in 1923 based on the theory of Fourier series. </p> <p> Subsequently, several great mathematical scientists including <a href="https://en.wikipedia.org/wiki/Andrey_Kolmogorov">Andrey Kolmogorov</a> (1903--1987), <a href="https://en.wikipedia.org/wiki/Paul_L%C3%A9vy_%28mathematician%29">Paul Levy</a> (1886--1971), and <a href="https://en.wikipedia.org/wiki/Kiyosi_It%C3%B4">Kiyoshi Itô</a> (1915--2008) made some major contributions to stochastic analysis for a further understanding of Brownian motion and diffusive processes. Based on Fourier transforms, Levy formulated a new and simpler approach to the general theory of stable probability distributions. From a thermodynamic point of view, Brownian motion is irreversible in the sense that, for a large number of diffusing particles, there exists a smoothing of the probability distributions describing the state of particles as time increases. On the other hand, Itô discovered a powerful method of stochastic calculus and stochastic differential equations for an explicit description of the Markov diffusion process. <p> </p> <div class="chapter"> <p class="custom">Some applications of Fourier analysis to diverse areas </p> <hr color="#cccccc" noshade="noshade" size="1" width="100%" /> </div> <p> We next give some general idea of the broader impact Fourier analysis has had on some selected areas of mathematics, mathematical physics and mathematical economics. Using <a href="https://en.wikipedia.org/wiki/Parseval's_theorem">Parseval’s formula</a> for Fourier series, <a href="https://en.wikipedia.org/wiki/Adolf_Hurwitz">Adolf Hurwitz</a> (1859--1919) gave a very elegant solution of the famous isoperimetric problem in geometry. He proved that the area <i>A</i> of the region enclosed by a simple closed curve Γ in ℝ² of length ℓ satisfies the so-called <a href="https://en.wikipedia.org/wiki/Isoperimetric_inequality">isoperimetric inequality</a> <div class="math"> \[ A \le \frac{\ell^2}{4\pi} , \] </div> where equality holds if Γ is a circle. Hurwitz made some comments praising the use of Fourier series in geometry: ‘Fourier series and analogous expansions intervene very naturally in the general theory of curves and surfaces. In effect, this theory, conceived from the point of view of analysis, deals obviously with the study of arbitrary functions. I was thus led to use Fourier series in several questions of geometry, and I have obtained in this direction a number of results which will be presented in this work. One notes that my considerations form only a beginning of a principal series of researches, which would without doubt give many new results’. </p> <p> Several great mathematicians including <a href="https://en.wikipedia.org/wiki/Bernhard_Riemann">Bernhard Riemann</a> (1826--1866), <a href="https://en.wikipedia.org/wiki/Karl_Weierstrass">Karl Weierstrass</a> (1815--1897) and <a href="https://en.wikipedia.org/wiki/G._H._Hardy">G.H. Hardy</a> (1877--1947) investigated the famous problem of everywhere continuous but nowhere differentiable function in real analysis. In 1861, Riemann suggested that the function <div class="math"> \[ R(x) = \sum_{n\ge 1} \frac{1}{n^2} \,\sin \left( n^2 x \right) , \qquad x \in \mathbb{R} \] </div> is everywhere continuous, but nowhere differentiable in one of his lectures without proof. In order to give a formal proof of the Riemann function, Weierstrass discovered another everywhere continuous but nowhere differentiable function defined by <div class="math"> \[ W(x) = \sum_{n\ge 1} b^n \,\cos \left( a^n x \right) , \] </div> where 0 < <i>b</i> < 1 and <i>a</i> is an integer greater than one. He proved that if <span class="math">\( ab > \left( 1 + \frac{3\pi}{2} \right) , \) </span> the Weierstrass function, <i>W</i>(<i>x</i>) is nowhere differentiable. In 1916, Hardy proved that <i>R</i>(<i>x</i>) is not differentiable at all irrational multiples of π and at some rational multiples of π. In 1969, Joseph Gerver proved that the Riemann function is differentiable at all rational multiples of π of the form (π<i>p</i>/<i>q</i>), where <i>p</i> and <i>q</i> are odd integers and is not differentiable at other points. In the context of Fourier series, there is another function <div class="math"> \[ f(x) = \sum_{n\ge 0} 2^{-n\alpha} \left[ \cos \left( 2^n x \right) + {\bf j} \, \sin \left( 2^n x \right) \right] , \] </div> where 0 < α < 1. This is another continuous but nowhere differentiable function with many vanishing Fourier coefficients. A Fourier series with many missing Fourier coefficients, like <i>W</i>(<i>x</i>) and <i>f</i>(<i>x</i>) is called <a href="http://mathworld.wolfram.com/LacunaryFourierSeries.html">lacunary Fourier series</a>. It can be shown that the real as well as the imaginary part of <i>f</i>(<i>x</i>) are both nowhere differentiable. </p> <p> With ideas from Fourier series, <a href="https://en.wikipedia.org/wiki/Hermann_Weyl">Hermann Weyl</a> (1885--1955) proved his famous criterion which states that a sequence of real numbers <span class="math">\( \left\{ x_n \right\}_1^{\infty} \) </span> in [0, 1] is equidistributed for all integers <i>k</i> ≠ 0 if and only if <div class="math"> \[ \frac{1}{N}\, \sum_{n=1}^N \exp \left\{ 2\pi {\bf j} k\,x_n \right\} \to 0 \qquad \mbox{as} \qquad N \to \infty . \] </div> This means that, to understand the equidistribution properties of a real sequence <span class="math">\( \left\{ x_n \right\}_1^{\infty} , \) </span> it is sufficient to estimate the size of the corresponding Fourier sum <span class="math">\( \sum_{n=1}^N \exp \left\{ 2\pi {\bf j} k\,x_n \right\} . \) </span> For example, it can be shown using the Weyl criterion that the sequence of the fractional parts <span class="math">\( \left\langle n^2 \gamma \right\rangle = n^2 \gamma -\left\lfloor n^2 \gamma \right\rfloor , \) </span> where ⌊A⌋ is the greatest integer ≤A, and is the integral part of <i>n</i>²γ and γ is irrational, is equidistributed in [0,1]. However, it can be proved that the sequence <span class="math">\( \left\{ x_n \right\}_1^{\infty} , \) </span> where <i>x</i><sub>n</sub> is the fractional part of <i>g</i><sup>n</sup>, where <span class="math">\( g = \left( 1 + \sqrt{5} \right) /2 \) </span> is the golden ratio, is not equidistributed in [0, 1]. </p> <p> As a remarkable application of the theory of finite Fourier series, we state Dirichlet’s theorem on primes in arithmetic progression. If <i>a</i> and <i>h</i> are positive integers with no common factor, then the arithmetic progression <div class="math"> \[ a, \quad a+h, \quad a+2h , \quad \cdots \quad a+ nh , \quad \cdots \] </div> contains infinitely many prime numbers. In this theorem, it is the theory of Fourier series on finite Abelian group ℤ<sup>*</sup>(<i>a</i>) that plays a fundamental role in the solution of the problem, where ℤ<sup>*</sup>(<i>a</i>) is defined as the set of all integers modulo a that have multiplicative inverses with multiplication modulo <i>a</i> as the group operation. </p> <p> In 1924, <a href="https://en.wikipedia.org/wiki/Werner_Heisenberg">Werner Heisenberg</a> (1901--1976) first formulated his celebrated uncertainty principle between the position and momentum of a moving particle in quantum mechanics which states that it is not possible to determine the position and momentum of a particle exactly and simultaneously. This principle has an important physical interpretation as an uncertainty of both position and momentum of a particle described by a wave function ψ ∈ 𝔏²(ℝ), so it is a <a href="https://en.wikipedia.org/wiki/Square-integrable_function">square-integrable function</a>. Mathematically, this principle can be formulated in terms of a function and its Fourier transform. If <i>f</i> ∈ 𝔏²(ℝ), then <i>f</i> and its Fourier transform cannot both be essentially localized at <i>x</i>=0, <i>k</i>=0. In other words, it is not possible that the widths of the graphs of <span class="math">\( | f(x) |^2 \) </span> and the square of its Fourier transform <span class="math">\( | F(k) |^2 \) </span> can both be made arbitrarily small. This fact underlines the Heisenberg uncertainty principle in quantum mechanics and the bandwidth theorem in signal analysis. All these can be explained more explicitly as follows: </p> <p> The function <span class="math">\( f_X (x) = |\psi (x)|^2 \) </span> can be interpreted as the probability density for the random position variable <i>X</i> with its probabilistic interpretation <div class="math"> \[ \| \psi \|^2 = \int_{-\infty}^{\infty} |\psi (x) |^2 {\text d} x = \int_{-\infty}^{\infty} f_X (x) \,{\text d} x = 1. \] </div> Similarly, <span class="math">\( f_P (k) = |{\cal F}\left\{ \psi (x) \right\} |^2 = \left\vert \hat{\psi} \right\vert^2 \) </span> is the corresponding probability density function for its momentum <i>P</i>. The Heisenberg uncertainty principle states in the form of an exact inequality that these two probability densities cannot be essentially localized. </p> <p> We then introduce the expectation of the random variable <i>X</i>² by <div class="math"> \[ E \left( X^2 \right) = \int_{-\infty}^{\infty} x^2 f_X (x)\, {\text d} x = \| x\,\psi (x) \|^2 . \] </div> This can be regarded as a measure of the horizontal spread of the wave function ψ(<i>x</i>). Similarly, the spread of its Fourier transform, <span class="math">\( \hat{\psi} (k) \) </span> over the <i>k</i>-axis can be introduced by the integral <div class="math"> \[ \int_{-\infty}^{\infty} k^2 f_P (k)\, {\text d} k = \int_{-\infty}^{\infty} k^2 \left\vert \hat{\psi} (k) \right\vert^2 {\text d} k = \left\| k\,\hat{\psi} (k) \right\|^2 . \] </div> In terms of these two quantities, the uncertainty principle can be formulated in the form of a precise inequality, known as the <a href="https://en.wikipedia.org/wiki/Uncertainty_principle">Heisenberg inequality</a>, <div class="math"> \[ \left\| k\,\hat{\psi} (k) \right\| \cdot \left\| x\,\psi (x) \right\| \ge \frac{1}{2}\,\| \psi \|^2 , \] </div> where the equality holds when ψ(<i>x</i>) is the Gaussian function, that is, <span class="math">\( \psi (x) = C\,\exp \left\{ -a x^2 \right\} , \quad a> 0, \) </span> and <i>C</i> is a constant. </p> <p> It is important to point out new and modern impact of Fourier analysis on mathematical economics and finance. Based on the principle of Brownian motion for random fluctuations of share prices, <a href="https://www.investopedia.com/terms/b/blackscholes.asp">Black and Scholes</a> first formulated a dynamic replication model of price and show that the option price function, <i>C</i>(<i>s, t</i>) satisfies a partial differential equation with an extra drift term in the form <div class="math"> \[ C_t + rs\,C_s + \frac{1}{2}\,\sigma^2 s^2 C_{ss} = r\,C, \qquad 0< t < T, \] </div> with the final boundary condition, <i>C</i>(<i>s, t</i>) = <i>f</i>(<i>s</i>), where <i>r</i> and σ² are constants. This equation is now known as the <a href="https://en.wikipedia.org/wiki/Black%E2%80%93Scholes_equation"><em>Black--Scholes</em></a> (or BS) equation (see <a href="../ch6/bs.html">section</a> in this tutorial). With a suitable change of variables, the BS equation can be reduced to the one-dimensional Fourier diffusion equation. Thus, the Fourier transform of C(s, t) with respect to s can be used to find solution of the BS model equation. Although the BS dynamic replication approach satisfies the principle of market completeness, one problem of this approach is that the price of a derivative does not depend on the drift of the share price. So, there is a need for an alternative approach to derivatives pricing theory, known as the <a href="https://en.wikipedia.org/wiki/Risk-neutral_measure">risk-neutral</a> expectation pricing theory which has also subsequently been developed. Thus, there are two major methods for pricing: the BS replication approach and the risk-neutral expectation approach. Surprisingly, the <a href="https://en.wikipedia.org/wiki/Feynman%E2%80%93Kac_formula">Feynman--Kac formula</a> describing the relation between the solution of the Fourier diffusion equation and probability theory helped combine the above two approaches together so that certain second-order linear partial differential equation can be solved using expectations of diffusive processes. However, there is considerable evidence that logarithm of the share price does not follow the principle of Brownian motion with drift. In particular, when stock market crashes that corresponds to a big jump in the share price leading to a failure of the BS replication strategics. So, the BS model fails to capture some important features of share-price evolution because all paths in Brownian motion are continuous but non-differentiable. This leads to many generalization of the BS model, but most generalizations do not retain the principle of market completeness which means that they give rise to multiple option prices rather than just one price. The analysis of real market data on small-scale share movements reveals that they do not resemble a diffusion process. In fact, they behave more or less like a series of jumps than a Brownian motion. So, rescaling time based on the number of trades that have occurred rather than on calendar time leads to returns that do become approximately normal distribution. Thus, the BS model can be generalized based on trading time. An example of such a model is known as the variance gamma model. In general, the theory of Levy processes can be employed to develop more general theories of price movements for shares and other assets. </p> <p> Since Fourier’s 1822 celebrated work on the theory of heat, the sine and cosine functions have played the fundamental basis functions in the subsequent development of Fourier series and Fourier transforms. However, in many applications, especially in the space-wavenumber (or time-frequency) analysis of a function (or signal), the standard Fourier analysis is not adequate because the Fourier transform of the function (or signal) does not contain any local information. In fact, it excludes the idea of wavenumbers changing with space (or frequencies changing with time), or equivalently, the notion of finding the wavenumber spectrum of a function with space (or frequency spectrum of a signal with time). So, the Fourier transform analysis of functions cannot be used for the analysis of functions (or signals) in a joint space wavenumber (or joint time-frequency) domain. So, in 1980s, a whole new idea of wavelets has been introduced as a family of functions constructed from translations and dilations of a single function called the <a href="https://www.sciencedirect.com/topics/computer-science/mother-wavelet">‘mother wavelet’</a> ψ(<i>x</i>). They are defined by <div class="math"> \[ \psi_{a,b} (x) = \frac{1}{\sqrt{|a|}}\,\psi \left( \frac{x-b}{a} \right) , \qquad a,b \in \mathbb{R} , \] </div> where <i>a</i> is called a <em>scaling parameter</em> which measures the degree of compression or scale, and <i>b</i> is a <em>translation parameter</em> which determines the space location of a wavelet. If ψ ∈ 𝔏² (ℝ) and ψ<sub>a,b</sub>(<i>x</i>) is given by the latter equation, the continuous wavelet transform is defined by <div class="math"> \[ \left( W_{\psi} (f)\right) (a,b) = \int_{-\infty}^{\infty} \overline{\psi_{a,b} (x)}\,f(x)\,{\text d}x , \] </div> where ψ<sub>a,b</sub>(<i>x</i>) plays the same role as the kernel exp(<b>j</b><i>kx</i>) in the Fourier transform. Similar the Fourier transform, the <a href="https://en.wikipedia.org/wiki/Continuous_wavelet_transform">continuous wavelet transform</a> is linear, and it is not a single transform, but any transform obtained in this way. The inverse wavelet transform can be defined so that <i>f</i>(<i>x</i>) can be reconstructed by a double integral formula. It is now well-known that, analogous to the Fourier analysis, the wavelet analysis has provided a famous new method for decomposing a function or a signal in 1980s. It is more useful and better suited to a wide variety of problems in applied and computational mathematics and also to many questions in real and harmonic analysis. The most significant contributions of both discrete and continuous wavelets can be found in signal processing and digital image compression. <p> </p> <div class="chapter"> <p class="custom">Historical development of Fourier series </p> <hr color="#cccccc" noshade="noshade" size="1" width="100%" /> </div> <p> Historically, in 1747, <a href="https://en.wikipedia.org/wiki/Jean_le_Rond_d'Alembert">d’Alembert</a> first formulated the wave equation with the first solution of a partial differential equation. In fact, Fourier series originated from the early study of the initial-boundary problem of the one-dimensional wave equation for vibration of an elastic string fixed at its ends in the eighteenth century and of the Fourier heat equation of a solid body in 1800s. <a href="https://en.wikipedia.org/wiki/Leonhard_Euler">Euler</a>, d’Alembert and <a href="https://en.wikipedia.org/wiki/Daniel_Bernoulli">Daniel Bernoulli</a> made earliest attempts to solve the initial-boundary value problem of the classical one-dimensional wave equation <div class="math"> \[ u_{tt} = c^2 u_{xx} , \qquad 0 < x < \ell , \quad t> 0, \] </div> where <i>u</i> = <i>u</i>(<i>x,t</i>) is the displacement of the string of length ℓ at the position <i>x</i> and time <i>t</i>, <i>c</i>² = (<i>T/ρ</i>), <i>T</i> is the tension and ρ s the mass density of the string. <p> The initial and boundary conditions of this problem are given by <div class="math"> \[ u(x,0) = d(x) \qquad\mbox{and}\qquad u_t (x,0) = v(x), \qquad 0 < x < \ell , \] </div> <div class="math"> \[ u(0,y) = 0 = u(\ell , t) , \qquad t> 0, \] </div> where <i>d</i>(<i>x</i>) is the initial displacement and <i>v</i>(<i>x</i>) is the initial velocity of the string. <p> Both Euler and d’Alembert made the first attempt to solve the initial value problem as the traveling wave solution in the functional form <div class="math"> \[ u(x,y) = \phi (x-ct) + \psi (x+ct) , \] </div> where φ(x-ct) represents a wave traveling to the right with constant velocity <i>c</i> and φ(x+ct) represents a wave moving to the left with constant speed <i>c</i>. The above solution is the general solution of the wave equation representing the superposition of two traveling wave solutions. <p> On the other hand, Daniel Bernoulli derived the solution of the wave equation with the boundary conditions satisfied in an infinite trigonometric series form <div class="math"> \[ u(x,y) = \sum_{n\ge 1} \left[ a_n \cos \left( \frac{n\pi ct}{\ell} \right) + b_n \sin \left( \frac{n\pi ct}{\ell} \right) \right] \sin \left( \frac{n\pi x}{\ell} \right) , \] </div> where <i>a</i><sub>n</sub> and <i>b</i><sub>n</sub> are to be determined from the initial conditions <div class="math"> \begin{align*} d(x) &= u(x,0) = \sum_{n\ge 1} a_n \sin \left( \frac{n\pi x}{\ell} \right) , \\ v(x) &= u_t (x,0) = \sum_{n\ge 1} b_n \left( \frac{n\pi c}{\ell} \right) \sin \left( \frac{n\pi x}{\ell} \right) . \end{align*} </div> This means that "arbitrary" functions <i>d</i>(<i>x</i>) and <i>v</i>(<i>x</i>) have been expanded in trigonometric (or Fourier) series. Since the above equations are infinite series, the question of convergence arose. But many mathematicians including Euler, d’Alembert and <a href="https://en.wikipedia.org/wiki/Pierre-Simon_Laplace">Laplace</a> were not entirely convinced about the existence of Bernoulli’s Fourier series solution. This led to a serious controversy as the solution could hold if an arbitrary function could be represented in Fourier series. In 1759, a young mathematical scientist, Lagrange joined in this controversial debate on the various forms of the solution of the vibrating string problem and published a paper which was criticized by both Euler and d’Alembert. In order to resolve this controversy, <a href="https://en.wikipedia.org/wiki/Joseph-Louis_Lagrange">Lagrange</a> gave an entirely new treatment of the problem. From 1768, d’Alembert wrote a series of notes in which he attacked Euler’s views on the solution of the problem. In 1779, Laplace entered into this debate and supported d’Alembert views. They all discovered solutions in several different forms and the merit of their solutions and their relations among these solutions extending over more than 25 years. Their major concerns were whether an arbitrary function can be expanded in Fourier series and convergence of such series. Unfortunately, the outcome of this controversial issues remained inconclusive for a long period of time. </p> <p> In 1807, it was Joseph Fourier who resolved this long-standing controversy by discovering the solution of the one-dimensional heat equation in terms of Fourier series. Indeed, he formulated and solved the initial-boundary value of the one-dimensional heat equation in a rod of length <div class="math"> \[ u_{t} = \kappa\, u_{xx} , \qquad 0 < x < \ell , \quad t> 0, \] </div> with the initial and boundary conditions <div class="math"> \[ u(x,0) = f(x) , \qquad 0 < x < \ell , \] </div> <div class="math"> \[ u(0,y) = 0 = u(\ell , t) , \qquad t> 0. \] </div> Using the method of separation of variables, Fourier obtained the most general and rapidly convergent solution for the temperature distribution <i>u</i>(<i>x, t</i>) in the form <div class="math"> \[ u(x,y) = \sum_{n\ge 1} c_n \exp \left\{ -\left( \frac{n\pi x}{\ell} \right)^2 \kappa t \right\} \sin \left( \frac{n\pi x}{\ell} \right) , \] </div> where the initial condition implies that <div class="math"> \[ f(x) = \sum_{n\ge 1} c_n \sin \left( \frac{n\pi x}{\ell} \right) . \] </div> The upshot of Fourier’s analysis led to expanding an arbitrary function in an infinite series of sine functions and then finding the infinite number of coefficients <div class="math"> \[ c_n = \frac{2}{\ell} \int_0^{\ell} f(x)\,\sin \left( \frac{n\pi x}{\ell} \right) {\text d} x, \qquad n= 1,2,\ldots . \] </div> <p> </p> <button type="button" class="btn btn-info collapse_panel" data-toggle="collapse" data-target="#ref"> References </button> <div id="ref" class="collapse"> <ol> <li> Debnath, L., A short biography of Joseph Fourier and historical development of Fourier series and Fourier transforms, <em>International Journal of Mathematical Education in Science and Technology</em>, 2012, Vol. 43, No. 5, pp. 58900612; https://doi.org/10.1080/0020739X.2011.633712 </li> <li> Gerver, J., The Differentiability of the Riemann Function at Certain Rational Multiples of π, <em>American Journal of Mathematics</em>, 1970, Vol. 92, No. 1, pp. 33--55. http://www.jstor.org/stable/2373496 </li> <li> Gibbs, J. W., Fourier Series, <em>Nature</em>, 1899, Vol. 59, 200 and 606. </li> <li> </li> <li> </li> <li> </li> </ol> </div> <p> </p> <p>Return to <a href="../index.html">Mathematica</a> page <br /> Return to the <a href="../../index.html">main page</a> (APMA0340) <br /> Return to the <a href="../ch1/part1.html">Part 1</a> Matrix Algebra <br /> Return to the <a href="../ch2/part2.html">Part 2</a> Linear Systems of Ordinary Differential Equations <br /> Return to the <a href="../ch3/part3.html">Part 3</a> Non-linear Systems of Ordinary Differential Equations <br /> Return to the <a href="../ch4/part4.html">Part 4 </a> Numerical Methods <br /> Return to the <a href="../ch5/part5.html">Part 5 </a> Fourier Series <br /> Return to the <a href="../ch6/part6.html">Part 6 </a> Partial Differential Equations<br /> Return to the <a href="../ch7/part7.html">Part 7 </a> Special Functions </p> <p> </p> </body> </html>