function [x,y]=njutn(f1,f2,x0,y0,tol)
%f1 i f2 anonimne funkcije
f1=sym(f1);
f2=sym(f2);
J=[diff(f1,'x'),diff(f1,'y');diff(f2,'x'),diff(f2,'y')]
X=[x0;y0];
while 1
    inverz=inv(subs(J,{'x','y'},{x0,y0}));
    X1=X-inverz*[subs(f1,{'x','y'},{x0,y0});subs(f2,{'x','y'},{x0,y0})];
    X1=double(X1);
    if norm(X-X1,'inf')<tol
        x=X1(1);
        y=X1(2);
        break
    end
    X=X1;
    x0=X(1);
    y0=X(2);
end





%f1=f1(x,y)-simbolicka funkcija dve prom (npr f1=@(x,y) sin(x+y))
%f1(5,1) se dobija ovako: subs(f1,{'x','y'},{5,1})