Неправильный подсчет элементов матрицы

Tab_dif Tab1={3, { {&sin, &cos},{&log, &obr}, {&atan, &atan_diff} },
            { {"sin(x)","cos(x)"},{"Ln(x)","1/x"}, {"Atan(x)","1/(1+pow(x))"} }   };
            
Tab_int Tab2={1, { {&ctg, &ctg_int} },
            {  {"ctg(x)","Ln(sin(x))"} }   };
 
double f1(int i,int j);
void print_arr(double **A,int n);
double **CreateArr(int n);
 
double f1(int i,int j)
{
    double sqr = pow((double)i,2)-pow((double)j,2);
    return pow(-1,sqr);
}
 
double obr(double x)
{
    return 1/x;
}
 
double atan_diff(double x)
{
    return 1/(1+pow(x,2));
}
 
double ctg(double x)
{
    return cos(x*M_PI/180)/sin(x*M_PI/180);
}
 
double ctg_int(double x)
{
    return log(sin(x*M_PI/180));
}
 
double Diff(fun_d ff, double x){
    fun_d f1; int i,ii;
    for(i=0;i<Tab1.ndif;i++)
    if (Tab1.Tab_fun[i][0]==ff){ ii=i;  f1=Tab1.Tab_fun[i][1]; }
    printf("\n производная %6s равна %6s \n", Tab1.name_fun[ii][0], Tab1.name_fun[ii][1]);
    return (*f1)(x);
}  
 
double Integ(fun_i ff, double a, double b){
    fun_i f1; int i,ii;
    for(i=0;i<Tab2.nint;i++)
         if (Tab2.Tab_fun[i][0]==ff){ ii=i; f1=Tab2.Tab_fun[i][1]; }
    printf("\n интеграл %6s равна %6s \n", Tab2.name_fun[ii][0], Tab2.name_fun[ii][1]);
    return (*f1)(b) - (*f1)(a);
}  
 
 
 
double **CreateArr(int n)
{
    int i,j;
    printf(" \n Creating matrix A : \n");
    // Создаётся матрица A
    double **A = new double *[n];
    for(i = 0; i < n; i++) A[i] = new double [n];
 
    for(i = 0; i < n; i++)
    {
        for(j = 0; j < n; j++)
        {
        printf("%f %f %f \n ",f1(i+1,j+1),Diff(atan, (i+1)/(j+1)),Integ(ctg, 2,(i+1)+(j+1)));
            A[i][j] = f1(i+1,j+1)*Diff(atan, (i+1)/(j+1))-Integ(ctg, 2,(i+1)+(j+1));
        }
    }
    return A;
}