Advertisement
apexsquirt

AAAAAAAAA

Apr 20th, 2018
213
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
PHP 0.53 KB | None | 0 0
  1. <?php
  2.  
  3. $sup = readline("Sum upper bound : ");
  4.  
  5. aaa:
  6. print "\n";
  7.     $s = readline("Real part > 0  : ");
  8.     $t = readline("Imaginary part : ");
  9.     bbb:
  10.     $re = 0;
  11.     $im = 0;
  12.     for ($i = 1; $i <= $sup; $i++) {
  13.         $re -= cos(pi()*$i)*cos($t*log($i))/pow($i,$s);
  14.         $im += cos(pi()*$i)*sin($t*log($i))/pow($i,$s);
  15.     }
  16.     if ($t == 0) { $im = 0; }
  17.     if ($im < 0) { $sgnim = "-"; } else { $sgnim = "+"; }
  18.     print "DirichletEta($s + {$t}i) ≈ $re ";
  19.     if ($im != 0) { print "$sgnim ".abs($im)."i\n\n"; }
  20.     else { print "\n\n"; }
  21.     goto aaa;
  22. stop:
  23.  
  24. ?>
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement