ex_fminsearch.m

d=[  0.51355  0.57819  0.60880  0.77224  0.80556  0.87892  0.91661
     5.1744   5.8805   6.2325   8.2532   8.6894   9.6982  10.2318];
x0=[2 2];
f = inline('sum(((1+x(1)*d(1,:)).^x(2) - d(2,:)).^2)', 'x', 'd');
[x fval, ok, output]= fminsearch(f,x0,[],d)