function y=phi(j,x,X,n);
    // licznik
    for i=0:(j-1);
        L(i+1)=x-X(i+1);
    end
    for i=(j+1):n;
        L(i)=x-X(i+1);
    end
    // dla mianownika
    for i=0:(j-1);
        M(i+1)=X(j+1)-X(i+1);
    end
    for i=(j+1):n;
        M(i)=X(j+1)-X(i+1);
    end
    // dzielenie
    y=prod(L)/prod(M);
endfunction;