meditationatae

Just another WordPress.com site

PARI/gp code related to RH criterion, Lagarias: sigma_1(n) asymptotics

The essential PARI/gp code to check my computation is below:

 

? primorial(1)
%462 = 2
? primorial(1)^6
%463 = 64
? m16
%464 = m16
? \u
harmonic =
  (Z)->Euler+psi(Z+1)

primorial =
  (W)->prod(X=1,W,prime(X))

? Qr(W) = sigma(W,1)/(harmonic(W)+log(harmonic(W))*exp(harmonic(W)))
%465 = (W)->sigma(W,1)/(harmonic(W)+log(harmonic(W))*exp(harmonic(W)))
? m16
%466 = m16
? m16 = primorial(1614)*primorial(37)*primorial(11)*primorial(6)*primorial(4)*primorial(3)*primorial(2)^4*primorial(1)^6 ;
? Qr(m16)
%468 = 0.99923459761253613955435440267376833156
?

Advertisements

Written by meditationatae

June 17, 2013 at 10:41 am

Posted in History

Tagged with ,

%d bloggers like this: