function D=JCHZH(n,x,u)
%使用基于基的運算的插值法求控制頂點
%%該算法取兩端點節點重度r=k+1
D=zeros(1,n+2);
DT=zeros(n+2,n+2);
for j=1:2
    DT(1,j)=(-1)^j;
end
    DT(n+2,n+2)=1;
for i=2:n+1
    for j=i-1:i
        DT(i,j)=Ji2(j-1,i,u);
    end
end

P=zeros(1,n+2);
P(1,1)=(u(4)-u(3))/2;
for i=2:n+2
    P(1,i)=x(i-1);
end

D=(inv(DT)*(P'))';