from scipy import optimize """ Minimizes function (x-a)^2+(y-b)^2 using powell's method """ def func((x,y),(a,b)) : return (x-a)**2 + (y-b)**2 print optimize.fmin_powell(func,(0,0),((5,6),))