Google translation:
Wrong output when starting the program
Code:
f = @(x) x - 5 - 1/3*cos(2*x+1) + (2*x)/(2+x^2);
g = @(x) x - f(x);
a = 0; b = 1; % interval on which we are looking for a solution
M1 = @(x) max(abs(4*del2(f(x)))); % estimate of the modulus of the second order derivative
eps = 1e-5; % specified accuracy
N_apr = @(x) ceil(log(M1(x)*(b-a)/eps)/log(2)); % a priori estimate of the number of iterations
x0 = (a+b)/2; % initial approximation
N_post = 0; % iteration counter
while true
x1 = g(x0);
N_post = N_post + 1;
if abs(x1-x0) < eps*(1-M1(x1))/M1(x1) % check the iteration termination condition
break;
end
x0 = x1;
end
disp(['Approximate solution: x = ', num2str(x1)]);
disp(['Number of iterations: N = ', num2str(N_post)]);
disp(['A priori estimate of the number of iterations: N_apr = ', num2str(N_apr(x1))]);
Conclusion:
Approximate solution: x = 4.4168
Number of iterations: N = 1
A priori estimate of the number of iterations: N_apr = -Inf
How to fix what's in N_apr = -inf?
