Processing math: 100%

First time here? Check out the FAQ!

Ask Your Question
0

How to calculate L(1,χ)/L(1,χ)?

asked 13 years ago

anonymous user

Anonymous

Question as in title, where L(s,χ) is the Dirichlet L-function associated with the nontrivial character modulo 3. Please provide complete SAGE code. Thank you in advance.

Preview: (hide)

Comments

Maybe you could start by showing us how much you have already?

niles gravatar imageniles ( 13 years ago )

Have you had a look at the [L-function tutorial](http://wiki.sagemath.org/days33/lfunction/tutorial) from Sage days 33? In particular, it has some fairly recent comments on the development status. I have no idea whether Sage yet implements, for example, the insights in Ihara, Murty & Shimura's paper from circa 2007, [On the Logarithmic Derivative of Dirichlet L-Functions at s=1](http://www.kurims.kyoto-u.ac.jp/~kenkyubu/emeritus/ihara/Publications-and-Recent-Preprints/RecentArticles/pdf-files/IMS.main.pdf).

bgins gravatar imagebgins ( 13 years ago )

Actually I managed to calculate the particular example "by hand", using that L(1,χ)=χ(n)/n and L(1,χ)=χ(n)log(n)/n. Of course it would be nicer to use some built-in function for that purpose. Note that I am a beginner at SAGE. Also, I looked at http://wiki.sagemath.org/days33/lfunction/tutorial, but the command "LSeries" did not work at http://www.sagenb.org/ (upon calling "L=LSeries(DirichletGroup(3).0)" I get the error message "name 'LSeries' is not defined").

Anonymous gravatar imageAnonymous ( 13 years ago )

1 Answer

Sort by » oldest newest most voted
0

answered 8 years ago

dan_fulea gravatar image

We can try to compute using bare hands. The following is thus not useful if the question wants more than this particular case. I cannot see how sage supports better this now.

First of all, we can compute L(χ,1)=n0(13n+113n+2)=n01(3n+1)(3n+2) exactly, for instance:

sage: sum( 1/(3*n+1)/(3*n+2), n, 0, oo )                           
1/9*pi*sqrt(3)

and the result is connected to "polylogarithmic computations".

Comment: We may start with 1+x3+x6+=1/(1x3) and integrate twice, first with x from 0 to y, then with y from 0 to 1. Fubini shows we can forget about polylogarithms, since L(χ,1)=10dyy0dx11x3=10dx1xdy11x3=10dx1+x+x2=13π .

The derivative is more complex. Pari/GP gave

? sum( n=0,10000000, -log(3*n+1)/(3*n+1)+log(3*n+2)/(3*n+2) )
%6 = 0.2226631782653383756620209560

but suminf( n=0,-log(3*n+1)/(3*n+1)+log(3*n+2)/(3*n+2) ) was testing my patience.

At any rate, the rest can be estimated by rewriting the sum as a sum over 1(3n+1)(3n+2)[ ln(1+13n+1)3n+2ln(3n+2) ] or over 1(3n+1)(3n+2)[ ln(1+13n+1)3n+1ln(3n+1) ] .

Preview: (hide)
link

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower

Stats

Asked: 13 years ago

Seen: 816 times

Last updated: Mar 13 '17