Сформировать массив из среднегеометрических значений положительных элементов каждого столбца матрицы - C++

int mnogomer (double (&arr)[5][5]) 
{
double sum0=0, sum1=0, sum2=0, sum3=0, sum4=0;
int k0=0, k1=0, k2=0, k3=0, k4=0;
double p0=1.0, p1=1.0,p2=1.0, p3=1.0, p4=1.0, k=0;
for (int i=0; i<5; i++)
   for(int j=0;j<5; j++)
       if (arr[i][j]>0)
           {if ((j==0)&&(arr[i][j]>0))
                 { 
                  sum0+=arr[i][j];
                  k0++;
                 }
           else if ((j==1)&&(arr[i][j]>0))
                 { 
                  sum1+=arr[i][j];
                  k1++;
                 }
           else if ((j==2)&&(arr[i][j]>0))
                 { 
                  sum2+=arr[i][j];
                  k2++;
                 }
           else if ((j==3)&&(arr[i][j]>0))
                 {
                  sum3+=arr[i][j];
                  k3++;
                 }
           else if ((j==4)&&(arr[i][j]>0))
                 { 
                  sum4+=arr[i][j];
                  k4++;
                 }    
           }
 
 
double mas0= pow( p0, 1.0/k0 );
double mas1= pow( p1, 1.0/k1 );
double mas2= pow( p2, 1.0/k2 );
double mas3= pow( p3, 1.0/k3 );
double mas4= pow( p4, 1.0/k4 );
 
return mas0, mas1, mas2, mas3, mas4;
 
}
 
int _tmain(int argc, _TCHAR* argv[])
{
double  arr[5][5]; 
double mas0;
double mas1;
double mas2;
double mas3;
double mas4;
 
for (int i = 0; i < 5; i++)
        for (int j = 0; j < 5; j++)
            arr [i][j] = rand() % 10;
                        
        
for (int i = 0; i <5; i++)
     {for(int j = 0; j <5; j++)
         cout<<arr[i][j];
         
     cout<<endl;
     }
mnogomer(arr);
double aray[5]={mas0, mas1, mas2, mas3, mas4};
for (int i=0; i<5; i++)
     cout<<aray[i]<<" ";
 
}