20.02.2013 Views

Atomic Structure Theory

Atomic Structure Theory

Atomic Structure Theory

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

262 Solutions<br />

From the above, it follows<br />

2 · l ···2 · 1<br />

Il =<br />

(2l +1)···5 · 3 2=<br />

2l l!<br />

3 · 5 ···(2l +1) 2<br />

<br />

l 2<br />

2 l!<br />

= 2 .<br />

(2l + 1)!<br />

Therefore,<br />

<br />

1<br />

c = =<br />

Il<br />

1<br />

2l <br />

(2l + 1)!<br />

.<br />

l! 2<br />

1.3 maple program to generate first 10 Legendre polynomials using Rodrigues’<br />

formula:<br />

for l from 0 to 9 do<br />

if l>0 then p := expand(diff((x^2-1)^l,x$l)/(2^l*l!))<br />

else p :=1 fi;<br />

print(P[l]=p);<br />

od;<br />

mathematica program,<br />

P = {1, Table[Expand[D[(x^2 - 1)^l, {x, l}]/(2^l l!)],<br />

{l, 1, 9}]} // Flatten // TableForm<br />

1.4 mathematica program to generate first 10 Legendre polynomials using<br />

the recurrence relation:<br />

For[P[0] = 1; P[1] = x; l = 2, l < 11,<br />

P[l] = Expand[((2l - 1)x P[l - 1] - (l - 1)P[l - 2])/l];<br />

Print["P[", l, "]=", P[l]]; l++]<br />

1.5 maple routine to generate P m<br />

l (x) forl≤4and1≤m≤ l:<br />

for l from 1 to 4 do<br />

for m from 1 to l do<br />

q:=(sqrt(1-x^2))^m*<br />

diff(expand(diff((x^2-1)^l,x$l)/(2^l * l!)),x$m);<br />

print(P[l,m]= q);<br />

od;<br />

od;<br />

1.6 mathematica routine to find j6(x) using the upward recurrence relation:<br />

j[0] = Sin[x]/x;<br />

j[1] = Sin[x]/x^2 - Cos[x]/x;<br />

Do[j[n + 1] = Expand[(2n + 1) j[n]/x - j[n - 1]], {n, 1, 5}]<br />

j[6]

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!