Add a loop index, e.g. ii. Also, you set your code as a function. This means you can call it, e.g., from the command window using
c = newton(1, 1e-3)
and thus it does not make sense to define x0 and delta in the function.
function c = newton(x0, delta)
format long e
% x0 = 5;
% delta = 10^-8;
c = x0;
fc = f(x0);
fprintf('Initial guess: c = %d, fc = %d\n',c,fc)
if abs(fc) <= delta
return;
end
ii = 0;
while abs(fc) > delta
ii = ii + 1;
fpc = fprime(c);
if (fpc==0)
error('fprime is 0')
end
c = c - fc/fpc;
fc = f(c);
fprintf('ii = %d: c = %10.7e, fc = %10.7e\n',ii,c,fc)
end
end
function fx = f(x)
fx = sin(x);
end
function fprimex = fprime(x)
fprimex = cos(x);
end
