function [P,T] = ReduceSplOrd(Q,k,n)

l = length(Q);
%Knotenvektor für f
T = (k-(n-1)):((l+1)+(n-1));
m = length(T)-n;

load omega
P=zeros(1,m);
%Bestimmung von p_j
for j = 1:m
    b=0;
    for mu = 0:n-1
        a=0;
        for i=1:l
            B_Wert = bspline(0:n+k,k-j+i-1,2*mu);
            a = B_Wert * Q(i) + a;
        end
        b=a*(-1)^mu*Omega(n,mu+1) +b;
    end
    P(j)=b;
end