{ "hash": "7de6019c5ee744a1d63f36de9098cd63", "result": { "markdown": "# Area under a curve\n\n\n\nThis section uses these add-on packages:\n\n``` {.julia .cell-code}\nusing CalculusWithJulia\nusing Plots\nusing QuadGK\nusing Roots\n```\n\n\n\n\n---\n\n\nThe question of area has long fascinated human culture. As children, we learn early on the formulas for the areas of some geometric figures: a square is $b^2$, a rectangle $b\\cdot h$ a triangle $1/2 \\cdot b \\cdot h$ and for a circle, $\\pi r^2$. The area of a rectangle is often the intuitive basis for illustrating multiplication. The area of a triangle has been known for ages. Even complicated expressions, such as [Heron's](http://tinyurl.com/mqm9z) formula which relates the area of a triangle with measurements from its perimeter have been around for 2000 years. The formula for the area of a circle is also quite old. Wikipedia dates it as far back as the [Rhind](http://en.wikipedia.org/wiki/Rhind_Mathematical_Papyrus) papyrus for 1700 BC, with the approximation of $256/81$ for $\\pi$.\n\n\nThe modern approach to area begins with a non-negative function $f(x)$ over an interval $[a,b]$. The goal is to compute the area under the graph. That is, the area between $f(x)$ and the $x$-axis between $a \\leq x \\leq b$.\n\n\nFor some functions, this area can be computed by geometry, for example, here we see the area under $f(x)$ is just $1$, as it is a triangle with base $2$ and height $1$:\n\n::: {.cell hold='true' execution_count=4}\n``` {.julia .cell-code}\nf(x) = 1 - abs(x)\nplot(f, -1, 1)\nplot!(zero)\n```\n\n::: {.cell-output .cell-output-display execution_count=5}\n![](area_files/figure-html/cell-5-output-1.svg){}\n:::\n:::\n\n\nSimilarly, we know this area is also $1$, it being a square:\n\n::: {.cell hold='true' execution_count=5}\n``` {.julia .cell-code}\nf(x) = 1\nplot(f, 0, 1)\nplot!(zero)\n```\n\n::: {.cell-output .cell-output-display execution_count=6}\n![](area_files/figure-html/cell-6-output-1.svg){}\n:::\n:::\n\n\nThis one, is simply $\\pi/2$, it being half a circle of radius $1$:\n\n::: {.cell hold='true' execution_count=6}\n``` {.julia .cell-code}\nf(x) = sqrt(1 - x^2)\nplot(f, -1, 1)\nplot!(zero)\n```\n\n::: {.cell-output .cell-output-display execution_count=7}\n![](area_files/figure-html/cell-7-output-1.svg){}\n:::\n:::\n\n\nAnd this area can be broken into a sum of the area of square and the area of a rectangle, or $1 + 1/2$:\n\n::: {.cell hold='true' execution_count=7}\n``` {.julia .cell-code}\nf(x) = x > 1 ? 2 - x : 1.0\nplot(f, 0, 2)\nplot!(zero)\n```\n\n::: {.cell-output .cell-output-display execution_count=8}\n![](area_files/figure-html/cell-8-output-1.svg){}\n:::\n:::\n\n\nBut what of more complicated areas? Can these have their area computed?\n\n\n## Approximating areas\n\n\nIn a previous section, we saw this animation:\n\n::: {.cell cache='true' hold='true' execution_count=8}\n\n::: {.cell-output .cell-output-display execution_count=9}\n```{=html}\n
\"A\n

The first triangle has area \\(1/2\\), the second has area \\(1/8\\), then \\(2\\) have area \\((1/8)^2\\), \\(4\\) have area \\((1/8)^3\\), ... With some algebra, the total area then should be \\(1/2 \\cdot (1 + (1/4) + (1/4)^2 + \\cdots) = 2/3\\).

\n
\n
\n
\n```\n:::\n:::\n\n\nThis illustrates a method of [Archimedes](http://en.wikipedia.org/wiki/The_Quadrature_of_the_Parabola) to compute the area contained in a parabola using the method of exhaustion. Archimedes leveraged a fact he discovered relating the areas of triangle inscribed with parabolic segments to create a sum that could be computed.\n\n\nThe pursuit of computing areas persisted. The method of computing area by finding a square with an equivalent area was known as *quadrature*. Over the years, many figures had their area computed, for example, the area under the graph of the [cycloid](http://en.wikipedia.org/wiki/Cycloid) (...Galileo tried empirically to find this using a tracing on sheet metal and a scale).\n\n\nHowever, as areas of geometric objects were replaced by the more general question of area related to graphs of functions, a more general study was called for.\n\n\nOne such approach is illustrated in this figure due to Beeckman from 1618 (from [Bressoud](http://www.math.harvard.edu/~knill/teaching/math1a_2011/exhibits/bressoud/))\n\n\n\nBeeckman actually did more than find the area. He generalized the relationship of rate $\\times$ time $=$ distance. The line was interpreting a velocity, the \"squares\", then, provided an approximate distance traveled when the velocity is taken as a constant on the small time interval. Then the distance traveled can be approximated by a smaller quantity - just add the area of the rectangles squarely within the desired area ($6+16+6$) - and a larger quantity - by including all rectangles that have a portion of their area within the desired area ($10 + 16 + 10$). Beeckman argued that the error vanishes as the rectangles get smaller.\n\n\nAdding up the smaller \"squares\" can be a bit more efficient if we were to add all those in a row, or column at once. We would then add the areas of a smaller number of rectangles. For this curve, the two approaches are basically identical. For other curves, identifying which squares in a row would be added is much more complicated (though useful), but for a curve generated by a function, identifying which \"squares\" go in a rectangle is quite easy, in fact we can see the rectangle's area will be a base given by that of the squares, and height depending on the function.\n\n\n### Adding rectangles\n\n\nThe idea of the Riemann sum then is to approximate the area under the curve by the area of well-chosen rectangles in such a way that as the bases of the rectangles get smaller (hence adding more rectangles) the error in approximation vanishes.\n\n\nDefine a partition of $[a,b]$ to be a selection of points $a = x_0 < x_1 < \\cdots < x_{n-1} < x_n = b$. The norm of the partition is the largest of all the differences $\\lvert x_i - x_{i-1} \\rvert$. For a partition, consider an arbitrary selection of points $c_i$ satisfying $x_{i-1} \\leq c_i \\leq x_{i}$, $1 \\leq i \\leq n$. Then the following is a **Riemann sum**:\n\n\n\n$$\nS_n = f(c_1) \\cdot (x_1 - x_0) + f(c_2) \\cdot (x_2 - x_1) + \\cdots + f(c_n) \\cdot (x_n - x_{n-1}).\n$$\n\n\nClearly for a given partition and choice of $c_i$, the above can be computed. Each term $f(c_i)\\cdot(x_i-x_{i-1})$ can be visualized as the area of a rectangle with base spanning from $x_{i-1}$ to $x_i$ and height given by the function value at $c_i$. The following visualizes left Riemann sums for different values of $n$ in a way that makes Beekman's intuition plausible – that as the number of rectangles gets larger, the approximate sum will get closer to the actual area.\n\n::: {.cell hold='true' execution_count=10}\n\n::: {.cell-output .cell-output-display execution_count=11}\n```{=html}\n
\"A\n

Illustration of left Riemann sum for increasing \\(n\\) values

\n
\n
\n
\n```\n:::\n:::\n\n\nTo successfully compute a good approximation for the area, we would need to choose $c_i$ and the partition so that a formula can be found to express the dependence on the size of the partition.\n\n\nFor Archimedes' problem - finding the area under $f(x)=x^2$ between $0$ and $1$ - if we take as a partition $x_i = i/n$ and $c_i = x_i$, then the above sum becomes:\n\n\n\n$$\n\\begin{align*}\nS_n &= f(c_1) \\cdot (x_1 - x_0) + f(c_2) \\cdot (x_2 - x_1) + \\cdots + f(c_n) \\cdot (x_n - x_{n-1})\\\\\n&= (x_1)^2 \\cdot \\frac{1}{n} + (x_2)^2 \\cdot \\frac{1}{n} + \\cdot + (x_n)^2 \\cdot \\frac{1}{n}\\\\\n&= 1^2 \\cdot \\frac{1}{n^3} + 2^2 \\cdot \\frac{1}{n^3} + \\cdots + n^2 \\cdot \\frac{1}{n^3}\\\\\n&= \\frac{1}{n^3} \\cdot (1^2 + 2^2 + \\cdots + n^2) \\\\\n&= \\frac{1}{n^3} \\cdot \\frac{n\\cdot(n-1)\\cdot(2n+1)}{6}.\n\\end{align*}\n$$\n\n\nThe latter uses a well-known formula for the sum of squares of the first $n$ natural numbers.\n\n\nWith this expression, it is readily seen that as $n$ gets large this value gets close to $2/6 = 1/3$.\n\n\n:::{.callout-note}\n## Note\nThe above approach, like Archimedes', ends with a limit being taken. The answer comes from using a limit to add a big number of small values. As with all limit questions, worrying about whether a limit exists is fundamental. For this problem, we will see that for the general statement there is a stretching of the formal concept of a limit.\n\n:::\n\n---\n\nThere is a more compact notation to $x_1 + x_2 + \\cdots + x_n$, this using the *summation notation* or capital sigma. We have:\n\n\n\n$$\n\\Sigma_{i = 1}^n x_i = x_1 + x_2 + \\cdots + x_n\n$$\n\n\nThe notation includes three pieces of information:\n\n\n * The $\\Sigma$ is an indication of a sum\n * The ${i=1}$ and $n$ sub- and superscripts indicate the range to sum over.\n * The term $x_i$ is a general term describing the $i$th entry, where it is understood that $i$ is just some arbitrary indexing value.\n\n\nWith this notation, a Riemann sum can be written as $\\Sigma_{i=1}^n f(c_i)(x_i-x_{i-1})$.\n\n\n### Other sums\n\n\nThe choice of the $c_i$ will give different answers for the approximation, though for an integrable function these differences will vanish in the limit. Some common choices are:\n\n\n * Using the right hand endpoint of the interval $[x_{i-1}, x_i]$ giving the right-Riemann sum, $R_n$.\n * The choice $c_i = x_{i-1}$ gives the left-Riemann sum, $L_n$.\n * The choice $c_i = (x_i + x_{i-1})/2$ is the midpoint rule, $M_n$.\n * If the function is continuous on the closed subinterval $[x_{i-1}, x_i]$, then it will take on its minimum and maximum values. By the extreme value theorem, we could take $c_i$ to correspond to either the maximum or the minimum. These choices give the \"upper Riemann-sums\" and \"lower Riemann-sums\".\n\n\nThe choice of partition can also give different answers. A common choice is to break the interval into $n+1$ equal-sized pieces. With $\\Delta = (b-a)/n$, these pieces become the arithmetic sequence $a = a + 0 \\cdot \\Delta < a + 1 \\cdot \\Delta < a + 2 \\cdot \\Delta < \\cdots a + n \\cdots < \\Delta = b$ with $x_i = a + i (b-a)/n$. (The `range(a, b, length=n+1)` command will compute these.) An alternate choice made below for one problem is to use a geometric progression:\n\n\n\n$$\na = a(1+\\alpha)^0 < a(1+\\alpha)^1 < a (1+\\alpha)^2 < \\cdots < a (1+\\alpha)^n = b.\n$$\n\n\nThe general statement allows for any partition such that the largest gap goes to $0$.\n\n---\n\n\nRiemann sums weren't named after Riemann because he was the first to approximate areas using rectangles. Indeed, others had been using even more efficient ways to compute areas for centuries prior to Riemann's work. Rather, Riemann put the definition of the area under the curve on a firm theoretical footing with the following theorem which gives a concrete notion of what functions are integrable:\n\n\n> **Riemann Integral**: A function $f$ is Riemann integrable over the interval $[a,b]$ and its integral will have value $V$ provided for every $\\epsilon > 0$ there exists a $\\delta > 0$ such that for any partition $a =x_0 < x_1 < \\cdots < x_n=b$ with $\\lvert x_i - x_{i-1} \\rvert < \\delta$ and for any choice of points $x_{i-1} \\leq c_i \\leq x_{i}$ this is satisfied:\n>\n> $$\n> \\lvert \\sum_{i=1}^n f(c_i)(x_{i} - x_{i-1}) - V \\rvert < \\epsilon.\n> $$\n>\n> When the integral exists, it is written $V = \\int_a^b f(x) dx$.\n\n\n\n:::{.callout-note}\n## History note\nThe expression $V = \\int_a^b f(x) dx$ is known as the *definite integral* of $f$ over $[a,b]$. Much earlier than Riemann, Cauchy had defined the definite integral in terms of a sum of rectangular products beginning with $S=(x_1 - x_0) f(x_0) + (x_2 - x_1) f(x_1) + \\cdots + (x_n - x_{n-1}) f(x_{n-1})$ (the left Riemann sum). He showed the limit was well defined for any continuous function. Riemann's formulation relaxes the choice of partition and the choice of the $c_i$ so that integrability can be better understood.\n\n:::\n\n### Some immediate consequences\n\n\nThe following formulas are consequences when $f(x)$ is integrable. These mostly follow through a judicious rearranging of the approximating sums.\n\n\nThe area is $0$ when there is no width to the interval to integrate over:\n\n\n> $$\n> \\int_a^a f(x) dx = 0.\n> $$\n\n\n\nEven our definition of a partition doesn't really apply, as we assume $a < b$, but clearly if $a=x_0=x_n=b$ then our only\"approximating\" sum could be $f(a)(b-a) = 0$.\n\n\nThe area under a constant function is found from the area of rectangle, a special case being $c=0$ yielding $0$ area:\n\n\n> $$\n> \\int_a^b c dx = c \\cdot (b-a).\n> $$\n\n\n\nFor any partition of $a < b$, we have $S_n = c(x_1 - x_0) + c(x_2 -x_1) + \\cdots + c(x_n - x_{n-1})$. By factoring out the $c$, we have a *telescoping sum* which means the sum simplifies to $S_n = c(x_n-x_0) = c(b-a)$. Hence any limit must be this constant value.\n\n\nScaling the $y$ axis by a constant can be done before or after computing the area:\n\n\n> $$\n> \\int_a^b cf(x) dx = c \\int_a^b f(x) dx.\n> $$\n\n\n\nLet $a=x_0 < x_1 < \\cdots < x_n=b$ be any partition. Then we have $S_n= cf(c_1)(x_1-x_0) + \\cdots + cf(c_n)(x_n-x_0)$ $=$ $c\\cdot\\left[ f(c_1)(x_1 - x_0) + \\cdots + f(c_n)(x_n - x_0)\\right]$. The \"limit\" of the left side is $\\int_a^b c f(x) dx$. The \"limit\" of the right side is $c \\cdot \\int_a^b f(x)$. We call this a \"sketch\" as a formal proof would show that for any $\\epsilon$ we could choose a $\\delta$ so that any partition with norm $\\delta$ will yield a sum less than $\\epsilon$. Here, then our \"any\" partition would be one for which the $\\delta$ on the left hand side applies. The computation shows that the same $\\delta$ would apply for the right hand side when $\\epsilon$ is the same.\n\n\nThe area is invariant under shifts left or right.\n\n\n> $$\n> \\int_a^b f(x - c) dx = \\int_{a-c}^{b-c} f(x) dx.\n> $$\n\n\n\nAny partition $a =x_0 < x_1 < \\cdots < x_n=b$ is related to a partition of $[a-c, b-c]$ through $a-c < x_0-c < x_1-c < \\cdots < x_n - c = b-c$. Let $d_i=c_i-c$ denote this partition, then we have:\n\n\n\n$$\nf(c_1 -c) \\cdot (x_1 - x_0) + f(c_2 -c) \\cdot (x_2 - x_1) + \\cdots + f(c_n -c) \\cdot (x_n - x_{n-1}) =\nf(d_1) \\cdot(x_1-c - (x_0-c)) + f(d_2) \\cdot(x_2-c - (x_1-c)) + \\cdots + f(d_n) \\cdot(x_n-c - (x_{n-1}-c)).\n$$\n\n\nThe left side will have a limit of $\\int_a^b f(x-c) dx$ the right would have a \"limit\" of $\\int_{a-c}^{b-c}f(x)dx$.\n\n\nSimilarly, reflections don't effect the area under the curve, they just require a new parameterization:\n\n\n> $$\n> \\int_a^b f(x) dx = \\int_{-b}^{-a} f(-x) dx\n> $$\n\n\n\nThe scaling operation $g(x) = f(cx)$ has the following:\n\n\n> $$\n> \\int_a^b f(c\\cdot x) dx = \\frac{1}{c} \\int_{ca}^{cb}f(x) dx\n> $$\n\n\n\nThe scaling operation shifts $a$ to $ca$ and $b$ to $cb$ so the limits of integration make sense. However, the area stretches by $c$ in the $x$ direction, so must contract by $c$ in the $y$ direction to stay in balance. Hence the factor of $1/c$.\n\n\nCombining two operations above, the operation $g(x) = \\frac{1}{h}f(\\frac{x-c}{h})$ will leave the area between $a$ and $b$ under $g$ the same as the area under $g$ between $(a-c)/h$ and $(b-c)/h$.\n\n\n---\n\nThe area between $a$ and $b$ can be broken up into the sum of the area between $a$ and $c$ and that between $c$ and $b$.\n\n\n> $$\n> \\int_a^b f(x) dx = \\int_a^c f(x) dx + \\int_c^b f(x) dx.\n> $$\n\n\n\nFor this, suppose we have a partition for both the integrals on the right hand side for a given $\\epsilon/2$ and $\\delta$. Combining these into a partition of $[a,b]$ will mean $\\delta$ is still the norm. The approximating sum will combine to be no more than $\\epsilon/2 + \\epsilon/2$, so for a given $\\epsilon$, this $\\delta$ applies.\n\n\nThis is due to the area on the left and right of $0$ being equivalent.\n\n\nThe \"reversed\" area is the same, only accounted for with a minus sign.\n\n\n> $$\n> \\int_a^b f(x) dx = -\\int_b^a f(x) dx.\n> $$\n\n\n\nA consequence of the last few statements is:\n\n\n> If $f(x)$ is an even function, then $\\int_{-a}^a f(x) dx = 2 \\int_0^a f(x) dx$. If $f(x)$ is an odd function, then $\\int_{-a}^a f(x) dx = 0$.\n\n\n\nIf $g$ bounds $f$ then the area under $g$ will bound the area under $f$, in particular if $f(x)$ is non negative, so will the area under $f$ be non negative for any $a < b$. (This assumes that $g$ and $f$ are integrable.)\n\n\n> If $0 \\leq f(x) \\leq g(x)$ then $\\int_a^b f(x) dx \\leq \\int_a^b g(x) dx.$\n\n\n\nFor any partition of $[a,b]$ and choice of $c_i$, we have the term-by-term bound $f(c_i)(x_i-x_{i-1}) \\leq g(c_i)(x_i-x_{i-1})$ So any sequence of partitions that converges to the limits will have this inequality maintained for the sum.\n\n\n### Some known integrals\n\n\nUsing the definition, we can compute a few definite integrals:\n\n\n> $$\n> \\int_a^b c dx = c \\cdot (b-a).\n> $$\n\n\n\n> $$\n> \\int_a^b x dx = \\frac{b^2}{2} - \\frac{a^2}{2}.\n> $$\n\n\n\nThis is just the area of a trapezoid with heights $a$ and $b$ and side length $b-a$, or $1/2 \\cdot (b + a) \\cdot (b - a)$. The right sum would be:\n\n\n\n$$\n\\begin{align*}\nS &= x_1 \\cdot (x_1 - x_0) + x_2 \\cdot (x_2 - x_1) + \\cdots x_n \\cdot (x_n - x_{n-1}) \\\\\n&= (a + 1\\frac{b-a}{n}) \\cdot \\frac{b-a}{n} + (a + 2\\frac{b-a}{n}) \\cdot \\frac{b-a}{n} + \\cdots (a + n\\frac{b-a}{n}) \\cdot \\frac{b-a}{n}\\\\\n&= n \\cdot a \\cdot (\\frac{b-a}{n}) + (1 + 2 + \\cdots n) \\cdot (\\frac{b-a}{n})^2 \\\\\n&= n \\cdot a \\cdot (\\frac{b-a}{n}) + \\frac{n(n+1)}{2} \\cdot (\\frac{b-a}{n})^2 \\\\\n& \\rightarrow a \\cdot(b-a) + \\frac{(b-a)^2}{2} \\\\\n&= \\frac{b^2}{2} - \\frac{a^2}{2}.\n\\end{align*}\n$$\n\n\n> $$\n> \\int_a^b x^2 dx = \\frac{b^3}{3} - \\frac{a^3}{3}.\n> $$\n\n\n\nThis is similar to the Archimedes case with $a=0$ and $b=1$ shown above.\n\n\n> $$\n> \\int_a^b x^k dx = \\frac{b^{k+1}}{k+1} - \\frac{a^{k+1}}{k+1},\\quad k \\neq -1\n> $$\n>\n> .\n\n\n\nCauchy showed this using a *geometric series* for the partition, not the arithmetic series $x_i = a + i (b-a)/n$. The series defined by $1 + \\alpha = (b/a)^{1/n}$, then $x_i = a \\cdot (1 + \\alpha)^i$. Here the bases $x_{i+1} - x_i$ simplify to $x_i \\cdot \\alpha$ and $f(x_i) = (a\\cdot(1+\\alpha)^i)^k = a^k (1+\\alpha)^{ik}$, or $f(x_i)(x_{i+1}-x_i) = a^{k+1}\\alpha[(1+\\alpha)^{k+1}]^i$, so, using $u=(1+\\alpha)^{k+1}=(b/a)^{(k+1)/n}$, $f(x_i) \\cdot(x_{i+1} - x_i) = a^{k+1}\\alpha u^i$. This gives\n\n\n\n$$\n\\begin{align*}\nS &= a^{k+1}\\alpha u^0 + a^{k+1}\\alpha u^1 + \\cdots + a^{k+1}\\alpha u^{n-1}\n&= a^{k+1} \\cdot \\alpha \\cdot (u^0 + u^1 + \\cdot u^{n-1}) \\\\\n&= a^{k+1} \\cdot \\alpha \\cdot \\frac{u^n - 1}{u - 1}\\\\\n&= (b^{k+1} - a^{k+1}) \\cdot \\frac{\\alpha}{(1+\\alpha)^{k+1} - 1} \\\\\n&\\rightarrow \\frac{b^{k+1} - a^{k+1}}{k+1}.\n\\end{align*}\n$$\n\n\n> $$\n> \\int_a^b x^{-1} dx = \\log(b) - \\log(a), \\quad (0 < a < b).\n> $$\n\n\n\nAgain, Cauchy showed this using a geometric series. The expression $f(x_i) \\cdot(x_{i+1} - x_i)$ becomes just $\\alpha$. So the approximating sum becomes:\n\n\n\n$$\nS = f(x_0)(x_1 - x_0) + f(x_1)(x_2 - x_1) + \\cdots + f(x_{n-1}) (x_n - x_{n-1}) = \\alpha + \\alpha + \\cdots \\alpha = n\\alpha.\n$$\n\n\nBut, letting $x = 1/n$, the limit above is just the limit of\n\n\n\n$$\n\\lim_{x \\rightarrow 0+} \\frac{(b/a)^x - 1}{x} = \\log(b/a) = \\log(b) - \\log(a).\n$$\n\n\n(Using L'Hopital's rule to compute the limit.)\n\n\nCertainly other integrals could be computed with various tricks, but we won't pursue this. There is another way to evaluate integrals using the forthcoming Fundamental Theorem of Calculus.\n\n\n### Some other consequences\n\n\n * The definition is defined in terms of any partition with its norm bounded by $\\delta$. If you know a function $f$ is Riemann integrable, then it is enough to consider just a regular partition $x_i = a + i \\cdot (b-a)/n$ when forming the sums, as was done above. It is just that showing a limit for just this particular type of partition would not be sufficient to prove Riemann integrability.\n * The choice of $c_i$ is arbitrary to allow for maximum flexibility. The Darboux integrals use the maximum and minimum over the subinterval. It is sufficient to prove integrability to show that the limit exists with just these choices.\n * Most importantly,\n\n\n> A continuous function on $[a,b]$ is Riemann integrable on $[a,b]$.\n\n\n\nThe main idea behind this is that the difference between the maximum and minimum values over a partition gets small. That is if $[x_{i-1}, x_i]$ is like $1/n$ is length, then the difference between the maximum of $f$ over this interval, $M$, and the minimum, $m$ over this interval will go to zero as $n$ gets big. That $m$ and $M$ exists is due to the extreme value theorem, that this difference goes to $0$ is a consequence of continuity. What is needed is that this value goes to $0$ at the same rate – no matter what interval is being discussed – is a consequence of a notion of uniform continuity, a concept discussed in advanced calculus, but which holds for continuous functions on closed intervals. Armed with this, the Riemann sum for a general partition can be bounded by this difference times $b-a$, which will go to zero. So the upper and lower Riemann sums will converge to the same value.\n\n\n * A \"jump\", or discontinuity of the first kind, is a value $c$ in $[a,b]$ where $\\lim_{x \\rightarrow c+} f(x)$ and $\\lim_{x \\rightarrow c-}f(x)$ both exist, but are not equal. It is true that a function that is not continuous on $I=[a,b]$, but only has discontinuities of the first kind on $I$ will be Riemann integrable on $I$.\n\n\nFor example, the function $f(x) = 1$ for $x$ in $[0,1]$ and $0$ otherwise will be integrable, as it is continuous at all but two points, $0$ and $1$, where it jumps.\n\n\n * Some functions can have infinitely many points of discontinuity and still be integrable. The example of $f(x) = 1/q$ when $x=p/q$ is rational, and $0$ otherwise is often used as an example.\n\n\n## Numeric integration\n\n\nThe Riemann sum approach gives a method to approximate the value of a definite integral. We just compute an approximating sum for a large value of $n$, so large that the limiting value and the approximating sum are close.\n\n\nTo see the mechanics, let's again return to Archimedes' problem and compute $\\int_0^1 x^2 dx$.\n\n\nLet us fix some values:\n\n::: {.cell execution_count=11}\n``` {.julia .cell-code}\na, b = 0, 1\nf(x) = x^2\n```\n\n::: {.cell-output .cell-output-display execution_count=12}\n```\nf (generic function with 1 method)\n```\n:::\n:::\n\n\nThen for a given $n$ we have some steps to do: create the partition, find the $c_i$, multiply the pieces and add up. Here is one way to do all this:\n\n::: {.cell execution_count=12}\n``` {.julia .cell-code}\nn = 5\nxs = a:(b-a)/n:b # also range(a, b, length=n)\ndeltas = diff(xs) # forms x2-x1, x3-x2, ..., xn-xn-1\ncs = xs[1:end-1] # finds left-hand end points. xs[2:end] would be right-hand ones.\n```\n\n::: {.cell-output .cell-output-display execution_count=13}\n```\n0.0:0.2:0.8\n```\n:::\n:::\n\n\nNow to multiply the values. We want to sum the product `f(cs[i]) * deltas[i]`, here is one way to do so:\n\n::: {.cell execution_count=13}\n``` {.julia .cell-code}\nsum(f(cs[i]) * deltas[i] for i in 1:length(deltas))\n```\n\n::: {.cell-output .cell-output-display execution_count=14}\n```\n0.24000000000000002\n```\n:::\n:::\n\n\nOur answer is not so close to the value of $1/3$, but what did we expect - we only used $n=5$ intervals. Trying again with $50,000$ gives us:\n\n::: {.cell hold='true' execution_count=14}\n``` {.julia .cell-code}\nn = 50_000\nxs = a:(b-a)/n:b\ndeltas = diff(xs)\ncs = xs[1:end-1]\nsum(f(cs[i]) * deltas[i] for i in 1:length(deltas))\n```\n\n::: {.cell-output .cell-output-display execution_count=15}\n```\n0.3333233333999998\n```\n:::\n:::\n\n\nThis value is about $10^{-5}$ off from the actual answer of $1/3$.\n\n\nWe should expect that larger values of $n$ will produce better approximate values, as long as numeric issues don't get involved.\n\n\nBefore continuing, we define a function to compute the Riemann sum for us with an extra argument to specifying one of four methods for computing $c_i$:\n\n``` {.julia .cell-code}\nfunction riemann(f::Function, a::Real, b::Real, n::Int; method=\"right\")\n if method == \"right\"\n meth = f -> (lr -> begin l,r = lr; f(r) * (r-l) end)\n elseif method == \"left\"\n meth = f -> (lr -> begin l,r = lr; f(l) * (r-l) end)\n elseif method == \"trapezoid\"\n meth = f -> (lr -> begin l,r = lr; (1/2) * (f(l) + f(r)) * (r-l) end)\n elseif method == \"simpsons\"\n meth = f -> (lr -> begin l,r=lr; (1/6) * (f(l) + 4*(f((l+r)/2)) + f(r)) * (r-l) end)\n end\n\n xs = range(a, b, n+1)\n pairs = zip(xs[begin:end-1], xs[begin+1:end]) # (x₀,x₁), …, (xₙ₋₁,xₙ)\n sum(meth(f), pairs)\n\nend\n```\n\n\n(This function is defined in `CalculusWithJulia` and need not be copied over if that package is loaded.)\n\n\nWith this, we can easily find an approximate answer. We wrote the function to use the familiar template `action(function, arguments...)`, so we pass in a function and arguments to describe the problem (`a`, `b`, and `n` and, optionally, the `method`):\n\n::: {.cell execution_count=16}\n``` {.julia .cell-code}\n𝒇(x) = exp(x)\nriemann(𝒇, 0, 5, 10) # S_10\n```\n\n::: {.cell-output .cell-output-display execution_count=16}\n```\n187.324835773627\n```\n:::\n:::\n\n\nOr with more intervals in the partition\n\n::: {.cell execution_count=17}\n``` {.julia .cell-code}\nriemann(𝒇, 0, 5, 50_000)\n```\n\n::: {.cell-output .cell-output-display execution_count=17}\n```\n147.42052988337647\n```\n:::\n:::\n\n\n(The answer is $e^5 - e^0 = 147.4131591025766\\dots$, which shows that even $50,000$ partitions is not enough to guarantee many digits of accuracy.)\n\n\n## \"Negative\" area\n\n\nSo far, we have had the assumption that $f(x) \\geq 0$, as that allows us to define the concept of area. We can define the signed area between $f(x)$ and the $x$ axis through the definite integral:\n\n\n\n$$\nA = \\int_a^b f(x) dx.\n$$\n\n\nThe right hand side is defined whenever the Riemann limit exists and in that case we call $f(x)$ Riemann integrable. (The definition does not suppose $f$ is non-negative.)\n\n\nSuppose $f(a) = f(b) = 0$ for $a < b$ and for all $a < x < b$ we have $f(x) < 0$. Then we can see easily from the geometry (or from the Riemann sum approximation) that\n\n\n\n$$\n\\int_a^b f(x) dx = - \\int_a^b \\lvert f(x) \\rvert dx.\n$$\n\n\nIf we think of the area below the $x$ axis as \"signed\" area carrying a minus sign, then the total area can be seen again as a sum, only this time some of the summands may be negative.\n\n\n##### Example\n\n\nConsider a function $g(x)$ defined through its piecewise linear graph:\n\n::: {.cell execution_count=18}\n\n::: {.cell-output .cell-output-display execution_count=18}\n![](area_files/figure-html/cell-19-output-1.svg){}\n:::\n:::\n\n\n * Compute $\\int_{-3}^{-1} g(x) dx$. The area comprised of a square of area $1$ and a triangle with area $1/2$, so should be $3/2$.\n * Compute $\\int_{-3}^{0} g(x) dx$. In addition to the above, there is a triangle with area $1/2$, but since the function is negative, this area is added in as $-1/2$. In total then we have $1 + 1/2 - 1/2 = 1$ for the answer.\n * Compute $\\int_{-3}^{1} g(x) dx$:\n\n\nWe could add the signed area over $[0,1]$ to the above, but instead see a square of area $1$, a triangle with area $1/2$ and a triangle with signed area $-1$. The total is then $1/2$.\n\n\n * Compute $\\int_{-3}^{3} g(x) dx$:\n\n\nWe could add the area, but let's use a symmetry trick. This is clearly twice our second answer, or $2$. (This is because $g(x)$ is an even function, as we can tell from the graph.)\n\n\n##### Example\n\n\nSuppose $f(x)$ is an odd function, then $f(x) = - f(-x)$ for any $x$. So the signed area between $[-a,0]$ is related to the signed area between $[0,a]$ but of different sign. This gives $\\int_{-a}^a f(x) dx = 0$ for odd functions.\n\n\nAn immediate consequence would be $\\int_{-\\pi}^\\pi \\sin(x) = 0$, as would $\\int_{-a}^a x^k dx$ for any *odd* integer $k > 0$.\n\n\n##### Example\n\n\nNumerically estimate the definite integral $\\int_0^e x\\log(x) dx$. (We redefine the function to be $0$ at $0$, so it is continuous.)\n\n\nWe have to be a bit careful with the Riemann sum, as the left Riemann sum will have an issue at $0=x_0$ (`0*log(0)` returns `NaN` which will poison any subsequent arithmetic operations, so the value returned will be `NaN` and not an approximate answer). We could define our function with a check:\n\n::: {.cell execution_count=19}\n``` {.julia .cell-code}\n𝒉(x) = x > 0 ? x * log(x) : 0.0\n```\n\n::: {.cell-output .cell-output-display execution_count=19}\n```\n𝒉 (generic function with 1 method)\n```\n:::\n:::\n\n\nThis is actually inefficient, as the check for the size of `x` will slow things down a bit. Since we will call this function 50,000 times, we would like to avoid this, if we can. In this case just using the right sum will work:\n\n::: {.cell execution_count=20}\n``` {.julia .cell-code}\nh(x) = x * log(x)\nriemann(h, 0, 2, 50_000, method=\"right\")\n```\n\n::: {.cell-output .cell-output-display execution_count=20}\n```\n0.38632208884775826\n```\n:::\n:::\n\n\n(The default is `\"right\"`, so no method specified would also work.)\n\n\n##### Example\n\n\nLet $j(x) = \\sqrt{1 - x^2}$. The area under the curve between $-1$ and $1$ is $\\pi/2$. Using a Riemann sum with 4 equal subintervals and the midpoint, estimate $\\pi$. How close are you?\n\n\nThe partition is $-1 < -1/2 < 0 < 1/2 < 1$. The midpoints are $-3/4, -1/4, 1/4, 3/4$. We thus have that $\\pi/2$ is approximately:\n\n::: {.cell hold='true' execution_count=21}\n``` {.julia .cell-code}\nxs = range(-1, 1, length=5)\ndeltas = diff(xs)\ncs = [-3/4, -1/4, 1/4, 3/4]\nj(x) = sqrt(1 - x^2)\na = sum(j(c)*delta for (c,delta) in zip(cs, deltas))\na, pi/2 # π ≈ 2a\n```\n\n::: {.cell-output .cell-output-display execution_count=21}\n```\n(1.629683664318002, 1.5707963267948966)\n```\n:::\n:::\n\n\n(For variety, we used an alternate way to sum over two vectors.)\n\n\nSo $\\pi$ is about `2a`.\n\n\n##### Example\n\n\nWe have the well-known triangle [inequality](http://en.wikipedia.org/wiki/Triangle_inequality) which says for an individual sum: $\\lvert a + b \\rvert \\leq \\lvert a \\rvert +\\lvert b \\rvert$. Applying this recursively to a partition with $a < b$ gives:\n\n\n\n$$\n\\begin{align*}\n\\lvert f(c_1)(x_1-x_0) + f(c_2)(x_2-x_1) + \\cdots + f(c_n) (x_n-x_1) \\rvert\n& \\leq\n\\lvert f(c_1)(x_1-x_0) \\rvert + \\lvert f(c_2)(x_2-x_1)\\rvert + \\cdots +\\lvert f(c_n) (x_n-x_1) \\rvert \\\\\n&= \\lvert f(c_1)\\rvert (x_1-x_0) + \\lvert f(c_2)\\rvert (x_2-x_1)+ \\cdots +\\lvert f(c_n) \\rvert(x_n-x_1).\n\\end{align*}\n$$\n\n\nThis suggests that the following inequality holds for integrals:\n\n\n> $\\lvert \\int_a^b f(x) dx \\rvert \\leq \\int_a^b \\lvert f(x) \\rvert dx$.\n\n\n\nThis can be used to give bounds on the size of an integral. For example, suppose you know that $f(x)$ is continuous on $[a,b]$ and takes its maximum value of $M$ and minimum value of $m$. Letting $K$ be the larger of $\\lvert M\\rvert$ and $\\lvert m \\rvert$, gives this bound when $a < b$:\n\n\n\n$$\n\\lvert\\int_a^b f(x) dx \\rvert \\leq \\int_a^b \\lvert f(x) \\rvert dx \\leq \\int_a^b K dx = K(b-a).\n$$\n\n\nWhile such bounds are disappointing, often, when looking for specific values, they are very useful when establishing general truths, such as is done with proofs.\n\n\n## Error estimate\n\n\nThe Riemann sum above is actually extremely inefficient. To see how much, we can derive an estimate for the error in approximating the value using an arithmetic progression as the partition. Let's assume that our function $f(x)$ is increasing, so that the right sum gives an upper estimate and the left sum a lower estimate, so the error in the estimate will be between these two values:\n\n\n\n$$\n\\begin{align*}\n\\text{error} &\\leq\n\\left[\nf(x_1) \\cdot (x_{1} - x_0) + f(x_2) \\cdot (x_{2} - x_1) + \\cdots + f(x_{n-1})(x_{n-1} - x_n) + f(x_n) \\cdot (x_n - x_{n-1})\\right]\\\\\n&-\n\\left[f(x_0) \\cdot (x_{1} - x_0) + f(x_1) \\cdot (x_{2} - x_1) + \\cdots + f(x_{n-1})(x_{n-1} - x_n)\\right] \\\\\n&= \\frac{b-a}{n} \\cdot (\\left[f(x_1) + f(x_2) + \\cdots f(x_n)\\right] - \\left[f(x_0) + \\cdots f(x_{n-1})\\right]) \\\\\n&= \\frac{b-a}{n} \\cdot (f(b) - f(a)).\n\\end{align*}\n$$\n\n\nWe see the error goes to $0$ at a rate of $1/n$ with the constant depending on $b-a$ and the function $f$. In general, a similar bound holds when $f$ is not monotonic.\n\n\nThere are other ways to approximate the integral that use fewer points in the partition. [Simpson's](http://tinyurl.com/7b9pmu) rule is one, where instead of approximating the area with rectangles that go through some $c_i$ in $[x_{i-1}, x_i]$ instead the function is approximated by the quadratic polynomial going through $x_{i-1}$, $(x_i + x_{i-1})/2$, and $x_i$ and the exact area under that polynomial is used in the approximation. The explicit formula is:\n\n\n\n$$\nA \\approx \\frac{b-a}{3n} (f(x_0) + 4 f(x_1) + 2f(x_2) + 4f(x_3) + \\cdots + 2f(x_{n-2}) + 4f(x_{n-1}) + f(x_n)).\n$$\n\n\nThe error in this approximation can be shown to be\n\n\n\n$$\n\\text{error} \\leq \\frac{(b-a)^5}{180n^4} \\text{max}_{\\xi \\text{ in } [a,b]} \\lvert f^{(4)}(\\xi) \\rvert.\n$$\n\n\nThat is, the error is like $1/n^4$ with constants depending on the length of the interval, $(b-a)^5$, and the maximum value of the fourth derivative over $[a,b]$. This is significant, the error in $10$ steps of Simpson's rule is on the scale of the error of $10,000$ steps of the Riemann sum for well-behaved functions.\n\n\n:::{.callout-note}\n## Note\nThe Wikipedia article mentions that Kepler used a similar formula $100$ years prior to Simpson, or about $200$ years before Riemann published his work. Again, the value in Riemann's work is not the computation of the answer, but the framework it provides in determining if a function is Riemann integrable or not.\n\n:::\n\n## Gauss quadrature\n\n\nThe formula for Simpson's rule was the *composite* formula. If just a single rectangle is approximated over $[a,b]$ by a parabola interpolating the points $x_1=a$, $x_2=(a+b)/2$, and $x_3=b$, the formula is:\n\n\n\n$$\n\\frac{b-a}{6}(f(x_1) + 4f(x_2) + f(x_3)).\n$$\n\n\nThis formula will actually be exact for any 3rd degree polynomial. In fact an entire family of similar approximations using $n$ points can be made exact for any polynomial of degree $n-1$ or lower. But with non-evenly spaced points, even better results can be found.\n\n\nThe formulas for an approximation to the integral $\\int_{-1}^1 f(x) dx$ discussed so far can be written as:\n\n\n\n$$\n\\begin{align*}\nS &= f(x_1) \\Delta_1 + f(x_2) \\Delta_2 + \\cdots + f(x_n) \\Delta_n\\\\\n &= w_1 f(x_1) + w_2 f(x_2) + \\cdots + w_n f(x_n).\n\\end{align*}\n$$\n\n\nThe $w$s are \"weights\" and the $x$s are nodes. A [Gaussian](http://en.wikipedia.org/wiki/Gaussian_quadrature) *quadrature rule* is a set of weights and nodes for $i=1, \\dots n$ for which the sum is *exact* for any $f$ which is a polynomial of degree $2n-1$ or less. Such choices then also approximate well the integrals of functions which are not polynomials of degree $2n-1$, provided $f$ can be well approximated by a polynomial over $[-1,1]$. (Which is the case for the \"nice\" functions we encounter.) Some examples are given in the questions.\n\n\n### The quadgk function\n\n\nIn `Julia` a modification of the Gauss quadrature rule is implemented in the `quadgk` function (from the `QuadGK` package) to give numeric approximations to integrals. The `quadgk` function also has the familiar interface `action(function, arguments...)`. Unlike our `riemann` function, there is no `n` specified, as the number of steps is *adaptively* determined. (There is more partitioning occurring where the function is changing rapidly.) Instead, the algorithm outputs an estimate on the possible error along with the answer. Instead of $n$, some trickier problems require a specification of an error threshold.\n\n\nTo use the function, we have:\n\n::: {.cell hold='true' execution_count=22}\n``` {.julia .cell-code}\nf(x) = x * log(x)\nquadgk(f, 0, 2)\n```\n\n::: {.cell-output .cell-output-display execution_count=22}\n```\n(0.38629436103070175, 4.856575104393709e-9)\n```\n:::\n:::\n\n\nAs mentioned, there are two values returned: an approximate answer, and an error estimate. In this example we see that the value of $0.3862943610307017$ is accurate to within $10^{-9}$. (The actual answer is $-1 + 2\\cdot \\log(2)$ and the error is only $10^{-11}$. The reported error is an upper bound, and may be conservative, as with this problem.) Our previous answer using $50,000$ right-Riemann sums was $0.38632208884775737$ and is only accurate to $10^{-5}$. By contrast, this method uses just $256$ function evaluations in the above problem.\n\n\nThe method should be exact for polynomial functions:\n\n::: {.cell hold='true' execution_count=23}\n``` {.julia .cell-code}\nf(x) = x^5 - x + 1\nquadgk(f, -2, 2)\n```\n\n::: {.cell-output .cell-output-display execution_count=23}\n```\n(3.9999999999999973, 8.881784197001252e-16)\n```\n:::\n:::\n\n\nThe error term is $0$, answer is $4$ up to the last unit of precision (1 ulp), so any error is only in floating point approximations.\n\n\nFor the numeric computation of definite integrals, the `quadgk` function should be used over the Riemann sums or even Simpson's rule.\n\n\nHere are some sample integrals computed with `quadgk`:\n\n\n\n$$\n\\int_0^\\pi \\sin(x) dx\n$$\n\n::: {.cell execution_count=24}\n``` {.julia .cell-code}\nquadgk(sin, 0, pi)\n```\n\n::: {.cell-output .cell-output-display execution_count=24}\n```\n(2.0, 1.7905676941154525e-12)\n```\n:::\n:::\n\n\n(Again, the actual answer is off only in the last digit, the error estimate is an upper bound.)\n\n\n\n$$\n\\int_0^2 x^x dx\n$$\n\n::: {.cell execution_count=25}\n``` {.julia .cell-code}\nquadgk(x -> x^x, 0, 2)\n```\n\n::: {.cell-output .cell-output-display execution_count=25}\n```\n(2.8338767448900546, 1.9481752001546115e-8)\n```\n:::\n:::\n\n\n$$\n\\int_0^5 e^x dx\n$$\n\n::: {.cell execution_count=26}\n``` {.julia .cell-code}\nquadgk(exp, 0, 5)\n```\n\n::: {.cell-output .cell-output-display execution_count=26}\n```\n(147.41315910257657, 2.6594506152832764e-8)\n```\n:::\n:::\n\n\nWhen composing the answer with other functions it may be desirable to drop the error in the answer. Two styles can be used for this. The first is to just name the two returned values:\n\n::: {.cell hold='true' execution_count=27}\n``` {.julia .cell-code}\nA, err = quadgk(cos, 0, pi/4)\nA\n```\n\n::: {.cell-output .cell-output-display execution_count=27}\n```\n0.7071067811865475\n```\n:::\n:::\n\n\nThe second is to ask for just the first component of the returned value:\n\n::: {.cell hold='true' execution_count=28}\n``` {.julia .cell-code}\nA = quadgk(tan, 0, pi/4)[1] # or first(quadgk(tan, 0, pi/4))\n```\n\n::: {.cell-output .cell-output-display execution_count=28}\n```\n0.3465735902799726\n```\n:::\n:::\n\n\n---\n\n\nTo visualize the choice of nodes by the algorithm, we have for $f(x)=\\sin(x)$ over $[0,\\pi]$ relatively few nodes used to get a high-precision estimate:\n\n\n\n::: {.cell hold='true' execution_count=30}\n\n::: {.cell-output .cell-output-display execution_count=30}\n![](area_files/figure-html/cell-31-output-1.svg){}\n:::\n:::\n\n\nFor a more oscillatory function, more nodes are chosen:\n\n::: {.cell hold='true' execution_count=31}\n\n::: {.cell-output .cell-output-display execution_count=31}\n![](area_files/figure-html/cell-32-output-1.svg){}\n:::\n:::\n\n\n##### Example\n\n\nIn probability theory, a *univariate density* is a function, $f(x)$ such that $f(x) \\geq 0$ and $\\int_a^b f(x) dx = 1$, where $a$ and $b$ are the range of the distribution. The [Von Mises](http://en.wikipedia.org/wiki/Von_Mises_distribution) distribution, takes the form\n\n\n\n$$\nk(x) = C \\cdot \\exp(\\cos(x)), \\quad -\\pi \\leq x \\leq \\pi\n$$\n\n\nCompute $C$ (numerically).\n\n\nThe fact that $1 = \\int_{-\\pi}^\\pi C \\cdot \\exp(\\cos(x)) dx = C \\int_{-\\pi}^\\pi \\exp(\\cos(x)) dx$ implies that $C$ is the reciprocal of\n\n::: {.cell execution_count=32}\n``` {.julia .cell-code}\nk(x) = exp(cos(x))\nA,err = quadgk(k, -pi, pi)\n```\n\n::: {.cell-output .cell-output-display execution_count=32}\n```\n(7.954926521012919, 3.9023298370466364e-8)\n```\n:::\n:::\n\n\nSo\n\n::: {.cell execution_count=33}\n``` {.julia .cell-code}\nC = 1/A\nk₁(x) = C * exp(cos(x))\n```\n\n::: {.cell-output .cell-output-display execution_count=33}\n```\nk₁ (generic function with 1 method)\n```\n:::\n:::\n\n\nThe *cumulative distribution function* for $k(x)$ is $K(x) = \\int_{-\\pi}^x k(u) du$, $-\\pi \\leq x \\leq \\pi$. We just showed that $K(\\pi) = 1$ and it is trivial that $K(-\\pi) = 0$. The quantiles of the distribution are the values $q_1$, $q_2$, and $q_3$ for which $K(q_i) = i/4$. Can we find these?\n\n\nFirst we define a function, that computes $K(x)$:\n\n::: {.cell execution_count=34}\n``` {.julia .cell-code}\nK(x) = quadgk(k₁, -pi, x)[1]\n```\n\n::: {.cell-output .cell-output-display execution_count=34}\n```\nK (generic function with 1 method)\n```\n:::\n:::\n\n\n(The trailing `[1]` is so only the answer - and not the error - is returned.)\n\n\nThe question asks us to solve $K(x) = 0.25$, $K(x) = 0.5$ and $K(x) = 0.75$. The `Roots` package can be used for such work, in particular `find_zero`. We will use a bracketing method, as clearly $K(x)$ is increasing, as $k(u)$ is positive, so we can just bracket our answer with $-\\pi$ and $\\pi$. (We solve $K(x) - p = 0$, so $K(\\pi) - p > 0$ and $K(-\\pi)-p < 0$.). We could do this with `[find_zero(x -> K(x) - p, (-pi, pi)) for p in [0.25, 0.5, 0.75]]`, but that is a bit less performant than using the `solve` interface for this task:\n\n::: {.cell hold='true' execution_count=35}\n``` {.julia .cell-code}\nZ = ZeroProblem((x,p) -> K(x) - p, (-pi, pi))\nsolve.(Z, (1/4, 1/2, 3/4))\n```\n\n::: {.cell-output .cell-output-display execution_count=35}\n```\n(-0.8097673745015915, 0.0, 0.809767374501629)\n```\n:::\n:::\n\n\nThe middle one is clearly $0$. This distribution is symmetric about $0$, so half the area is to the right of $0$ and half to the left, so clearly when $p=0.5$, $x$ is $0$. The other two show that the area to the left of $-0.809767$ is equal to the area to the right of $0.809767$ and equal to $0.25$.\n\n\n## Questions\n\n\n###### Question\n\n\nUsing geometry, compute the definite integral:\n\n\n\n$$\n\\int_{-5}^5 \\sqrt{5^2 - x^2} dx.\n$$\n\n::: {.cell hold='true' execution_count=36}\n\n::: {.cell-output .cell-output-display execution_count=36}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nUsing geometry, compute the definite integral:\n\n\n\n$$\n\\int_{-2}^2 (2 - \\lvert x\\rvert) dx\n$$\n\n::: {.cell hold='true' execution_count=37}\n\n::: {.cell-output .cell-output-display execution_count=37}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nUsing geometry, compute the definite integral:\n\n\n\n$$\n\\int_0^3 3 dx + \\int_3^9 (3 + 3(x-3)) dx\n$$\n\n::: {.cell hold='true' execution_count=38}\n\n::: {.cell-output .cell-output-display execution_count=38}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nUsing geometry, compute the definite integral:\n\n\n\n$$\n\\int_0^5 \\lfloor x \\rfloor dx\n$$\n\n\n(The notation $\\lfloor x \\rfloor$ is the integer such that $\\lfloor x \\rfloor \\leq x < \\lfloor x \\rfloor + 1$.)\n\n::: {.cell hold='true' execution_count=39}\n\n::: {.cell-output .cell-output-display execution_count=39}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nUsing geometry, compute the definite integral between $-3$ and $3$ of this graph comprised of lines and circular arcs:\n\n::: {.cell hold='true' execution_count=40}\n\n::: {.cell-output .cell-output-display execution_count=40}\n![](area_files/figure-html/cell-41-output-1.svg){}\n:::\n:::\n\n\nThe value is:\n\n::: {.cell hold='true' execution_count=41}\n\n::: {.cell-output .cell-output-display execution_count=41}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nFor the function $f(x) = \\sin(\\pi x)$, estimate the integral for $-1$ to $1$ using a left-Riemann sum with the partition $-1 < -1/2 < 0 < 1/2 < 1$.\n\n::: {.cell hold='true' execution_count=42}\n\n::: {.cell-output .cell-output-display execution_count=42}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nWithout doing any *real* work, find this integral:\n\n\n\n$$\n\\int_{-\\pi/4}^{\\pi/4} \\tan(x) dx.\n$$\n\n::: {.cell hold='true' execution_count=43}\n\n::: {.cell-output .cell-output-display execution_count=43}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nWithout doing any *real* work, find this integral:\n\n\n\n$$\n\\int_3^5 (1 - \\lvert x-4 \\rvert) dx\n$$\n\n::: {.cell hold='true' execution_count=44}\n\n::: {.cell-output .cell-output-display execution_count=44}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nSuppose you know that for the integrable function $\\int_a^b f(u)du =1$ and $\\int_a^c f(u)du = p$. If $a < c < b$ what is $\\int_c^b f(u)du$?\n\n::: {.cell hold='true' execution_count=45}\n\n::: {.cell-output .cell-output-display execution_count=45}\n```{=html}\n
\n
\n
\n
\n
\n
\n \n
\n
\n \n
\n
\n \n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nWhat is $\\int_0^2 x^4 dx$? Use the rule for integrating $x^n$.\n\n::: {.cell hold='true' execution_count=46}\n\n::: {.cell-output .cell-output-display execution_count=46}\n```{=html}\n
\n
\n
\n
\n
\n
\n \n
\n
\n \n
\n
\n \n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nSolve for a value of $x$ for which:\n\n\n\n$$\n\\int_1^x \\frac{1}{u}du = 1.\n$$\n\n::: {.cell hold='true' execution_count=47}\n\n::: {.cell-output .cell-output-display execution_count=47}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nSolve for a value of $n$ for which\n\n\n\n$$\n\\int_0^1 x^n dx = \\frac{1}{12}.\n$$\n\n::: {.cell hold='true' execution_count=48}\n\n::: {.cell-output .cell-output-display execution_count=48}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nSuppose $f(x) > 0$ and $a < c < b$. Define $F(x) = \\int_a^x f(u) du$. What can be said about $F(b)$ and $F(c)$?\n\n::: {.cell hold='true' execution_count=49}\n\n::: {.cell-output .cell-output-display execution_count=49}\n```{=html}\n
\n
\n
\n
\n
\n
\n \n
\n
\n \n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nFor the right Riemann sum approximating $\\int_0^{10} e^x dx$ with $n=100$ subintervals, what would be a good estimate for the error?\n\n::: {.cell hold='true' execution_count=50}\n\n::: {.cell-output .cell-output-display execution_count=50}\n```{=html}\n
\n
\n
\n
\n
\n
\n \n
\n
\n \n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nUse `quadgk` to find the following definite integral:\n\n\n\n$$\n\\int_1^4 x^x dx .\n$$\n\n::: {.cell hold='true' execution_count=51}\n\n::: {.cell-output .cell-output-display execution_count=51}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nUse `quadgk` to find the following definite integral:\n\n\n\n$$\n\\int_0^3 e^{-x^2} dx .\n$$\n\n::: {.cell hold='true' execution_count=52}\n\n::: {.cell-output .cell-output-display execution_count=52}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nUse `quadgk` to find the following definite integral:\n\n\n\n$$\n\\int_0^{9/10} \\tan(u \\frac{\\pi}{2}) du. .\n$$\n\n::: {.cell hold='true' execution_count=53}\n\n::: {.cell-output .cell-output-display execution_count=53}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nUse `quadgk` to find the following definite integral:\n\n\n\n$$\n\\int_{-1/2}^{1/2} \\frac{1}{\\sqrt{1 - x^2}} dx\n$$\n\n::: {.cell hold='true' execution_count=54}\n\n::: {.cell-output .cell-output-display execution_count=54}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\n\n```{=html}\n
\n```\n\n```{ojs}\n//| echo: false\n//| output: false\nJXG = require(\"jsxgraph\");\n\nb = JXG.JSXGraph.initBoard('jsxgraph', {\n boundingbox: [-0.5,0.3,1.5,-1/4], axis:true\n});\n\ng = function(x) { return x*x*x*x + 10*x*x - 60* x + 100}\nf = function(x) {return 1/Math.sqrt(g(x))};\n\ntype = \"right\";\nl = 0;\nr = 1;\nrsum = function() {\n return JXG.Math.Numerics.riemannsum(f,n.Value(), type, l, r);\n};\nn = b.create('slider', [[0.1, -0.05],[0.75,-0.05], [2,1,50]],{name:'n',snapWidth:1});\n\ngraph = b.create('functiongraph', [f, l, r]);\nos = b.create('riemannsum',\n [f,\n function(){ return n.Value();},\n type, l, r\n ],\n {fillColor:'#ffff00', fillOpacity:0.3});\n\nb.create('text', [0.1,0.25, function(){\n return 'Riemann sum='+(rsum().toFixed(4));\n}]);\n```\n\n\nThe interactive graphic shows the area of a right-Riemann sum for different partitions. The function is\n\n\n\n$$\nf(x) = \\frac{1}{\\sqrt{ x^4 + 10x^2 - 60x + 100}}\n$$\n\n\nWhen $n=5$ what is the area of the Riemann sum?\n\n::: {.cell hold='true' execution_count=56}\n\n::: {.cell-output .cell-output-display execution_count=56}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\nWhen $n=50$ what is the area of the Riemann sum?\n\n::: {.cell hold='true' execution_count=57}\n\n::: {.cell-output .cell-output-display execution_count=57}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\nUsing `quadgk` what is the area under the curve?\n\n::: {.cell hold='true' execution_count=58}\n\n::: {.cell-output .cell-output-display execution_count=58}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nGauss nodes for approximating the integral $\\int_{-1}^1 f(x) dx$ for $n=4$ are:\n\n::: {.cell execution_count=59}\n``` {.julia .cell-code}\nns = [-0.861136, -0.339981, 0.339981, 0.861136]\n```\n\n::: {.cell-output .cell-output-display execution_count=59}\n```\n4-element Vector{Float64}:\n -0.861136\n -0.339981\n 0.339981\n 0.861136\n```\n:::\n:::\n\n\nThe corresponding weights are\n\n::: {.cell execution_count=60}\n``` {.julia .cell-code}\nwts = [0.347855, 0.652145, 0.652145, 0.347855]\n```\n\n::: {.cell-output .cell-output-display execution_count=60}\n```\n4-element Vector{Float64}:\n 0.347855\n 0.652145\n 0.652145\n 0.347855\n```\n:::\n:::\n\n\nUse these to estimate the integral $\\int_{-1}^1 \\cos(\\pi/2 \\cdot x)dx$ with $w_1f(x_1) + w_2 f(x_2) + w_3 f(x_3) + w_4 f(x_4)$.\n\n::: {.cell hold='true' execution_count=61}\n\n::: {.cell-output .cell-output-display execution_count=61}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\nThe actual answer is $4/\\pi$. How far off is the approximation based on 4 points?\n\n::: {.cell hold='true' execution_count=62}\n\n::: {.cell-output .cell-output-display execution_count=62}\n```{=html}\n
\n
\n
\n
\n
\n
\n \n
\n
\n \n
\n
\n \n
\n
\n \n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n###### Question\n\n\nUsing the Gauss nodes and weights from the previous question, estimate the integral of $f(x) = e^x$ over $[-1, 1]$. The value is:\n\n::: {.cell hold='true' execution_count=63}\n\n::: {.cell-output .cell-output-display execution_count=63}\n```{=html}\n
\n
\n
\n
\n
\n
\n
\n \n
\n\n \n
\n
\n
\n
\n
\n
\n\n\n```\n:::\n:::\n\n\n", "supporting": [ "area_files" ], "filters": [], "includes": { "include-in-header": [ "\n\n\n" ] } } }