syms xo; syms x;
syms yo; syms y;
f_zerowe = [2*yo+4*xo*yo; 1+2*xo+2*xo^2-2*yo^2];
f_pierwsze = [diff(1+2*xo+2*xo^2-2*yo^2+j*(2*yo+4*xo*yo),xo), diff(1+2*xo+2*xo^2-2*yo^2+j*(2*yo+4*xo*yo),yo)]*[x;y];
f_drugie = 0.5*[x,y]*hessian(1+2*xo+2*xo^2-2*yo^2+j*(2*yo+4*xo*yo), [xo,yo])*[x;y];
cale = f_zerowe + f_pierwsze + f_drugie

% wstawiamy zmienne
pom = subs(cale,xo,0);
pom = subs(pom, x, 0);
pom = subs(pom, yo, 0);
pom = subs(pom, y, 0);
pom'