c - c中的递归多项式族

Question

我有一个递归定义的多项式族：

T0(x)=1
T1(x)=2x−1,
Tk(x)=2(2x−1)T(k−1)(x)−T(k−2)(x),∀k= 2,3,...

我要做的是定义一个函数来计算多项式的系数并将其存储在矩阵中。

示例：+1.0000e+00

-1.0000e+00 +2.0000e+00

+1.0000e+00 -8.0000e+00 +8.0000e+00

我的尝试是：

if(k==0){
    c[0][0]=1.;
}
if(k==1){
    c[0][0]=1.;
    c[1][0]=-1.;
    c[1][1]=2.;
}
if(k>1){ 
    f(k-1,c);
    c[k][k-1]=-2*c[k-1][k-1];
    for(i=0; i<k-1; i++){
        c[k][i]=-(c[k-2][i]+2*c[k-1][i]);
    }
    c[k][k]=0;
    for(i=0; i<k; i++){
        c[k][i+1]+=4*c[k-1][i];
    }
}

但我有一个分段错误。我做错了什么？

标签： crecursionpolynomials