Correction > > > stefan a b = a - a * sqrt (1 - b*b / a*a) > should be: stefan a b = a - a * sqrt (1 - b*b / (a*a)) > > *Main> stefan 10 8 > 4.0 > (0.01 secs, 524896 bytes) > Thanks @@i