bisection method, iteration does not stop although hand calculation were diffrent
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
while running this scrept it keeps loopin on the same iteration for some reason
while doing the bisection by hand error reached below zero after 5 iterations
end
% Setting x as symbolic variable
syms q w e r t y u i o p l k j h g f d s a z x c v b n m ;
% Input Section
y = input('Enter the equation equations: ');
a = input('Enter A value of the interval: ');
b = input('Enter B value of the interval: ');
e = input('Entre the minimum tolerable error: ');
% Finding Functional Value
fa = eval(subs(y,x,a));
fb = eval(subs(y,x,b));
% Implementing Bisection Method
if fa*fb > 0
disp('Given initial values do not bracket the root.');
else
c = (a+b)/2;
fc = eval(subs(y,x,c));
fprintf('\n\na\t\t\tb\t\t\tc\t\t\tf(c)\n');
while abs(fc)>e
fprintf('%f\t%f\t%f\t%f\n',a,b,c,fc);
if fa*fc< 0
b =c;
else
a =c;
end
c = (a+b)/2;
fc = eval(subs(y,x,c));
end
fprintf('\nRoot is: %f\n', c);
end
my inputs are
>> bisection2
Enter the equation equations: (sqrt(19.62.*x))*tanh((sqrt(x*19.62)/8)*2.5)
Enter A value of the interval: 0
Enter B value of the interval: 2
Entre the minimum tolerable error: 0.00005
the thing afterwards goes on forever unless i stop the script
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1217137/image.gif)
this the function i am trying to the find the root of
3 comentarios
Respuestas (1)
Steven Lord
el 2 de Dic. de 2022
Let's look at your function.
syms x
f = sqrt(19.62*x);
fun = f*tanh(2.5*f/8);
fplot(fun)
yline(0)
It looks like you have a zero crossing at x = 0 but you never check for that in your code. You only check that the value of the function at the midpoint of your interval is 0. Let's check numerically what we observed graphically.
value = subs(fun, x, 0)
By the way, stop using eval on symbolic expressions. If you need to convert the result of the subs call into a double precision value, call double on it.
d = double(value)
0 comentarios
Ver también
Categorías
Más información sobre Number Theory en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!