[x,fval] = fminsearch(@(x)100*(x(2)-x(1)^2)^2+(1-x(1))^2,[-1.2, 1]);
save('C:\STIMATH\Output\MyOut18.dat','x','fval','-ascii')
