function y=fun(x);
y=x(1)^4+x(2)^4+2*(x(1)^2)*(x(2)^2)-4*x(1)+3;
%â îñíîâíîé ïðîãðàììå
%x0=[0;1];
%[x,feval]=fminsearch('fn',x0)