Считать размерность матрицы из файла

#include <iostream>
#include <string>
 
using namespace std;

#include <stdio.h>
#include <string.h>
#include <istream>
using namespace std;
void main(void)
{
    int N=3,M=0; 
    int a[50][50],b[50][50]; //исходные матрицы
    int c[50][50],d[50][50],e[50][50]; //результат
    FILE *in=fopen("file.txt", "r");
    FILE *out=fopen("out.txt", "w");
    char str[100],str2[100];
    fgets(str,100,in);
    for (int i = 0; i < strlen (str); i++) //алгоритм подсчета ширины
    {
        if ( str[i] == ' ' )
        {
            M++;
        }
    }
    M++;
    printf("Hor: %d\n",M);
    delete (str);
    rewind(in);
 
    while (str!=NULL) //алгоритм подсчета высоты - не пашет
    {
        //fgets (str2,100,in);
        getline(in,str2,'\n');
        N++;
        printf("!");
    }
    N++;
    printf("Vert: %d\n",N);
    delete (str2);
    rewind(in);     //написано по подобию этого кода:
    
                                      ////        int i=0;
                                      ////        ifstream base("name.txt");
                                      ////        while (!base.eof())
                                      ////        {
                                      ////                base.getline(str, 1024, '\n');
                                      ////                i++;
                                      ////        }
                                      ////        base.close();
                                      ////        delete str;
                                      ////        cout << i << '\n';
                                      ////        system ("pause");
                                      ////}                                 
    //fscanf(in,"%d",&N);//vert
    //fscanf(in,"%d",&M);//hor
    for(int i=0;i<N;i++) //считывание матриц уже известных размеров
    {
        for(int j=0;j<M;j++)
        {
            fscanf(in,"%d",&a[i][j]);
        }
    }
    for(int i=0;i<N;i++)
    {
        for(int j=0;j<M;j++)
        {
            fscanf(in,"%d",&b[i][j]);
            c[i][j]=a[i][j]+b[i][j]; //действия над ними, дальше вывод результатов
            d[i][j]=a[i][j]-b[i][j];
            e[i][j]=b[i][j]-a[i][j];
        }
    }
    fprintf(out,"\na+b\n\n");
    for(int i=0;i<N;i++)
    {
        for(int j=0;j<M;j++)
        {
            fprintf(out,"%d ",c[i][j]);
        }
        fprintf(out,"\n");
    }
    fprintf(out,"\na-b\n\n");
    for(int i=0;i<N;i++)
    {
        for(int j=0;j<M;j++)
        {
            fprintf(out,"%d ",d[i][j]);
        }
        fprintf(out,"\n");
    }
    fprintf(out,"\nb-a\n\n");
    for(int i=0;i<N;i++)
    {
        for(int j=0;j<M;j++)
        {
            fprintf(out,"%d ",e[i][j]);
        }
        fprintf(out,"\n");
    }
}