#include <iostream>
using namespace std;
float X[10]={200, 500, 500, 300, 150, 0, 0, 0, 0, 0};
float Q[10];
float H[10];
float Hr[5][10];
const int dt=1;
int I;
int t;
class Vodosbor
{
public:
       Vodosbor (){tau=2; Kst=0.16;};
       ~Vodosbor() {};
       void Rez_vod ();
private:
        int tau;
        float Kst;
};
class Vodoem
{
public:
       Vodoem (){Km=0.46; F=19.2; Q1=19; Q2=21;};
       ~Vodoem() {};
       void Rez_vodoem ();
private:
        float Km;
        float F;
        int Q1;
        int Q2;
};
class Ruslo
{
public:
       Ruslo (){a=1760; dx=5000; dL=20000; Kf=0.1; Hg=1.0; Hx=0; NaUr=1;};
       ~Ruslo(){};
       void Rez_Ruslo ();
private:
        float NaUr;
        float a;
        int dx;
        int dL;
        float Kf;
        float Hg;
        float Hx;
};
int main()
{
    enum Avaria {MinH=1, SrH, MaxH};
    cout<<"Vodosbor \n";
    Vodosbor Vodosbor_kr;
    Vodosbor_kr.Rez_vod();
    cout<<"\nVodoem \n";
    Vodoem Vodoem_kr;
    Vodoem_kr.Rez_vodoem();
    cout<<"\nRuslovoi stok \n";
    Ruslo Ruslo_kr;
    Ruslo_kr.Rez_Ruslo();
for (t=1; t<10; t++)
{
    if (Hr[4][t]>MaxH)
    cout << t<< ": Failure!!!\n";
}
system ("pause");
return 0;
}
void Vodosbor :: Rez_vod () 
{
     cout << "Initial charge: ";
     cin >> Q[0];
     for (t=1; t<10; t++)
     {
     Q[t]=Q[t-1]+dt*(Kst*X[t-1]-Q[t-1])/tau;
     cout << t<< ": " << Q[t] <<"\n";
     }
}
void Vodoem :: Rez_vodoem ()
{
     cout << "Initial level: ";
     cin >> H[0];
     for (t=1; t<10; t++)
     {
     H[t]=H[t-1]+dt*((Q[t-1]+Q1-Q2)/F-Km*H[t-1]);
     cout << t<< ": " << H[t]<< " \n";
     }
}
void Ruslo :: Rez_Ruslo()
{
if (a*dt/dx<1)
{
     cout << "Initial level: ";
     cin >> NaUr;
     for (I=0; I<5; I++)
     Hr[I][0]=NaUr;
     for (t=0; t<10; t++)
     Hr[0][t]=H[t];
     for (t=1; t<10; t++)
     {
         cout << t<<": ";
         for (I=1; I<5; I++)
         {
             if (I==2)
             Hx=0.5;
             if (I==3)
             Hx=-Kf*(Hr[3][t-I]-Hg);
             if (I==4)
             Hx=0;
             Hr[I][t]= Hr[I][t-1]+dt*(Hx-a*(Hr[I][t-1]-Hr[I-1][t-1])/dx);
             cout << Hr[I][t]<<"\t";
         }
         cout <<"\n";
     }
}
else
cout << "The condition of stability is not carried out!\n";
}