//Input the intial levels of concentration of u 
#include <iostream>
#include<fstream>
#include <cmath>
using namespace std;

int main()
{
 //double u_0[8] ;
 //double u_1[8];
  double u[11][8]; // T+1=11 and # of unknowns=8;
   //intialize the time 0 concentration values of u
   //for (int i = 0 ; i < )
      //  u[0][i] = 0.0 ; 
      u[0][0] =6.25 ; // 0
     for (int i = 0 ; i < 3; i++ )
       u[0][1+i] =0.4 ; // 1, 2, 3
     for (int j = 0; j < 4 ; j ++)
      u[0][j + 4] = 0.0 ;  // 4, 5 ,6, 7   
      cout << "Intial concentrations are: "<< endl;
     for (int k =0 ; k < 8 ; ++k)
     cout << u[0][k]<< endl;  
return 0 ;
 }