gsl_matrix *d = gsl_matrix_calloc(n,n);
    for(i = 0; i < n; i++)
    {
        for(j = 0; j < n; j++)
        {
                value = 0.0;
                gsl_matrix_set(d, i, j, 0.0);
                for(int k = 0; k < n; k++)
                {
                        value += 1./(i+k+2.) * gsl_matrix_get(c, k, j);
                }
                gsl_matrix_set(d, i, j, value);
        }
    }
    fp = fopen("d.dat", "w");
    for(i = 0; i < n; i++)
    {
        for(j = 0; j < n; j++) fprintf(fp, "%12.4g", gsl_matrix_get(d, i, j));
        fprintf(fp, "\n");
    }
    fclose(fp);