# Plotting double series

I am new and hoping a warm welcome from this platform.

I am trying to plot the graph of double series in one variable

$$\sum_{k=0}^{\infty} \sum_{j=0}^{\infty} C_{k,j} \exp(- 3^k 1.5^{j} x)$$

where coefficients are defined as \begin{align} C_{0,0} &= 4.55672\\ C_{k,0} &= \frac{(-1)^k 3^k}{2^k \prod_{s=1}^{k}(3^s -1)}C_{0,0}\\ C_{0,j} &= \frac{(-1)^j 1.5^j}{2^j \prod_{s=1}^{j}(1.5^s-1)}C_{0,0}\\ C_{k,j} &= \frac{(-1) (3 C_{k-1,j}+1.5 C_{k,j-1})}{2(3^k 1.5^j -1)} \end{align}

So please tell me how to define these coefficients and can be recalled in the double sum.

• Can you meet us halfway and write all these expressions in correct Mathematica syntax? As a next step, read defining functions and look up Product and Sum. Try small examples and become comfortable with these functions before tackling the big problem! Please try on your own first, based on the resources I pointed to, and then make the question specific: ask about the first speific issue you encountered. If you've never used Mathematica, then there will be a learning curve and this is too big for your very first task. – Szabolcs Mar 15 at 12:40
• I should note that I do not expect this to be a completely trivial exercise. Things you are likely to hit: defining recursive functions, memoization, a need for NumericQ in function definition, dealing with numerical underflow and truncating the series appropriately to get a good approximation. Explaining all this is too much for one answer so we either need to know that you understand these topic, or proceed step by step with them. – Szabolcs Mar 15 at 12:53
• Maybe I would proceed step by step, as am not much familiar about these topics. – John Mar 16 at 3:48
• Note that Product[a^s - 1, {s, k}] gives (-1)^k QPochhammer[a, a, k]. – Αλέξανδρος Ζεγγ Mar 17 at 5:48
• Just curious: is this a probability density function? – JimB Mar 17 at 6:24

As an example for beginners, I give code using a block. The series converges quickly, so we can take a finite number of members (10 in this example)

F[n0_, x0_] :=
Block[{n = n0, x = x0},
f[n_, x_] := Sum[c[k, j]*Exp[-3^k*(3/2)^j*x], {k, 0, n}, {j, 0, n}];
c[0, 0] = 4.55679;
c[k_, 0] := (-1)^k*3^k*c[0, 0]/2^k/Product[3^s - 1, {s, 1, k}];
c[0, j_] := (-1)^j*(3/2)^j*
c[0, 0]/2^j/Product[(3/2)^s - 1, {s, 1, j}];
Table[c[k,
j] = (-1)*(3*c[k - 1, j] + 3/2*c[k, j - 1])/
2/(3^k*(3/2)^j - 1), {k, 1, n}, {j, 1, n}];
f[n, x]]


In fig. function F[n,x] is shown for various n=1,2,3,4,5,6,10

{Plot[Evaluate[Table[F[k, x], {k, 1, 6}]], {x, 0, 5},
PlotRange -> All, PlotLegends -> Automatic],
Plot[F[10, x], {x, 0, 5}, PlotRange -> All]} 