22 lines
787 B
Plaintext
22 lines
787 B
Plaintext
print "First estimate is: ", using( "#.###############", NthRoot( 125, 5642, 0.001 ));
|
|
print " ... and better is: ", using( "#.###############", NthRoot( 125, 5642, 0.00001))
|
|
print "125'th root of 5642 by LB's exponentiation operator is "; using( "#.###############", 5642^(1 /125))
|
|
|
|
print "27^(1 / 3)", using( "#.###############", NthRoot( 3, 27, 0.00001))
|
|
print "2^(1 / 2)", using( "#.###############", NthRoot( 2, 2, 0.00001))
|
|
print "1024^(1 /10)", using( "#.###############", NthRoot( 10, 1024, 0.00001))
|
|
|
|
wait
|
|
|
|
function NthRoot( n, A, p)
|
|
x( 0) =A
|
|
x( 1) =A /n
|
|
while abs( x( 1) -x( 0)) >p
|
|
x( 0) =x( 1)
|
|
x( 1) =( ( n -1.0) *x( 1) +A /x( 1)^( n -1.0)) /n
|
|
wend
|
|
NthRoot =x( 1)
|
|
end function
|
|
|
|
end
|