Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- <?php
- $sup = readline("Sum upper bound : ");
- aaa:
- print "\n";
- $s = readline("Real part > 0 : ");
- $t = readline("Imaginary part : ");
- bbb:
- $re = 0;
- $im = 0;
- for ($i = 1; $i <= $sup; $i++) {
- $re -= cos(pi()*$i)*cos($t*log($i))/pow($i,$s);
- $im += cos(pi()*$i)*sin($t*log($i))/pow($i,$s);
- }
- if ($t == 0) { $im = 0; }
- if ($im < 0) { $sgnim = "-"; } else { $sgnim = "+"; }
- print "DirichletEta($s + {$t}i) ≈ $re ";
- if ($im != 0) { print "$sgnim ".abs($im)."i\n\n"; }
- else { print "\n\n"; }
- goto aaa;
- stop:
- ?>
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement