blob: f4031a234872b4082e2483c56b078d769b49b1bb (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
function [x]=fixedp(x0,f)
//fixed-point iteration
N = 100; eps = 1.e-5; // define max. no. iterations and error
maxval = 10000.0; // define value for divergence
a = x0;
while (N>0)
xn = f(a);
if(abs(xn-a)<eps)then
x=xn
disp(100-N);
return(x);
end;
if (abs(f(x))>maxval)then
disp(100-N);
error('Solution diverges');
abort;
end;
N = N - 1;
xx = xn;
end;
error('No convergence');
abort;
//end function
|
