# # Written by: # -- # John L. Weatherwax 2009-04-21 # # email: wax@alum.mit.edu # # Please send comments and especially bug reports to the # above email address. # #----- a = 0.8 b = 1 print( c( a, b ) ) for( n in 1:5 ){ aprime = 0.5 * ( a + b ) bprime = ( 2*a*b ) / ( a+b ) print( c( aprime, bprime ) ) a = aprime b = bprime }