#include <iostream>
#include <fstream>
#include <math.h>
/* 
 * This program is free software; you can redistribute it and/or modify 
 * it under the terms of the GNU General Public License as published by 
 * the Free Software Foundation; either version 2 of the License, or 
 * (at your option) any later version.
 * 
 * This program is distributed in the hope that it will be useful, but 
 * WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 
 * or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 
 * for more details.
 * 
 * You should have received a copy of the GNU General Public License along 
 * with this program; if not, write to the Free Software Foundation, Inc., 
 * 59
 */

# include "CholeraRelaySystemHeader.h"
#include <iostream>

using namespace std;

int main() {

// Initilaizes ODE initial values
double dy[var];
for (int i= 0; i < var; i++) {
            dy[i] = 0.0;
            }

ofstream datfile;
datfile.open ("CholeraTest.dat");

double y [var] = {0.0, // LuxOmRNA
                        0.0, // LuxO
                        0.0, // LuxO-P
                        0.0, // QrrsRNA
                        0.0, // HapRmRNA
                        0.0, // QrrsRNAHapRmRNA
                        0.0, // HapR
                        0.0, // HapAmRNA
                        0.0}; // HapA

            for (int j = 0; j < N; j++) {
            RungeKutta(y, dy);
            datfile << j * 0.05 << "        " <<  y[0] << "         " << y[1]<< "        " << y[2] << "        " << y[3] << "        " << y[4] << "        " << y[5] << "        " << y[6] << "        " << y[7] << "        " << y[8]  << endl;
            }

datfile.close();

return 0;

}