SheriKabashi

#include <iostream>
#include <cpgplot.h>

const double N = 7900000; // Total population
const double dt = 1.0;

double ds_dt(double s, double i, double b) {
    return -b * s * i;
}

SheriKabashi

#include <iostream>
#include "ode_def.h"
#include <cpgplot.h>
#include <cmath>

using namespace std;

double dTdt(double b, double Tt, double Vt) {
    return -b * Tt * Vt;
}

SheriKabashi

#include <iostream>
#include <cpgplot.h>

const double N = 7900000; // Total population
const double dt = 1.0;

double ds_dt(double s, double i, double b) {
    	return -b * s * i;
}

SheriKabashi

import matplotlib.pyplot as plt  
import numpy as np  
%matplotlib inline

##Initial Conditions
t=[0]
T=[4*10**8]
I=[0.0001]
V=[3.5*10**(-1)]