Das Finden lokaler Extremwerte erfordert die Angabe von Suchbereichen. Deshalb ist es zweckmäßig zuerst die Funktion graphisch darzustellen, um daraus die Suchbereiche zu entnehmen.
Remove["Global`*"] f= Plot3D[f,{x,-10,10},{y,-10,10}]
Aus der Abbildung werden die ungefähren Lagen der Extremwerte entnommen.
xa:= (* Start- und Endwerte des Suchbereiches für x *) xe:= ya:= (* Start- und Endwerte des Suchbereiches für y *) ye:= FindMinimum[f,{x,{xa,xe}},{y,{ya,ye}}] (* Minimumsuche *) FindMinimum[-f,{x,{xa,xe}},{y,{ya,ye}}] (* Maximumsuche *) (* Achtung! Der erste Wert bei der Maximumsuche muß mit -1 multipliziert, um das Maximum zu erhalten *)