30 lines
339 B
Plaintext
30 lines
339 B
Plaintext
maxit = 13;
|
|
maxitj = 10;
|
|
a1 = 1.0;
|
|
a2 = 0.0;
|
|
d1 = 3.2;
|
|
a = 0.0;
|
|
Table[
|
|
a = a1 + (a1 - a2)/d1;
|
|
Do[
|
|
x = 0.0;
|
|
y = 0.0;
|
|
Do[
|
|
y = 1.0 - 2.0 y x;
|
|
x = a - x x;
|
|
,
|
|
{k, 1, 2^i}
|
|
];
|
|
a = a - x/y
|
|
,
|
|
{j, maxitj}
|
|
];
|
|
d = (a1 - a2)/(a - a1);
|
|
d1 = d;
|
|
a2 = a1;
|
|
a1 = a;
|
|
{i, d}
|
|
,
|
|
{i, 2, maxit}
|
|
] // Grid
|