4.10 Searching for Minima and Maxima of Functions

The minimise and maximise commands can be used to find the minima or maxima of algebraic expressions. In each case, a single algebraic expression should be supplied for optimisation, together with a comma-separated list of the variables with respect to which it should be optimised. In the following example, a minimum of the sinusoidal function $\cos (x)$ is sought:

pyxplot> set numerics real
pyxplot> x=0.1
pyxplot> minimise cos(x) via x
pyxplot> print x/pi
1

Note that this particular example doesn’t work when complex arithmetic is enabled, since $\cos (x)$ diverges to $-\infty $ at $x=\pi +\infty i$.

Various caveats apply both to the minimise and maximise commands, as well as to the solve command. All of these commands operate by searching numerically for optimal sets of input parameters to meet the criteria set by the user. As with all numerical algorithms, there is no guarantee that the locally optimum solutions returned are the globally optimum solutions. It is always advisable to double-check that the answers returned agree with common sense.

These commands can often find solutions to equations when these solutions are either very large or very small, but they usually work best when the solution they are looking for is roughly of order unity. PyXPlot does have mechanisms which attempt to correct cases where the supplied initial guess turns out to be many orders of magnitude different from the true solution, but it cannot be guaranteed not to wildly overshoot and produce unexpected results in such cases. To reiterate, it is always advisable to double-check that the answers returned agree with common sense.

Finding the maximum of a blackbody curve.

When a surface is heated to any given temperature $T$, it radiates thermally. The amount of electromagnetic radiation emitted at any particular frequency, per unit area of surface, per unit frequency of light, is given by the Planck Law:

  \[  B_\nu (\nu ,T)=\left(\frac{2h^3}{c^2}\right)\frac{\nu ^3}{\exp (h\nu /kT)-1}  \]    

The visible surface of the Sun has a temperature of approximately $5800\, \mathrm{K}$ and radiates in such a fashion. In this example, we use the solve, minimise and maximise commands to locate the frequency of light at which it emits the most energy per unit frequency interval. This task is simplified as PyXPlot has a system-defined mathematical function Bv(nu,T) which evaluates the expression given above.

Below, a plot is shown of the Planck Law for $T=5800\, \mathrm{K}$ to aid in visualising the solution to this problem:

\includegraphics[width=\textwidth ]{examples/eps/ex_eqnsolve}

To search for the maximum of this function using the maximise command, we must provide an initial guess to indicate that the answer sought should have units of Hz:

pyxplot> nu = 1e14*unit(Hz)
pyxplot> maximise Bv(nu,5800*unit(K)) via nu
pyxplot> print nu
340.9781 THz

This maximum could also be sought be searching for turning points in the function $B_\nu (\nu ,T)$, i.e. by solving the equation

  \[  \frac{\mathrm{d}B_\nu (\nu ,T)}{\mathrm{d}\nu }=0.  \]    

This can be done as follows:

pyxplot> nu = 1e14*unit(Hz)
pyxplot> solve diff_dnu(Bv(nu,5800*unit(K)),nu) = $\backslash $
.......>          unit(0*W/Hz**2/m**2) via nu
pyxplot> print nu
340.9781 THz

Finally, this maximum could also be found using PyXPlot’s built-in function Bvmax(T):
pyxplot> print Bvmax(5800*unit(K))
340.9781 THz