laguerre.m

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

function L = laguerre(k, n, x);
    
    if n < 0
        error('n nie moze byc < zero');
    end

    L0 = ones(size(x));
	if n == 0
    	L = L0;
    end

    L1 = k + 1 - x;
    if  n == 1
    	L = L1;
    end

    Lm2 = L0; 
    Lm1 = L1;

	for j = 2 : n
        L = 2 * Lm1 - Lm2 + ((k - 1 - x) .* Lm1 - (k - 1) * Lm2) / j;
	    Lm2 = Lm1; 
        Lm1 = L;
	end