import ode

#-------------------Enzyme initial condition setup-------------------
# Hexokinase initial condition

Hexokinase_Vmax = 
Hexokinase_Rate = 
Hexokinase_Km = 
Hexokinase_Gamma = 

# Phosphoglucoisomerase initial condition

Phosphoglucoisomerase_Vmax = 
Phosphoglucoisomerase_Rate = 
Phosphoglucoisomerase_Km = 
Phosphoglucoisomerase_Gamma = 

# PhosphoglucoisomeraseReverse initial condition

PhosphoglucoisomeraseReverse_Vmax = 
PhosphoglucoisomeraseReverse_Rate = 
PhosphoglucoisomeraseReverse_Km = 
PhosphoglucoisomeraseReverse_Gamma = 

# Phosphofructokinase initial condition

Phosphofructokinase_Vmax = 
Phosphofructokinase_Rate = 
Phosphofructokinase_Km = 
Phosphofructokinase_Gamma = 

#-------------------Reactant intial condition setup-------------------

Glucose_Vmax = 
Glucose_Gamma =
G6P_Gamma = 
F6P_Gamma = 
F16P_Gamma = 



#-------------------Concentration of reactants and enzymes setup-------------------
y = range()
 
#Enzyme side

#Reactant side
 
 
#-------------------Definition of reactants/enzymes concentration function with proper rate law-------------------

#Enzyme side


def Hexokinase(t,y):

def Phosphoglucoisomerase(t,y):

def PhosphoglucoisomeraseReverse(t,y):

def Phosphofructokinase(t,y):



#Reactant side

def 

Create F16P.py