Re: Cube root

Hmm.. Good question.  I can find the C function SQRT, but that calculates
the square root of a number.

I guess you could do in mathematically, but the only way I could figure out
how to calculate the square root would take a bit of time.  I guess it depends
on how accurate you need it to be, and how fast you need it to be.  An algorithm
would be something like:

If Number = 0 return 0
If Number < 0 return -1 error (can't return root of negative number)
Answer = Number
LastAnswer = Number
HighAnswer = Number
LowAnswer = 0

dow   Answer * 3 <> Number

If    Answer * 3 < Number
eval  Answer = Answer + (HighAnswer - Answer) / 2  (Split difference)
Eval  HighAnswer = LastAnswer
Else
Eval  Answer = Answer - (Answer - LowAnswer) / 2  (Split Difference)
Eval  LowAnswer = LastAnswer
EndIf

If   LastAnswer = Answer return Answer (Close as we can get with precision)

Eval  LastAnswer = Answer

endDo

Return Answer

This is off the top of my head and I'm sure this pseudo code contains bugs, but
you should get the idea.  Also, you can find the arbitrary root of a number this
way also.  Just replace the * 3 with * Root and pass it the root too.

I am not sure how many loops will occur, but using this method a number can be
found from 1 to 100 in a maximum of 7 splits.  That's 3 digits.  But it's a
geometric progression, not a linear one, so from 1 to 1000 takes only 10 splits,
1 to 10000 takes 15 splits, etc..

I am not suggesting you do it this way, by any means, just that this way should
work if you can't find another solution.

Regards,

Jim Langston

> Walter Hesius wrote:
> 
> Hi all,
> 
> I would like to know a way to calculate the cube root of a number. We are on 
>V4R4.
> 
> Thanx.
+---
| This is the RPG/400 Mailing List!
| To submit a new message, send your mail to RPG400-L@midrange.com.
| To subscribe to this list send email to RPG400-L-SUB@midrange.com.
| To unsubscribe from this list send email to RPG400-L-UNSUB@midrange.com.
| Questions should be directed to the list owner/operator: david@midrange.com
+---