Дана квадратная матрица A порядка n. Получить матрицу AB; элементы матрицы B вычисляются по формуле

#include <stdlib.h>
#include <stdio.h>
#include <time.h>
#include <signal.h>
#include <math.h>
#include <malloc.h>
#include <cstring>
#include <limits>
 
 
 
int new_array(double *** matrix, unsigned long int w,unsigned long int l){
    *matrix = new double* [w];
    if (!*matrix) return 6;
    for (unsigned long int i=0; i<w; i++){
        *(*matrix+i)=new double [l];
        if (!*(*matrix+i)) return 6;
    }
    return 0;
}
 
int del_array (double *** matrix,unsigned long int w){
    for (unsigned long int i=0;i<w; i++)
        delete[] *(*matrix+i);
    delete[] *matrix;
 
    return 0;
}
 
int print_array_console(double ** matrix, unsigned long int w,unsigned long int l){
    for (unsigned long int i=0;i<w;i++){
        for (unsigned long int j=0;j<l;j++){
            printf("|%lf|",*(*(matrix+i)+j));
        }
        printf("\n");
    }
    return 0;
}
 
int read_array_console(double *** matrix, unsigned long int* w,unsigned long int* l){
    int temp,debugger;
    printf("enter height:\n");
    temp=scanf("%lu",w);
    if (temp!=1) return 11;
    printf("enter length:\n");
    temp=scanf("%lu",l);
    if (temp!=1) return 11;
    debugger=new_array(matrix,*w,*l);
    if (debugger) return debugger;
    for (long unsigned int i=0;i<*w;i++){
        for (long  unsigned int j=0;j<*l;j++){
            printf("Enter element with i=%lu, j=%lu\n",i,j);
            temp = scanf("%lf",(*(*matrix+i)+j));
            if (temp!=1) return 11;
        }
    }
    return 0;
}
 
 
 
int read_array_file(double *** matrix, unsigned long int* w,unsigned long int* l){
    int temp,debugger;
    char fn[255];
    FILE *inp=NULL;
    printf("Enter filename: \n");
    temp=scanf("%s",fn);
    if (temp!=1) return 11;
    inp =fopen(fn,"r");
    if (!inp) return 11;
    temp=fscanf(inp,"%lu %lu",w,l);
    if (temp!=2) return 11;
    debugger=new_array(matrix,*w,*l);
    if (debugger) return debugger;
    for (unsigned long int i=0;i<*w;i++){
        for (unsigned long int j=0;j<*l;j++){
            temp = fscanf(inp,"%lf",(*(*matrix+i)+j));
            if (temp!=1) return 11;
        }
    }
    fclose(inp);
    return 0;
}
 
int print_array_file(double ** matrix, unsigned long int w,unsigned long int l){
    int temp;
    char fn[255];
    FILE *outp=NULL;
    printf("Enter filename: \n");
    temp=scanf("%s",fn);
    if (temp!=1) return 11;
    outp =fopen(fn,"w");
    if (!outp) return 11;
    for (unsigned long int i=0;i<w;i++){
        for (unsigned long int j=0;j<l;j++){
            fprintf(outp,"%lf ",*(*(matrix+i)+j));
        }
        fprintf(outp,"\n");
    }
    fprintf(outp,"\n");
    fclose(outp);
    return 0;
}
 
inline int multiple_element(double ** matrix,double ** matrix2,unsigned long int w,unsigned long int i,unsigned long int j,double* result){
    for (unsigned long int k=0;k<w;k++){
        *result+=(*(*(matrix+i)+k))*(*(*(matrix2+k)+j));
    }
    return 0;
}
 
int matrix_fun (double ** matrix,unsigned long int w,double ** result){
    double **b=0;
    int debugger=0;
    debugger = new_array(&b,w,w);
    if (debugger) return debugger;
    for (unsigned long int i=0;i<w;i++){
        for (unsigned long int j=0;j<w;j++){
            b[i][j]=1/(abs(i+j)+1);
        
        }
    }
 
    for (unsigned long int i=0;i<w;i++){
        for (unsigned long int j=0;j<w;j++){
            multiple_element(matrix,b,w,i,j,(*(result+j)+i));
        }
    }
    del_array(&b,w);
    return 0;
}
 
 
int main_fun(){
    unsigned long int w,s,temp=0,debugger;
    double **matrix=0;
    double **result=0;
    printf("What mode are you want? (type 0 to console (default),1 to file)");
    temp = scanf("%lu",&s);
    if (temp!=1) return 11;
    if (s==1) {
        debugger=read_array_file(&matrix,&w,&temp);
    }
    else {
        debugger=read_array_console(&matrix,&w,&temp);
    }
    if (debugger) return debugger;
    if (w!=temp) return 11;
    debugger = new_array(&result,w,w);
    if (debugger) return debugger;
    matrix_fun(matrix,w,result);
    if (s==1) {
        debugger=print_array_file(result,w,temp);
    }
    else {
        debugger=print_array_console(result,w,temp);
    }
    if (debugger) return debugger;
    return 0;
}
 
int main(){
    int a;
    a=main_fun();
    if (a==6) printf("Memory allocation error");
    if (a==11) printf("Data error");
    return a;
}