Обработка двумерного числового массива

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <windows.h>
double Ymnojenie(double s1[10][10],int n);
double Summa(double s1[10][10],int n);
void OutPutMatrica(double s1[10][10],int n);
void Transponirovat1(double s1[10][10],int n);
using namespace std;
 
int main()
{
 setlocale(LC_CTYPE,"Russian");
 SetConsoleCP(1251);
 SetConsoleOutputCP(1251);
 
 int n=3,i,j,m=3;
 double Ymnojenie1;
 double Sum1;
 double s1[10][10];
 
 s1[0][0]=1.6;   s1[0][1]=1.9;  s1[0][2]=-4.6;  s1[0][3]=0.9;  s1[0][4]=6.1
 s1[1][0]=9.7;   s1[1][1]=-10.0;  s1[1][2]=1.1;  s1[1][3]=5.0; s1[1][4]=-0.6
 s1[2][0]=6.6;    s1[2][1]=7.1;  s1[2][2]=-1.0;   s1[2][3]=-4.0; s1[2][4]=5.9
 s1[3][0]=4.3;    s1[3][1]=-4.5;  s1[3][2]=0.6;   s1[3][3]=-8.5; s1[3][4]=5.2
 
 OutPutMatrica(s1,n);
 cout<<endl;
 
 Sum1=Summa(s1,n);
 cout<<"Сумма всех положительных элементов = "<<Sum1<<endl<<endl;
 
}
 
double Summa(double s1[10][10],int n)
{
    double Summ=0;
    int i,j;
    for (i=0;i<=n;i++) for(j=0;j<=n;j++) if(s1[i][j]>0) Summ=Summ+s1[i][j];
    return Summ;
}
 
    for (r = 0; r <= n; r++, printf("\n"))
        for (c = 0; c < n; c++)
            printf("%5.1f ",s1[r][c]);
}