Tu napisałam program dla funkcji o 1 zmiennej.
import scipy
from scipy.optimize import fmin
def f(x): return -2 * x**2 + 4 * x
max_x = scipy.optimize.fmin(lambda x: -f(x), 0)
print max_x
A teraz próbuję napisać to samo dla funkcji -x^2 - y^2, tylko trudno jest mi zrozumieć jak mam to zrobić.