Казус при поиске максимального элемента массива.

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <ctype.h>
#include <stdlib.h>
 
using namespace std;
 
int main()
{
    double lamda,delta,alpha,Emax,max,C,eps,Z,C0;
    float W_p,D_p,W_d,D_d,hx,hy;
    long i,j,m,n,p,di,dj,q,Y,Imin,Imax,Jmin,Jmax;
    lamda=1.9;
    delta=1e-3;
    max=0;
    Emax=delta+1;
    C=0;
    C0=0;
    eps=2.2;
 
    //Ввод входных параметров
 
    printf("Hello!  We are welcome program for calculation of wave resistian micro \n");
    printf("W_p=");     //Ширина проводника
    scanf("%lf", &W_p);
    printf("D_p=");     //Длинна проводника
    scanf("%lf", &D_p);
    printf("W_d=");     //Ширина диэлекрика
    scanf("%lf", &W_d);
    printf("D_d=");     //Длинна диэлектрика
    scanf("%lf", &D_d);
    printf("eps=");
    scanf("%lf", &eps);  // Диэлектрическая проницаемость
 
    q=1;
    di=6;
    dj=1;
    hx=D_p/di;      //Шаги сетки по координатам х и у
    hy=W_p/dj;
    n=round((D_d/hx)*q);  //кол-во узлов матрицы по х
    m=round((W_d/hy)*q);  //кол-во узлов матрицы по y
    Y=round(m/2);
    Imin=((n-di)/2)+1;Imax=Imin+di; //координаты проводника
    Jmin=Y;Jmax=Y+1;
    alpha=n/m;
 
    float Ex[n][m],Ey[n][m],E[n][m],B[n][m],P[n][m],D[n][m],A[n][m];
 
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;j++)
        A[i][j]=0;
    }
 
    for(i=Jmin;i<=Jmax;i++)
    {
        for(j=Imin;j<Imax;j++)
         A[i][j]=1;
    }
    //
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;j++)
        {
        P[i][j]=A[i][j];
        B[i][j]=A[i][j];
        D[i][j]=A[i][j];
        }
    }
 
//
 
    //Lаплас
    while (Emax>delta)
    {
        for(i=2;i<n;i++)
        {
            for(j=2;j<m;j++)
            {
                if(P[i][j]==0)
                {
                    D[i][j]=(1-lamda)*D[i][j];
                    D[i][j]=D[i][j]+(lamda*(D[i+1][j]+D[i-1][j]+alpha*D[i][j+1]+alpha*D[i][j-1]))/(2.*(1+alpha));  //D[i][j]+(lamda*(D[i+1][j]+D[i-1][j]+alpha*D[i][j+1]+alpha*D[i][j-1]))/(2*(1+alpha));
                }
                else D[i][j]=1;
 
                B[i][j]=fabs(D[i][j]-A[i][j]);
            }
        }
        max=0;
        for(i=1;i<n;i++)    //Максимальный элемент массива
        {
            for(j=1;j<m;j++)
            {
            if (B[i][j]>max)
            {
                max=B[i][j];
                printf("%e", max);
            }
            }
        }
 
        Emax=max;
        for(i=1;i<=n;i++)        //Изменение параметров начальной матрицы
        {
            for(j=1;j<=m;j++)
            {
                A[i][j]=D[i][j];
            }
        }
 
        printf("\n" "%e", Emax);
 
        for(i=1;i<=n;i++)
        {
        printf("\n");
        for(j=1;j<=m;j++)
        printf("% .2e", B[i][j]);
        }
        getchar();
    }
 
    //расчет напряженности поля
 
    for(i=1;i<=(n-1);i++)
    {
        for(j=1;j<=(m-1);j++)
        {
            Ex[i][j]=0.5*(D[i][j]-D[i+1][j]+D[i][j+1]-D[i+1][j+1]);
            Ey[i][j]=alpha*0.5*(D[i][j]-D[i][j+1]+D[i+1][j]-D[i+1][j+1]);
            E[i][j]=Ex[i][j]+Ey[i][j];
        }
    }
     //Расчет погонной емкости.
     for(i=1;i<n;i++)
    {
        for(j=1;j<m;j++)
        {
            C+=eps*pow(fabs(E[i][j]),2);
        }
    }
    for(i=1;i<n;i++)
    {
        for(j=1;j<m;j++)
        {
            C0=C0+pow(fabs(E[i][j]),2);
        }
    }
    //Расчет погонных параметров
    C=C/alpha;
    C0=C0/alpha;
    //
    Z=(120*3.14)/pow(C*C0,1./2);
 
 
 
 
    //Вывод матрицы на экран   /*
      for(i=1;i<=n;i++)
    {
        printf("\n");
        for(j=1;j<=m;j++)
        printf("% .2e", A[i][j]);
    }
 
    printf("\n");
  //  printf("%f,% f,% f,% f", max, Z, C, C0);
    printf("Z=");
    printf("%f", Z);
    printf("\n");
    printf("C=");
    printf("%f", C);
    printf("\n");
    printf("C0=");
    printf("%f", C0);
    printf("\n");
 
    return 0;
}