/* 
 * 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 <iostream>
#include "BiStableSwitch31st.h"
#include <math.h>

using namespace std;

            double Volume                        =pow(10,12) ; //Total Volume of Container

void CellGrowth(double y[], double dy []){
		
		double p = y[10] ; //Stores cell number for start of increment
            	y[10] *= 1.0000288815496 ; //Adds new cells to colony at exponential rate
            	y[11] = y[10] ; 
            	y[12] = Volume - y[11] ; // VolumeOutside
            	dy[10] =  y[10] - p ; // dyCellNumber
            	dy[11] = y[11] - p ; //dyVolumeInside
            	dy[12] = - dy[10] ; //dyVolumeOutside

		}