Riemann sums

This is one of the more complex examples. It draws Riemann partitions from point a to point b. While computing the partitions, their areas are summed up to show how the Riemann approximation compares to the true area under the given curve.

In this example, some extra markup is used outside of Mafs to provide some inputs into the Mafs visualization. This is common, and recommended! Movable points are not the only way to provide inputs to Mafs.

Partitions:
import * as React from "react" import { Mafs, useMovablePoint, Plot, Polygon, Text, CartesianCoordinates } from "mafs" import sumBy from "lodash/sumBy" import range from "lodash/range" interface Partition { polygon: [number, number][] area: number } function RiemannSum() { const maxNumPartitions = 200 // Inputs const [numPartitions, setNumPartitions] = React.useState(40) const lift = useMovablePoint([0, -1], { constrain: "vertical", }) const a = useMovablePoint([1, 0], { constrain: "horizontal", }) const b = useMovablePoint([11, 0], { constrain: "horizontal", }) // The function const wave = (x: number) => Math.sin(3 * x) + x ** 2 / 20 - 2 + lift.y + 2 const integral = (x: number) => (1 / 60) * (x ** 3 - 20 * Math.cos(3 * x)) + lift.y * x // Outputs const exactArea = integral(b.x) - integral(a.x) const dx = (b.x - a.x) / numPartitions const partitions: Partition[] = range( a.x, b.x - dx / 2, dx ).map((x) => { const yMid = wave(x + dx / 2) return { polygon: [ [x, 0], [x, yMid], [x + dx, yMid], [x + dx, 0], ], area: dx * yMid, } }) const areaApprox = sumBy(partitions, "area") return ( <> <Mafs viewBox={{ x: [-1, 12], y: [-3, 10] }} > <CartesianCoordinates subdivisions={2} /> <Plot.OfX y={wave} color="#358CF1" /> {partitions.map((partition, index) => ( <Polygon key={index} points={partition.polygon} fillOpacity={numPartitions / maxNumPartitions} color={ partition.area >= 0 ? "hsl(112, 100%, 47%)" : "hsl(0, 100%, 47%)" } /> ))} <Text attach="e" x={1.2} y={5.5} size={20}> Midpoint Riemann sum: </Text> <Text attach="e" x={1.2} y={4.5} size={30}> {areaApprox.toFixed(4)} </Text> <Text attach="e" x={1.2} y={3.5} size={20}> True area: </Text> <Text attach="e" x={1.2} y={2.5} size={30}> {exactArea.toFixed(4)} </Text> {lift.element} {a.element} {b.element} </Mafs> {/* These classnames are part of the Mafs docs website—they won't work for you. */} <div className="p-4 border-gray-700 border-t bg-black text-white"> Partitions:{" "} <input type="range" min={20} max={200} value={numPartitions} onChange={(event) => setNumPartitions(+event.target.value) } /> </div> </> ) }