%program do aproksymowania szeregiem Teylora na dziedzinie dwuwymiarowej

%tworzenie zmeinnych symbolicznych
syms xo; syms x;
syms yo; syms y;
%implementacja poszczególnych funkcji
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
%wypisanie całego wzoru na zmiennych symbolicznych
pretty(cale)
% przypisanie wartości zmiennym x0 i y0
a = 0; %dla xo
b = 0; %dla yo
pom = subs(cale,xo,a);
pom = subs(pom, yo, b);
%wypisanie całego wektora na 2 zmiennych symbolicznych
pretty(pom)