{

 double khi[25][10];
 double f[25];
 double khi11p[25];
 double segexp[10]; 
 double esegmaexp[10];
 ofstream fcoef("../3-root/val.txt");
 ifstream f1(Form("../1-code-dr/fort.0303"));
 ifstream f2(Form("../1-code-dr/fort.0307"));
 ifstream f3(Form("../1-code-dr/fort.0311"));
 ifstream f4(Form("../1-code-dr/fort.0315"));
 ifstream f5(Form("../1-code-dr/fort.0319"));
 ifstream f6(Form("../1-code-dr/fort.0703"));
 ifstream f7(Form("../1-code-dr/fort.0707"));
 ifstream f8(Form("../1-code-dr/fort.0711"));
 ifstream f9(Form("../1-code-dr/fort.0715"));
 ifstream f10(Form("../1-code-dr/fort.0719"));
 ifstream f11(Form("../1-code-dr/fort.1103"));
 ifstream f12(Form("../1-code-dr/fort.1107"));
 ifstream f13(Form("../1-code-dr/fort.1111"));
 ifstream f14(Form("../1-code-dr/fort.1115"));
 ifstream f15(Form("../1-code-dr/fort.1119"));
 ifstream f16(Form("../1-code-dr/fort.1503"));
 ifstream f17(Form("../1-code-dr/fort.1507"));
 ifstream f18(Form("../1-code-dr/fort.1511"));
 ifstream f19(Form("../1-code-dr/fort.1515"));
 ifstream f20(Form("../1-code-dr/fort.1519"));
 ifstream f21(Form("../1-code-dr/fort.1903"));
 ifstream f22(Form("../1-code-dr/fort.1907"));
 ifstream f23(Form("../1-code-dr/fort.1911"));
 ifstream f24(Form("../1-code-dr/fort.1915"));
 ifstream f25(Form("../1-code-dr/fort.1919"));
 ifstream fexp(Form("../../cross-section-q2-05-oct2013/Results/section1.txt"));
 
 double qcm1[10][25], epsilon1[10][25], qprimcm1[10][25], thetacm1[10][25], phideg1[10][25], segmaDR1[10][25],costheta1[10][25];


   for(int j=0; j<25; j++){
     f1>>f[0];
     f1>>f[1];
     f1>>f[2];
     f1>>f[3];
     f1>>f[4];
     f1>>f[5];
     f1>>f[6];
     f1>>f[7];
     f1>>f[8];
     f1>>f[9];
     f1>>f[10];
     f1>>f[11];
     f1>>f[12];
     f1>>f[13];
     f1>>f[14];
     f1>>f[15];
     f1>>f[16];
     f1>>f[17];
     f1>>f[18];
     f1>>f[19];
     f1>>f[20];
     f1>>f[21];
     f1>>f[22];
     f1>>f[23];
     f1>>f[24];

   }

    for(j=0; j<25; j++){
    for(int i=0; i<10; i++){

     //lecture des fichiers txt
    f[j]>>qcm1[i][j]>>epsilon1[i][j]>>qprimcm1[i][j]>>thetacm1[i][j]>>phideg1[i][j]>>segmaDR1[i][j];
    fexp>>segexp[i]>>esegmaexp[i];

   //calcul de Khi2 pour chaque points de sigma
   khi[i][j]=((segexp[i][j]-segmaDR1[i][j])*(segexp[i][j]-segmaDR1[i][j]))/(esegmaexp[i][j]*esegmaexp[i][j]);

   }}



 // calcul de Khi2 somme
   for(int j=0; j<10; j++){
    for(int i=0; i<25; i++){
    khi11p[i] +=khi[i][j];
    }}

//remplissage de fechier txt "out"

 for(int i=0; i<25; i++){ fcoef<<khi11p[i]<<endl;}

}