for(0..@c-1) {
#   $lambda[$_] = abs($c[$_]-($rho*sqrt($c[$_])));
    $lambda[$_] = abs($c[$_]-($rho * Math::BigFloat->new($c[$_])->bsqrt()));
}