Simple fixed-point iteration method
49 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
John Smith
el 22 de Sept. de 2019
Respondida: Tsega'ab
el 13 de Dic. de 2023
My task is to implement (simple) fixed-point interation.
So far, I've got the following and I keep receiving error Undefined function 'fixedpoint' for input arguments of type 'function_handle'.
(I'm new in Matlab, so there may be both syntactical or semantical errors...)
function [ x ] = fixedpoint(g,I,y,tol,m)
% input: g, I, y, tol, max
% g - function
% I - interval
% y - starting point
% tol - tolerance (error)
% m - maximal number of iterations
% x - approximate solution
a=I(1);b=I(2);
if(y<a | y>b)
error('The starting iteration does not lie in I.')
end
x=y;
gx=g(y);
while(abs(x-gx)>tol & m>0)
if(gx<a | gx>b)
error('The point g(x) does not lie in I.')
end
y=x;
x=g(y);
m=m-1;
end
0 comentarios
Respuesta aceptada
Dimitris Kalogiros
el 22 de Sept. de 2019
Dear John
Put your function into the same folder with the program (m-file) that calls it.
4 comentarios
Dimitris Kalogiros
el 24 de Sept. de 2019
Hm hm, I don't know if it is necessary , but I always follow this rule.
Más respuestas (4)
emmanuel john Lavarias
el 27 de Sept. de 2021
- Solve one real root of ex−2x−5=0ex−2x−5=0 with x0=−2x0=−2 using the Fixed-Point Iteration Method accurate to four decimal places.
0 comentarios
Ahteshamul Hoque Tareq
el 8 de En. de 2022
a=I(1);b=I(2); if(y<a | y>b) error('The starting iteration does not lie in I.') end x=y; gx=g(y); while(abs(x-gx)>tol & m>0) if(gx<a | gx>b) error('The point g(x) does not lie in I.') end y=x; x=g(y); m=m-1; end
0 comentarios
Tsega'ab
el 13 de Dic. de 2023
function [ x ] = fixedpoint(g,I,y,tol,m)
% input: g, I, y, tol, max
% g - function
% I - interval
% y - starting point
% tol - tolerance (error)
% m - maximal number of iterations
% x - approximate solution
a=I(1);b=I(2);
if(y<a | y>b)
error('The starting iteration does not lie in I.')
end
x=y;
gx=g(y);
while(abs(x-gx)>tol & m>0)
if(gx<a | gx>b)
error('The point g(x) does not lie in I.')
end
y=x;
x=g(y);
m=m-1;
end
0 comentarios
Ver también
Categorías
Más información sobre Loops and Conditional Statements 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!