# Numerical Nonlinear Global Optimization Examples

## Finding Multiple Optima, Method 1

Here is one way to get multiple minima: call NMinimize multiple times with different random seeds, which will cause different optimization paths to be taken.

This defines a function with a whole ring of minima.
 In[31]:=
 Out[32]=
This makes a table of solutions by using different random seeds.
 In[3]:=
 Out[3]=
This turns the solutions into points and plots them on the contour plot of the function.
 In[4]:=
 Out[7]=

## Finding Multiple Optima, Method 2

Here is another way to get multiple minima: write the objective function in such a way as to make a list of every point that is visited, then select the points that have objective function values close to the final solution.

Define a function with a whole ring of minima.
 In[50]:=
 Out[52]=
Reap returns the solution from NMinimize and the points sown by EvaluationMonitor.
 In[21]:=
 Out[22]=
Find all the visited points that are close in objective function value to the final solution. Then color them white, color the final solution black, and show them on the contour plot of the function.
 In[23]:=
 Out[27]=

## Finding a Nonlinear Fit of Data

This defines a model based on five random parameters.
 In[4]:=
 Out[6]=
Create a function from the model and parameters, and generate sample points over the interval .
 In[7]:=
This plots the points and the solution from FindFit. The solution gets trapped by a local minimum due to the trigonometric functions.
 In[9]:=
 Out[11]=
This generates a sum of squares from the data, and uses NMinimize to find the minimum.
 In[12]:=
 Out[13]=
This plots the points and the solution from NMinimize.
 In[14]:=
 Out[15]=

## Solve Example

Solve cannot work with this system of equations because they are highly nonalgebraic.
 In[44]:=
 Out[44]=
Give NMinimize a constant objective function, and the equations to be solved as constraints. It finds the solution.
 In[43]:=
 Out[43]=

## Queens on a Chessboard

 In[46]:=
attackQ[pos1,pos2] is True if and only if pos1 is attacking pos2.
 In[47]:=
countAttacks[vec] converts the vector of real numbers into a permutation of the queens and counts the number of attacks.
 In[48]:=
Given a permutation, this shows the arrangement.
 In[49]:=
Use DifferentialEvolution to fit all the queens on the chessboard so that no queen is attacking another queen. Postprocessing is turned off because it is unlikely to help, given the discrete nature of the problem.
 In[50]:=
 Out[52]=
This shows the solution.
 In[53]:=
 Out[53]=