radial.m

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

function  R = radial(n, L, r);

    if L >= n
        error('N musi byc mniejsze niz L');
    end
    if L < 0 || n < 0
        error('N oraz L musza byc dodatnie');
    end

    % odwracamy jesli dane sa typu spektralnego aby bylo mozna dobrze mnozyc
    [nw nk] = size(r);
    if nk > nw
        r = r';
    end

    y = 2 * r ./ n;
    z = exp(-0.5 * y);
    
    k = 2 * L + 1; 
    p = n - L - 1;
    
    wsp_norm = 2 * sqrt(prod(1 : n - L - 1) / prod(1 : n + L)) / n^2;
    R = wsp_norm * (y.^L) .* z .* laguerre(k, p, y);