% Funkcija ciju nulu trazimo
f = [sym('1/3*(2*x^2-1)-y'); 
     sym('2*y^3/(3*y-1)-x')];
 
% Pocetna iteracija, tacnost i maksimalni broj iteracija
x0 = [-1; 0.3]; 
epsilon = 1e-5; 
maxiter = 10;

% Poziv iterativnog algoritma
r2_newton(f, x0, epsilon, maxiter);