banach.m

来自「Nonlinear Problem Files Collection」· M 代码 · 共 53 行

function [] = banach(f,x0,maxite);
%
% banach berechnet und visualisiert eine Fixpunktiteration
% Aufruf: banach(f,x0,maxite)
% Parameter: Funktionsstring f, Startwert x0, 
%            maximale Iterationsanzahl maxite
%
t=0:0.05:10;
plot(t,t);

axis([0 5 0 5]);
hold on;

plot(t,feval(f,t));

x(1)=x0; k=0;
y(1)=feval(f,x0);
disp(sprintf('x_%d=%f',k,x(k+1)));

plot(x(1),y(1),'*');

for k=2:2:2*(maxite-1)
y(k)=y(k-1);
x(k)=y(k);
disp(sprintf('x_%d=%f',k/2,x(k)));
plot(x,y,'*-');
waitforbuttonpress;
y(k+1)=feval(f,x(k));
x(k+1)=x(k);
plot(x,y,'*-');
waitforbuttonpress;
end;

k=2*maxite;
y(k)=y(k-1);
x(k)=y(k);
disp(sprintf('x_%d=%f',k/2,x(k)));
plot(x,y,'*-');

hold off;