view release on metacpan or  search on metacpan

examples/krzysztof_reformulated_zeta_function.pl  view on Meta::CPAN

use 5.018;
use strict;
use warnings;

use experimental qw(signatures);

use lib qw(../lib);
use Math::AnyNum qw(:overload);    # can be commented out
use Math::AnyNum qw(gamma binomial zeta factorial);

sub Ak($k) {
    my $sum = 0;
    foreach my $j (0 .. $k) {
        $sum += (-1)**$j * binomial($k, $j) * (2 * $j + 1) * zeta(2 * $j + 2);
    }
    $sum;
}

sub krzysztof_zeta ($s, $r = 100) {
    my $sum = 0;
    foreach my $k (0 .. $r) {