test_radial.m

来自「Three dimensional hydrogen atom spherica」· M 代码 · 共 18 行

r = [0 : 50];

N = 3;
L = 0;
M = 0;

rrA = rrr(N, r')
rrA_val = rrA(:, L + 1)';

rrB = radial(N, L, r);

plot(r, rrA_val, r, rrB);
legend(...
    ['Pierwotna [', num2str(N), ' ', num2str(L), ' ', num2str(M), ']'], ...
    ['Zmieniona [', num2str(N), ' ', num2str(L), ' ', num2str(M), ']'] ...
);