11
As I understood you start from the completeness relation
$$\sum_{\ell=0}^\infty \frac{2\ell + 1}{2} P_\ell(x)P_\ell(y) = \delta(x-y)$$
and use that
$$
P_n(0) =
\begin{cases}
\frac{(-1)^{m}}{4^m} \tbinom{2m}{m}
= \frac{(-1)^{m}}{2^{2m}} \frac{(2m)!}{\left(m!\right)^2}
& \text{for} \quad n = 2m
\\
0 & \text{for} \quad n = 2m+1 \,....
8
I think you're misinterpreting Mathematica output. First of all, notice that the first output you get from Mathematica is a message that precludes the subsequent numerical results being sensible (or at the very least must make you very cautious of them):
Sum::div: Sum does not converge.
If you omit the N[] wrapper, you'll get the actual sum simplified to ...
5
We need some numerical model for comparison, so this is one of them based on FEM. First we make sufficient mesh for this problem:
Needs["NDSolve`FEM`"];Needs["MeshTools`"];
L = .90; l = 1.80; w = 0.0003; bh = 17.394;
bc = 22.151; ph = 8.6;
pc = 13.93; pa = 10; n = 10;
thi = 460; tci = 300; Ta = 380; region = Rectangle[{0, 0}, {L, l}];
...
5
Your T function is on the left-hand side dependent on {x,y,z} but on the right-hand side not a y in the MathML code. You got confused by the name of the functions in special states of the solutions process and forget to use them consequently. The solution of the Subscript[C,1], Subscript[C,2] depends in length on the given parameters but all are not set in ...
3
I'm not sure if the definitions you posted are different from the ones used in your plot, or you have lingering definitions that are interfering, but the plot I get from your code is slightly different than the one you show.
Z1D[\[Beta]_, \[Xi]_, \[Theta]_] = (1/\[Xi] (2/\[Beta]^2 -
1) + ((\[Theta] Exp[-\[Beta]])/2 - 1/2) +
Sum[((-\[Beta])^n \[...
2
As a starting point, I suggest to look at this post
https://math.stackexchange.com/questions/1108246/double-sum-and-zeta-function
We can at least verify numerically that the sum is equal to
$$S(3)=\frac1{\Gamma(3)}\int_0^{\infty}\!\! t^2 \big(\theta_3(0,e^{-t})^2-1\big)\, \mathrm{d}t$$
With $\Gamma(3)=2$, we have
NIntegrate[t^2(EllipticTheta[3,0,Exp[-t]]^2-1)...
1
Clear["Global`*"]
f[x_] = x^2;
a[0] = 2*Integrate[f[x], {x, 0, 1}];
a[k_] = Assuming[Element[k, Integers],
2*Integrate[f[x]*Cos[2*k*Pi*x], {x, 0, 1}]]
(* 1/(k^2 π^2) *)
b[k_] = Assuming[Element[k, Integers],
2*Integrate[f[x]*Sin[2*k*Pi*x], {x, 0, 1}]]
(* -(1/(k π)) *)
F[x_, n_Integer?Positive] :=
a[0] + Sum[a[k - 1]*Cos[2*k - 1*Pi*x] + b[k]*Sin[2*...
Only top voted, non community-wiki answers of a minimum length are eligible
