function iter=stef(func,x0,accur,num_iter) iter = []; iter(1) = x0; if(func(iter(1)) == 0) return; end for(i = 1:(num_iter)) iter(i+1) = iter(i) - ((func(iter(i)))^2)/(func(iter(i)+func(iter(i))) - func(iter(i))); if(abs(iter(i+1)-iter(i)) < 10^(-accur)) return; end if(func(iter(i+1)) == 0) return; end end end