p>>
plot3d ([r * cos (theta), r * sin (theta), (r * R (d) * Theta (d)) ^ 2], < br> r = 0 .. 20, theta = 0 .. 2 * Pi, axes = boxed, orientation = [0,0],
shading = z, style = patchcontour, scaling = constrained,
title = ` Контурна проекція квадрата норми ` ); В
Визначимо процедуру, яка дозволяє побудувати всі розглянуті вище графіки для якої-небудь з гармонік:
> HydrogenPlots: = proc (n, l, m) global a, p1, p2, p3, p4, p5; local txt;
a: = 1; txt: = `nlm =` | | n | | l | | m:
p1: = sphereplot (Theta (d) ^ 2, phi = Pi/2..2 * Pi, theta = 0 .. Pi, axes = boxed,
scaling = constrained, grid = [15,100], title = `txt`); print (p1);
p2: = polarplot ([Theta (d) ^ 2, theta + Pi/2, theta = 0 .. 2 * Pi], scaling = constrained,
title = `txt`); print (p2);
p3: = plot ((r * R (d)) ^ 2, r = 0 .. 30, title = `txt`); print (p3);
p4: = plot3d ([r * cos (theta), r * sin (theta), (r * R (d) * Theta (d)) ^ 2],
r = 0 .. 30, theta = Pi/2..2 * Pi, axes = boxed, title = `txt`); print (p4);
p5: = plot3d ([r * cos (theta), r * sin (theta), (r * R (d) * Theta (d)) ^ 2],
r = 0 .. 30, theta = 0 .. 2 * Pi, axes = boxed, orientation = [0,0], style = patchcontour, scaling = constrained, shading = z, title = `txt`);
end:
Наприклад, нехай
> n: = 4: l: = 2: m: = 1:
HydrogenPlots (n, l, m);
В В В В В
Звичайно, вигляд графіків можна змінити, наприклад змінивши стиль
> replot (p1, style = patch, shading = z, orientation = [56,70]);
В
Можна подивитися, як змінюється розподіл ймовірності залежно від номера гармоніки і без анімації, наприклад при одній і тій же енергії
> for n from 4 to 4 do
for l from 0 to n-1 do
for m from 0 to l do
txt: = `nlm =` | | n | | l | | m:
p | | n | | l | | m: = plot3d ([r * cos (theta), r * sin (theta), (r * R (d) * Theta (d)) ^ 2],
r = 0 .. 50, theta = 0 .. 2 * Pi, axes = boxed, orientation = [0,0],
style = contour, scaling = constrained, shading = z,
title = `txt`, numpoints = 1000);
print (p | | n | | l | | m);
od; od; od;
В В В В В В В В В В В
Список літератури
В
Мінкін В.І., Сімкін Б.Я., Миняев PM p> Теорія буд...
