двумерный массив

#include"stdafx.h"
#include<iostream>
#include<iomanip>
#include<fstream>
#include<math.h>
using namespace std;
void read_matr(double Y[4][4]){
ifstream f("c://yes.txt");
for(int i=0; i<4; i++){
for(int j=0; j<4; j++){
f>>Y[i][j];}
 }
 f.close();}
void print_matr(double Y[4][4]){   //печатьматрицывконсоль
 cout<<"matrica:"<<endl;
for(int i=0; i<4; i++){
for(int j=0; j<4; j++){
cout<<setw(10)<<Y[i][j];}
 cout<<endl;}
 }
void print_values(char*str, double v){ // выводвфайлполученныхданных //(дляпунктов 1-4, 6)
fstream f1;
f1.open("values.xls", ios::app);
f1<<str<<"\t"<<v<<endl;
f1.close();
}
void print_matr_file(char*f_name, double Y[4][4]){// выводматрицыв //файлсименем f_name
ofstream f2(f_name);
for(int i=0; i<4; i++){
for(int j=0; j<4; j++){
f2<<Y[i][j]<<"\t";
}
f2<<endl;
}
f2.close();
f2.clear();}
 
double function1(double Y[4][4]){
double s=0;
for (int i = 0; i < 4; i++)
{s+=Y[i][i];
}
print_values("1)\tсумма элементов главной диагонали = ",s);
return s;
}
void function2(double Y[4][4]){
double max=-1000; int j_max;
for(int j=0; j<4; j++)
{
if(Y[2][j]>max){
max=Y[2][j];
j_max=j;}}
print_values("2)\tнаибольшийэлемент 3 строки = ",max);
print_values("\tномерстолбца = ",j_max+1);
cout<<" max = "<<max;
cout<<"; nomer stolbca = "<<j_max+1<<endl;
}
void function3(double Y[4][4], double ss){
Y[2][j_max]=ss;
print_matr(Y);
print_matr_file("3zadanie.xls",Y);
}
 
 
 
int main(){
double A[4][4];
 read_matr(A);
 cout<<"ishodnaja matrica:"<<endl;
 print_matr(A);
cout<<"vipolnenie:"<<endl;
// 1 - нахождение суммы всех элементов в матрице
double s1=function1(A);
cout<<"1) summa glavnoy diagomali = "<<s1<<endl;
// 2 - наибольший элемент вглавной диагонали
cout<<"2) naibolshii element v 3 stroke:"<<endl;
function2(A);
// 3 - поменять А[3][j_max] на s1
cout<<"3) pomenjat A[3][j_max] na s1"<<endl;
function3(A, s1);
return 0;
 }