using System;

public class Program

{

    public static void Main()

    {

        Console.WriteLine("Bienvenido a Pokemania");

        Console.WriteLine("Esta es la prueba 3 pero hay vamos xdxd");

        Console.WriteLine();

        Console.WriteLine("En ella se igualaran los casos de el COVID-19 (coronavirus SARS-TIPO II xdxd (>_<) y el SARS-TIPO I (x_x) ");

        Console.WriteLine();

        Console.WriteLine("Para asi conocer el dia en el cual se tendra el mismo numero de infectados xdxd >:v ");

        Console.WriteLine();    

            //Formulamos la ecuación que nos proporcionara un estimado del comportamiento del número de infectados xd

            //Definimos las variables >:v

            //Ecuacion Virus SARS-TIPO II (COVID-19) ;-;

            double Slnyx=0, Slny=0, Sx=0, Sx2=0, Px=0, Plny=0, a=0, b=0;                        

            int n = 30;//Definimos la variable que guardará el número de datos con los que trabajaremos, en este caso son 30 :v

            //Definimos las variables de tiempo (X) y el numero de infectados con respecto al tiempo (Y)

            //Nota: En este punto decidi usar los datos de los Casos por Covid-19 en Mexico desde el dia 1 al dia 30 xdxd       

            double[] x = new double[30] { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30}; 

            double[] y = new double[30] { 1, 3, 4, 5, 5, 5, 5, 5, 6, 7, 7, 7, 7, 11, 15, 26, 41, 53, 82, 93, 118, 164, 203, 251, 316, 367, 405, 475, 585, 717};

        // Se calculan las sumatorias correspondientes utilizando el ciclo for :v

            for (int i = 0; i < n; i++) 

            {

                Slny = Slny + Math.Log(y[i]); //La sumatoria de Slny se obtendrá sumando los logaritmos naturales de cada valor de y

                Slnyx = Slnyx + (Math.Log(y[i]) * x[i]); //De la misma manera la suma de lny por x

                Sx = Sx + x[i]; //La sumatoria de las x

                Sx2 = Sx2 + (x[i] * x[i]); //Y la sumatoria de x^2 (x al cuadrado :v)

            }

            Plny = Slny / n;  // Se calculan los promedios >:v definiendo el valor actual entre el numero de elementos al igual que el siguiente xdxd

            Px = Sx / n;    

            b = ((n * Slnyx) - (Sx * Slny)) / ((n * Sx2) - (Sx * Sx));  // Finalmente sustituiremos los valores correspondientes para calcular a y b

            a = Math.Exp(Plny - (b * Px));

           //El procedimiento que se utilizo en la primera ecuacion se vuelve a utilizar en la segunda  :v

           //Ecuacion VIRUS SARS TIPO I (2002)      

            double Slnyx2=0, Slny2=0, Sx3=0, Sx23=0, Px2=0, Plny2=0, a2=0, b2=0;

            int n2 = 30;

            //Decidi usar los casos a partir del dia 21 al dia 50 del virus SARS TIPO I xd ya que fueron los unicos datos que encontre sobre infectados :v

            double[] x2 = new double[30] {21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50};

            double[] y2 = new double[30] {2875,3135,3232,3399,3582,3724,3998,4133,4278,4409,4542,4376,4483,4614,4865,5098,5387,5547,5735,5904,6135,6387,6598,6823,7135,7484,7668,7750,7975,8403};

​

            for (int j = 0; j < n2; j++)

            {

                Slny2 = Slny2 + Math.Log(y2[j]);

                Slnyx2 = Slnyx2 + (Math.Log(y2[j]) * x2[j]);

                Sx3 = Sx3 + x2[j];

                Sx23 = Sx23 + (x2[j] * x2[j]);

            }

            Plny2 = Slny2 / n2;

            Px2 = Sx3 / n2;

            b2 = ((n2 * Slnyx2) - (Sx3 * Slny2)) / ((n2 * Sx23) - (Sx3 * Sx3));

            a2 = Math.Exp(Plny2 - (b2 * Px2));

            // Se muestran los resultados en pantalla

            Console.WriteLine("Las ecuaciones de virus para igualar son");

            Console.WriteLine("f(COVID-19)=" + a + "e^(" + b + "x)");

            Console.WriteLine("f(SARS-1)=" + a2 + "e^(" + b2 + "x)");

            // ae^(bx) = a2e^(b2x)

            Console.WriteLine();

            double dia = (Math.Log(a2-a)) / (b-b2);// Se define la variable la cual nos indicara exactamente (creo) en que las dos enfermedades tendán la misma cantidad de infectados >:v          //Las ecuaciones se igualan a 0 = "ae^(bx) - a2e^(b2x) = 0"

            //Para esta ecuacion se realizo lo siguiente

            //1- Se factoriza el término que tenga la "b" más pequeña, en este caso seria b2 (b = 0.185 , b2 = 0.033)

            //2. Se factoriza de la siguiente manera "(e^(b2x) (ae^((b-b2)x) - a2) = 0)"            

            //3. Se igualan a cero, creo jeje.. 

            //4.Se despeja y elimina epsilion y de nuevo despejamos x ;-;

            //5. resultado xd =  "x = ln(a2-a)/(b-b2)"

            Console.WriteLine("El día aproximado en que tendrán la misma cantidad de infectados es el dia " + dia + " >:v"); // Mostramos el resultado pedido en pantalla >:v

            Console.WriteLine();

            Console.WriteLine();

​

            Console.WriteLine("Presione Enter para seguir con la segunda igualacion con el virus COVID-19 >:v");

            Console.ReadLine();

            Console.WriteLine();

            Console.WriteLine();

​

            Console.WriteLine("En esta otra comparacion se igualaran los casos de el COVID-19 (>_<) y La Influencia Tipo A H1N1 (x_x) ");

            Console.WriteLine();

            Console.WriteLine("Para asi conocer el dia en el cual se tendra el mismo numero de infectados xdxd >:v ");

            Console.WriteLine();  

            //Se realiza el mismo procedimiento asi que no añadi comentarios xdxd

            double Slnzw=0, Slnz=0, Sw=0, Sw2=0, Pw=0, Plnz=0, f=0, g=0;        

            int d = 30; 

            double[] w = new double[30] { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30}; 

            double[] z = new double[30] { 1, 3, 4, 5, 5, 5, 5, 5, 6, 7, 7, 7, 7, 11, 15, 26, 41, 53, 82, 93, 118, 164, 203, 251, 316, 367, 405, 475, 585, 717};

            for (int i = 0; i < d; i++) // Utilizaremos el ciclo for, el cual irá aumentando de uno en uno para avanzando entre los valores establecidos

            {

                Slnz = Slnz + Math.Log(z[i]); 

                Slnzw = Slnzw + (Math.Log(z[i]) * w[i]);

                Sw = Sw + w[i]; 

                Sw2 = Sw2 + (w[i] * w[i]); 

            }           

            Plnz = Slnz / d;

            Pw = Sw / d; 

            g = ((d * Slnzw) - (Sw * Slnz)) / ((d * Sw2) - (Sw * Sw));

            f = Math.Exp(Plnz - (g * Pw));

            //Ecuacion Virus Influenza Tipo A H1N1 X_X

            double Slnzw2=0, Slnz2=0, Sw3=0, Sw23=0, Pw2=0, Plnz2=0, f2=0, g2=0;

            int d2 = 27;

            //Como datos use el numero de infectados a la semana ya que de igual manera fue lo unico que encontre :"v

            double[] w2 = new double[27] {7,14,21,28,35,42,49,56,63,70,77,91,105,112,126,133,140,147,154,161,168,175,182,189,196,203,210}; 

            double[] z2 = new double[27] {25,367,2500,7520,11168,15510,21940,29690,44787,59814,89921,134503,162308,182166,209438,254206,277707,296471,318925,343298,378223,399232,414945,441661,482300,503536,525060};

​

            for (int j = 0; j < d2; j++)

            {

                Slnz2 = Slnz2 + Math.Log(z2[j]);

                Slnzw2 = Slnzw2 + (Math.Log(z2[j]) * w2[j]);

                Sw3 = Sw3 + w2[j];

                Sw23 = Sw23 + (w2[j] * w2[j]);

            }

            Plnz2 = Slnz2 / d2;

            Pw2 = Sw3 / d2;

            g2 = ((d2 * Slnzw2) - (Sw3 * Slnz2)) / ((d2 * Sw23) - (Sw3 * Sw3));

            f2 = Math.Exp(Plnz2 - (g2 * Pw2));

            Console.WriteLine("Las ecuaciones de virus para igualar son");

            Console.WriteLine("f(COVID-19)=" + f + "e^(" + g + "x)");

            Console.WriteLine("f(A H1N1)=" + f2 + "e^(" + g2 + "x)");

            Console.WriteLine();        

            //De igual manera seria b2 (g = 0.185 , g2 = 0.032                  

            double noche = (Math.Log(f2-f)) / (g-g2);

            Console.WriteLine("El día aproximado en que tendrán la misma cantidad de infectados es el dia: " + noche + " >:v");

            Console.WriteLine();

            Console.WriteLine();

            Console.WriteLine("Eso seria todo para esta parte xdxd >:v");

            Console.WriteLine("02878534");

            Console.ReadLine();

​

    }

}