clear;

syms x f

f = x^2 - 3;
x0 = 1;
tol = 10^(-5);
n = 10;

xs = newt(f, x0, tol, n)