Feigenbaum constant calculation: Difference between revisions

m
→‎{{header|REXX}}: optimized the inner two DO loops.
m (→‎{{header|REXX}}: added the number of correct decimal digits in the on-going calculation of Feigenbaum's constant.)
m (→‎{{header|REXX}}: optimized the inner two DO loops.)
Line 483:
do i=2 for maxi-1
a= a1 + (a1 - a2) / d1
do j=1 for maxj
x= 0; y= 0
do k=1 for 2**i
y= 1 - 2 * x * y
x= a - x**2
end /*k2**i*/
a= a - x / y
end /*jmaxj*/
d= (a1 - a2) / (a - a1)
t= max(0, compare(d, #) - 2) /*# true digs so far, ignore dec. point*/