## 184 Reputation

16 years, 167 days

## Some progress...

OK, I'm not all that brilliant with Maple, here is what I have so far:

Rx := x -> sum( (ln(x))^k/(k*k!*evalf(Zeta(k+1))) ,k=1..100 ): Tk := (x,k) -> - ( Rx( x^( 0.5 + I*Zeroes[k] ) ) + Rx( x^( 0.5 - I*Zeroes[k] ) ) ): Rn:= (x,n) -> Rx(x) - sum( Rx( x^( -2*m ) ), m = 1..100 ) + sum( Tk(x,j), j=1..n ):

Now,  I trimmed infinity to 100 for computing R(x).  'Zeroes' is a list of zeroes. How would I go about plotting this?  I tried

complexplot(Rn(x,20), x=1..40);

But it is taking ages (yes, i used with(plots)).  Could someone who knows what they're doing help me out please?

## Some progress...

OK, I'm not all that brilliant with Maple, here is what I have so far:

Rx := x -> sum( (ln(x))^k/(k*k!*evalf(Zeta(k+1))) ,k=1..100 ): Tk := (x,k) -> - ( Rx( x^( 0.5 + I*Zeroes[k] ) ) + Rx( x^( 0.5 - I*Zeroes[k] ) ) ): Rn:= (x,n) -> Rx(x) - sum( Rx( x^( -2*m ) ), m = 1..100 ) + sum( Tk(x,j), j=1..n ):

Now,  I trimmed infinity to 100 for computing R(x).  'Zeroes' is a list of zeroes. How would I go about plotting this?  I tried

complexplot(Rn(x,20), x=1..40);

But it is taking ages (yes, i used with(plots)).  Could someone who knows what they're doing help me out please?

## Thanks...

I haven't tried it no.   The obvious way to try it is to define an array of zeroes then iterate through a for loop.  Could I simply define R(x) in Maple as I have done in my previous post, or would it be overkill?  I suppose I could use stirling's approx for the factorial.

Yes, I don't see me purchasing Mathematica just to do a few graphs.  It's annoying since my university has Maple, Matlab, Minitab, etc but I haven't come across systems with Mathematica installed.  And it really does seem ideal for this.