Вычислить локальные и глобальные минимумы многопараметрической функции в Python

Question

Я хотел бы вычислить все минимумы (локальные и глобальные) следующей функции с начальной точкой (2, 20):

I пробовал с помощью следующего кода, но я не уверен, правильно ли я рассчитал и возможно ли найти все локальные минимумы в Python.

import numpy as np
from scipy.optimize import fmin
import math

def f(x):
#     return x[0]**2/2 + 5*x[1]**2 + x[0]**3/3
    return (1/2)*np.power((x[0]), 2) + 5*np.power((x[1]), 2) + (1/3)*np.power((x[0]), 3)

fmin(f, np.array([2, 20]))

Вывод:

Optimization terminated successfully.
         Current function value: 0.000000
         Iterations: 58
         Function evaluations: 112
array([-1.7438898e-05,  2.0049784e-06])

Пожалуйста, помогите мне, спасибо.